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ABSTRACT 


Aerodynamics  and  control  characteristics  of  helicopters 
in  level  flight  and  steep  descents  at  low  speeds  are  analyzed. 
Single-rotor,  tandem-rotor,  and  compound  helicopters  are 
considered  and  are  specifically  represented  by  the  S-58, 

AH-56A,  and  YHC-1A,  for  each  of  which  derivatives  and  transfer 
functions  are  presented.  New  analytic  methods  are  used  to 
describe  the  flow  about  descending  rotors.  Using  these  methods, 
predictions  are  made  of  the  boundaries  of  the  vortex-ring  state 
and  of  tip  loss  factors.  Good  agreement  with  experiment  is 
obtained.  A  new  modular  stability  derivative  program, MOSTAB, 
is  described  and  is  used  to  calculate  derivatives  for  the  S-58 
and  AH-56A,  including  lateral- longitudinal  cross-coupling 
derivatives.  It  is  shown  that,  for  these  helicopters,  control 
of  flight  path  by  collective  pitch  becomes  difficult  for  low- 
speed  steep  descents  even  with  stability  augmentation,  due  to  a 
nonminimum  phase  transfer  function  and  the  associated  wrong-way 
step  response  characteristic.  The  effect  is  identified  as  an 
important  factor  in  the  observed  degradation  of  flying  qualities 
in  steep  approaches. 
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CHAPTER  I.  INTRODUCTION 


1.1  THE  SCOPE  OF  THIS  REPORT 

This  introduction  is  divided  into  two  parts.  In  the  first, 
we  describe  the  scope  of  the  report  and  indicate  how  it  advances 
the  state  of  the  art.  The  second  part  of  the  introduction  is  a 
summary  of  the  specific  contents  of  the  report,  chapter  by  chapter. 
The  reader  who  does  not  wish  to  read  the  entire  report  should  use 
Section  2  of  the  introduction  as  a  guide  to  areas  of  specific 
interest. 

Approach  and  landing  are  generally  the  most  demanding 
portions  of  a  flight,  for  both  fixed-wing  aircraft  and  helicopters. 
Proximity  to  the  ground  demands  increased  precision  in  controlling 
the  flight  path,  and  the  pilot's  task  may  be  extremely  difficult, 
particularly  in  IFR  conditions.  Pilot  comments  on  helicopters  have 
indicated  that  the  difficulties  become  more  pronounced  as  the  desired 
approach  angle  increases  (see,  e.g.,  Ref.  1).  This  suggests  that  at 
least  part  of  the  difficulty  stems  from  changes  in  the  dynamics  and 
aerodynamic  characteristics  of  the  helicopter  associated  with  change 
in  the  trim  condition  from  level  flight  to  descent.  This  is  the 
area  of  most  interest  to  the  present  report.  The  main  purpose  of 
this  report  is  to  identify  the  effects  of  descent  angle  on  the 
aerodynamics  and  dynamics  of  typical  helicopters,  and  to  interpret 
the  significance  of  these  effects  for  human  and  automatic  control. 

To  achieve  this  goal  it  was  necessary  to  go  back  to  funda¬ 
mentals.  As  will  be  described,  it  was  difficult  to  obtain  reliable 
stability  derivatives  for  small  perturbations  from  steep  descents, 
and  these  had  to  be  calculated.  Such  calculations  required  a  better 
understanding  of  the  effect  of  descent  angle  on  the  flow  about  the 
rotor  in  the  unperturbed  condition.  This  required  some  advances 
in  aerodynamic  theory.  These  included  delineation  of  the  boundaries 
of  the  vortex-ring  state,  plus  improved  analytic  models  for  the 
aerodynamic  tip  losses  of  rotors  in  level  flight  and  descent. 

Having  obtained  improved  flow  models  for  the  aerodynamics 
of  rotors  in  steep  descents,  it  was  necessary  to  calculate  stability 
derivatives  incorporating  these  models.  Stability  derivatives  are 
presented  here  for  three  specific  helicopters:  the  Sikorsky  S-58, 
the  Boeing-Vertol  YHC-1A,  and  the  Lockheed  AH-56A  Cheyenne.  These 
are  employed  to  represent  the  general  classes  of  single-rotor, 
tandem-rotor,  and  compound  helicopters.  Inferences  are  drawn 
about  the  behavior  of  each  class  from  the  dynamics  of  these  individual 
configurations . 
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The  variety  of  configurations  and  the  novel  flow  represen¬ 
tations  made  it  impossible  to  use  handbook  methods  of  computing 
derivatives,  or  even  existing  computer  programs.  The  available 
programs  were  tied  to  particular  classes  of  helicopter  and/or 
certain  flow  representations.  We  therefore  developed  a  new  computer 
program  for  calculating  derivatives.  This  program,  called  MOSTAB 
(modular  stability  derivative  program) ,  calculates  trim  and  stability 
derivatives  for  any  flying  vehicle,  described  as  an  assemblage  of 
rigid  and  flexible  rotors,  bodies,  and  fixed  lifting  surfaces  (wings, 
tails) .  The  MOSTAB  program  is  an  important  product  of  this  study. 

It  was  used  to  calculate  trim  conditions  and  stability  derivatives  for 
speeds  of  0,  20,  40,  60,  and  100  knots  for  the  S-58  and  AH-56A  in  level 
flight  and  descent.  For  the  YHC-1A  published  derivatives  were  used. 

Having  obtained  the  derivatives, it  was  noticed  that  strong 
cross-coupling  Detween  lateral  and  longitudinal  perturbations 
occurred,  as  evidenced  by  large  magnitudes  of  derivatives  such  as 
Lq,  Mp  as  compared  with  Lp,  M^.  In  many  previous  analytic  studies, 
these  cross -coupling  derivatives  were  not  available.  Such  an 
incomplete  treatment,  with  arbitrary  decoupling  of  longitudinal 
and  lateral  perturbations,  did  not  seem  appropriate  for  a  fundamental 
study,  and  so  all  the  cross-coupling  derivatives  were  retained; 
the  transfer  functions  thus  indicated  the  responses  in  six 
degrees  of  freedom  to  each  control  input. 

To  validate  the  above  transfer  functions ,  several  sets  of 
responses  to  step  control  inputs  calculated  using  MOSTAB  derivatives 
were  compared  with  flight  test  time  histories,  showing  good  agreement. 

Finally,  the  implications  of  the  transfer  functions  for 
control  of  each  helicopter  were  considered.  The  philosophy  was 
a  sort  of  "reverse  optimization";  instead  of  attempting  to  produce 
the  best  system,  by  modifying  the  characteristics  of  the  aircraft, 
we  searched  for  these  characteristics  which  would  cause  difficulties 
for  human  and  automatic  control.  Special  interest  centered  on 
adverse  characteristics  which  appeared  in  descent  but  were  absent 
from  level  flight.  It  was  found  that  some  transfer  functions  did 
exhibit  significant  differences  between  level  flight  and  descent. 

In  particular,  at  low  speeds,  control  of  the  normal  deviation 
from  the  desired  flight  path  by  collective  pitch  becomes  more 
difficult  as  the  desired  descent  angle  increases.  This  is  demonstrated 
by  the  appearance  of  a  right-half  plane  zero  in  the  appropriate 
transfer  functions. 
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The  major  advances  in  the  state  of  the  art  that  are  presented 
in  the  report  are: 

(i)  a  method  of  predicting  the  boundaries  of  the  vortex¬ 
ring  state. 

(ii)  improved  models  for  the  tip  losses  of  helicopter 
rotors  and  methods  of  calculating  their  effect  on 
derivatives,  including  cyclic  variations  of  tip 
losses. 

(iii)  the  MOSTAB  program,  which  calculates  trim  and 

stability  derivatives  of  any  helicopter  configuration. 

(iv)  presentation  of  derivatives  and  transfer  functions 
for  representative  single-rotor,  tandem-rotor,  and 
compound  helicopters,  including  lateral-longitudinal 
cross -coupling  effects. 

(v)  identification  of  some  significant  changes  in  closed- 
loop  control  characteristics  between  level  flight 
and  descents. 

1.2.  SUMMARY  OF  SPECIFIC  CONTENTS  OF  THIS  REPORT 

Chapter  II  presents  a  brief  summary  of  the  kinematics  of 
steep  approaches,  demonstrating  the  equivalence  between  descent  and 
deceleration  capability.  The  importance  of  descent/deceleration 
capability  to  Army  missions  is  explained.  A  description  is  also 
given  of  some  of  the  problems  that  have  been  observed  in  steep- 
approach  flight  tests.  This  chapter  contains  nothing  new;  it 
merely  sets  the  stage  for  the  detailed  technical  analyses  that 
follow,  by  explaining  their  relevance  to  practical  problems. 

Chapter  III  presents  a  theory  for  predicting  the  boundaries 
of  the  vortex-ring  state.  This  is  believed  to  be  the  first  published 
theory  which  predicts  the  major  features  of  these  unsteady  flow 
boundaries.  The  analysis  is  very  simple,  and  employs  momentum 
theory  and  actuator  disc  concepts.  Despite  this  simplicity,  the 
agreement  with  experiment  is  good. 

Chapter  IV  uses  flow  models  developed  in  Chapter  III  to 
calculate  the  derivatives  of  an  isolated  rotor  in  vertical  descent. 

For  purposes  of  calculating  derivatives,  it  is  customary  to  represent 
the  tip  losses  by  reducing  the  rotor  radius  from  R  to  an  effective 
radius  BR.  Standard  formulas  exist  for  B  in  hover,  but  in  descent 
these  disagree  with  experiment,  yielding  inaccurate  derivatives. 

It  was  therefore  necessary  to  produce  modified  formulas  for  B;  these 
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are  derived  in  Chapter  IV,  and  are  shown  to  give  good  agreement  with 
test  data  on  model  rotors  in  vertical  descent.  Data  were  not  available 
for  inclined  descent,  for  which  the  theory  predicts  that  the  tip  losses 
should  vary  cyclically.  This  cyclic  tip  loss  factor  was  included  in 
subsequent  calculations  of  derivatives.  In  addition,  some  cases  were 
re-run  without  tip  losses,  to  illustrate  their  importance  by  comparing 
transfer  functions  (given  in  Chapter  VI)  with  and  without  tip  losses. 

Chapter  V  presents  a  general  description  of  the  MOSTAB  modular 
stability  derivative  program.  This  briefly  summarizes  MOSTAB,  avoiding 
technical  detail  (which  is  given  in  Volume  II),  Chapter  V  also 
presents  a  discussion  of  the  accuracy  of  the  MOSTAB  derivatives.  Flight 
test  data  on  the  S-58  are  compared  with  predictions  from  MOSTAB.  The 
results  are  generally  in  good  agreement. 

The  derivatives  and  transfer  functions  for  the  YHC-1A,  AH-56A 
and  S-58  at  each  flight  condition  are  presented  in  Volumes 

II  -  IV,  Eigenvectors  and  residues  of  partial  fraction 
expansions  of  selected  transfer  functions  are  given  in  the  main  text 
of  the  report,  where  appropriate.  These  data  characterize  the  transient 
response  characteristics  of  each  helicopter  for  several  forward  speeds 
and  descent  angles. 

The  main  section  of  the  report  continues  in  Chapter  VI  with  the 
discussion  of  the  open- loop  and  closed- loop  dynamics  of  each  class  of 
helicopter,  as  revealed  by  the  above  data.  Chapter  VI  discusses  the 
tandem-rotor  configuration,  exemplified  by  the  YHC-1A.  Chapter  VII 
discusses  the  dynamics  of  the  S-58,  representing  single-rotor 
conf igurations .  Chapter  VIII  discusses  some  control  problems  of  the 
compound  configuration,  typified  by  the  AH-56A.  The  principal  con¬ 
clusions  of  the  report  are  given  in  Chapter  IX. 

Volume  II  provides  an  extensive  account  of  the  MOSTAB  program. 

It  should  be  noted  that  the  version  of  MOSTAB  used  in  this  report 
(MOSTAB-B)  does  not  include  rotor  stall  or  compressibility  characteris¬ 
tics,  as  these  were  not  significant  for  the  approach  flight  conditions 
of  interest  here.  Volume  II  describes  the  coordinate  systems  and  the 
rationale  for  finding  trim  and  derivatives,  presents  the  equations  used 
in  MOSTAB  and  explains  the  underlying  assumptions,  describes  the  rotor 
analysis  method  used  to  represent  the  first  flapping  mode,  includes  a 
listing  of  the  MOSTAB  program,  and  explains  computational  aspects  of 
MOSTAB. 
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Volume  III  presents  derivatives  and  transfer  functions  for 
the  YHC-1A  tandem-rotor  helicopter  at  forward  speeds  of  0,  60,  and 
80  knots,  and  several  descent  rates.  It  also  presents  derivatives 
for  the  S-58  at  0,  20,  40,  60,  and  100  knots,  and  four  descent  rates 
(0,  7.5,  15.0,  22.5  fps)  at  each  speed.  Transfer  functions  are  also 
tabulated  for  all  the  level  flight  and  maximum  descent  rate  conditions. 

Volume  IV  presents  similar  data  to  Volume  III  for  the  AH-56A 
and  reviews  the  literature  on  turbulence  representation  for  low 
altitudes. 

Reader's  Guide 


For  a  quick  overview  of  the  main  points  of  this  report,  read 
all  of  Chapter  II;  the  first  half  of  Chapter  V,  which  describes  the 
M0STAB  program;  the  summaries  at  the  ends  of  Chapters  III,  IV,  VI, 
VII,  and  VIII:  and  all  of  Chapter  X. 
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CHAPTER  II.  OPERATIONAL  ASPECTS  OF  STEEP  APPROACHES 


II. 1  THE  IMPORTANCE  OF  STEEP  APPROACH  CAPABILITY 

In  this  chapter  we  review  some  analytic  and  experimental 
work  on  helicopter  steep  approaches  in  order  to  pinpoint  the 
vehicle  parameters  of  prime  importance .  As  will  be  shown,  one  of 
the  most  important  parameters  is  maximum  obtainable  drag/lift 
ratio  expressed  as  a  function  of  airspeed .  This  determines  the 
helicopter's  descent  and  deceleration  capability. 

Figure  1  illustrates  the  balance  of  forces  in  a  straight- 
line  descent  in  still  air  with  acceleration  V.  (Note  that  $  is  negative 
for  deceleration  and  that  the  flight  path  angle  7  is  also  negative 
for  descent.)  The  relationship  between  y,  $,  and  the  drag/lift 
ratio  D/L  is 

0 

tan (-7)  =  —7—  + - ^  \  (l) 

y  "  L  g  cos(-7) 

The  second  term  is  merely  the  ratio  of  the  acceleration 
along  the  flight  path  to  the  component  of  gravitational  acceleration 
along  the  flight  path.  Since  $  is  negative  for  deceleration, 
increased  deceleration  at  a  given  V  and  y  requires  more  d/l,  which 
implies  increased  D  with  constant  L.  In  most  instances  the  value 
of  D/L  that  can  be  obtained  at  a  given  V  is  limited;  thus  descent 
angle  may  be  traded  off  against  deceleration,  but  the  sum  of  descent 
angle  and  deceleration  is  limited.  This  limit  is  fundamental  for 
helicopter  steep  approaches.  For  simplicity,  the  limit  is  defined 
by  considering  a  constant-speed  descent  for  which  the  maximum 
descent  angle  is 

(-7)  =  tan"1  (d/l)  „  (2) 

'  max  '  'max 

Assuming  that  the  landing  is  made  at  essentially  zero  forward 

speed,  the  terminal  phase  of  flight  must  involve  descent  and 

deceleration.  The  way  in  which  these  are  combined  can  greatly 

affect  the  effectiveness  of  the  mission.  This  is  illustrated 

by  Figure  2,  taken  from  Reference  2,  which  compares  two  alternative 

approach  profiles.  One  is  a  10- degree  straight-line  descent  with 

constant  deceleration  of  0.088  g.  The  other  approach  involves  two 

straight  segments:  the  first  at  l4 . 8  degrees  with  no  deceleration, 

followed  by  a  level  segment  with  0.264  g  deceleration.  In  both 

approaches  the  aircraft  is  flying  at  (d/l)  =  0.264  throughout. 

m&x 
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Figure  1 .  Forces  Acting  on  an  Aircraft  Flying  a  Straight-Line 
Accelerating  Descent. 


DESCENT  LEVEL  FLIGHT  TOTAL 

APPROACH  TAN  Y  V  At  TAN  T  V  At  TIME 

-  .176.088  91  — -  91  lie 

- .264  0  30  0  .264  31  61  sac 


Figure  2.  Alternative  Approach  Flight  Profiles. 
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From  Figure  2,  it  will  be  seen  that  the  total  approach 
time  is  reduced  from  91  seconds  to  6l  seconds  by  the  two-segment 
approach.  Undoubtedly,  further  savings  are  possible  by  means  of 
more  complicated  approach  profiles.  Much  work  has  been  done  on 
the  optimization  of  flight  paths  within  given  constraints  on  maximum 
D/L.  Interesting  though  this  is,  it  is  somewhat  outside  the  scope 
of  this  report.  Our  objective  is  to  determine  the  constraints 
limiting  (d/l) ma3C  for  given  helicopters.  This  forms  the  main  topic 
of  this  and  the  next  chapter. 

To  show  the  operational  benefits  of  increased  ( D/L) max  for 
a  given  approach  profile,  consider  Figure  5,  also  taken  from  Reference  2. 
(D/L)max  must  be  at  least  0.17  in  order  to  achieve  the  10-degree 
descent.  However,  raising  (D/L)max  from  0.21  to  0.35  reduces  the 
descent  time  from  3  to  1.5  minutes,  measured  from  the  3-wile 
point  to  touchdown. 

Further  benefits  of  increased  (D/L)max  include  enhanced 
capability  to  land  in  confined  areas  such  as  valleys  and  forest 
clearings.  In  addition,  a  high  ( D /L) max  permits  steep  downwind 
approaches.  This  may  be  important  for  forward  operations  or  where 
a  normal  into- the -wind  approach  may  bring  the  helicopter  uncom¬ 
fortably  close  to  the  enemy.  The  effects  of  head-  and  tailwinds 
are  illustrated  in  Figure  4. 

Using  Figure  4  to  construct  a  simple  example:  for  (D/L)  =  0.2, 
with  60  knots  airspeed  and  a  20- knot  headwind,  a  descent  angle  of 
17  degrees  can  be  achieved.  This  reduces  to  11.5  degrees  in  still 
air  and  to  8.3  degrees  for  a  20-knot  tailwind.  The  implications  for 
wind  shears  are  profound. 

The  net  conclusion  that  emerges  from  this  brief  review  is 
that  a  high  (D/L)max  is  desirable 

(i)  to  minimize  unproductive  approach  time. 

(ii)  to  permit  approaches  with  adverse  wind  directions  and 
to  cope  with  wind  shears . 

(iii)  to  permit  approaches  to  confined  areas. 

II. 2  FACTORS  LIMITING  (D/L)^ 

For  most  helicopters,  (D/L)max  is  limited  by  the  flow 
conditions  about  the  main  rotor(s).  This  is  illustrated  in  Figure  5, 
which  shows  the  permissible  regions  of  horizontal  versus  vertical 
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Figure  3*  Effect  of  Descent/Deceleration  Capability  on 
Unproductive  Mission  Time. 
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Figure  4.  Effect  of  Head-  and  Tailwinds  on  Descent  Angle. 
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VORTEX  RING  STATE 
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velocity  fora  single-rotor  helicopter.*  The  boundaries  shown  stem 
from  two  causes: 

(i}  the  vortex-ring  state 
(ii)  autorotation 

The  vortex-ring  state  is  a  region  of  unsteady  flow  which  occurs 
on  rotors  operating  at  high  (D/L)  at  low  speeds.  It  limits  the 
maximum  achievable  steady  (D/L) .  By  contrast,  the  autorotation 
boundary  is  not  associated  with  unsteady  flow;  it  occurs  because 
the  torque  on  the  rotor  is  zero.  Steeper  descents  would  be  possible 
if  a  braking  torque  could  be  applied  to  the  rotor.  Current  practice 
precludes  engine  braking,  as  it  is  usual  to  include  a  freewheel 
or  override  device  to  prevent  stoppage  of  the  rotor  following 
engine  failure. 

Autorotation  is  not  widely  used  in  IFR  situations  because 
the  rates  of  descent  are  usually  excessive,  particularly  for  the 
higher  disc  loading  helicopters.  The  principal  operational  limitation 
on  (d/l) max  therefore  is  due  to  the  vortex-ring  state.  This  is 
discussed  at  length  in  Chapter  III. 

11,3  OTHER  FACTORS  LIMITING  STEEPNESS  OF  HELICOPTER  APPROACHES 

The  foregoing  discussion  has  emphasized  the  importance  of 
(D/L) max*  However,  this  is  by  no  means  the  only  factor  limiting 
the  achievable  steepness  of  helicopter  approaches.  Let  us  define 
the  term  "nominal  flight  profile"  as  a  combination  of  airspeeds, 
descent  angles,  and  decelerations  which  is  within  the  limits 
permitted  by  (D/L)max.  In  practice,  a  nominal  flight  profile 
may  be  unflyable.  The  reasons  for  this  include  gust  response, 
displays,  handling  qualities,  loss  of  ground  reference,  etc. 

In  particular,  indications  exist  that  handling  qualities  deteriorate 
with  increased  steepness  of  the  nominal  flight  path.  For  example, 
Reference  1,  describing  tests  on  a  CH-5^C  helicopter,  notes  that, 
when  following  a  5-degree  nominal  flight  path,  the  rate  of  descent 
varied  from  100  to  700  ft/min.  At  y  =  -55  degrees  it  varied 
from  400  to  5000  ft/min,  which  was  regarded  as  unacceptable.  In 
addition  to  the  records  of  longitudinal  and  lateral  deviations 
presented  in  Reference  1,  pilot  comments  also  confirmed  the  deteri¬ 
oration  in  handling  qualities  with  increased -7.  A  similar  trend 
is  noted  in  flight  tests  on  a  HUP-1  tandem- rotor  helicopter  described 
in  References  5  and  4. 


In  Figure  5  the  horizontal  and  vertical  velocities  have  been  , 
normalized  through  division  by  =  (l/2p)V2  x  (Disc  Loading)1/2 

Figure  5  is  illustrative  only  and  should  not  be  scaled.  More 
precise  boundaries  are  given  in  Chapter  III. 
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To  form  a  complete  assessment  of  the  problems  involved 
in  flying  steep  approaches,  it  is  necessary  to  consider  both  the 
limitations  on  nominal  flight  profiles  and  also  the  limitations 
associated  with  small  perturbations  from  nominal  flight  profiles. 
Hence,  in  Chapter  III  we  analyze  the  limits  on  (D/L)max.  The 
results  are  used  to  determine  nominal  flight  profiles  for  the 
YHC-1A,  S-58,  and  AR-56 A,  and  subsequent  chapters  study  the 
behavior  of  these  helicopters  in  small  perturbations  from  these 
nominal  flight  profiles. 

Summary 

In  this  chapter  we  have  tried  to  "set  the  stage"  for  the 
detailed  analyses  that  follow  by  demonstrating  that: 

(i)  High  (D/L) is  required  to  minimize  unproductive 
mission  time,  to  permit  steep  descents  into  confined 
areas,  and  to  allow  approaches  from  any  direction, 
irrespective  of  the  wind  vector. 

(ii)  At  approach  speeds,  (D/L)max  is  limited  by  autorotation 
and  by  the  vortex-ring  state.  The  latter  limit  is 
more  serious  since  it  occurs  at  smaller  rates  of 
descent. 

(iii)  Within  the  bounds  set  by  (D/L)max, the  helicopter's 
steep  descent  capability  may  be  limited  by  some 
handling  qualities  factors,  which  appear  to  deteriorate 
with  increased  steepness  of  the  approach. 
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CHAP2SR  ZZZ.  THE  VORUX-RDrO  STATE 


zzz.i  ZHTRflDtmar 

In  performing  steep  descents  at  low  speeds,  helicopters 
encounter  an  unsteady  flow  condition  known  as  the  "vortex- ring  state". 

This  state  is  characterized  by  severe  thrust  fluctuations  and  difficulty  of 
control.  This  chapter  presents  a  simple  method  of  predicting  the 
combination  of  rate  of  descent  and  angle  of  descent  at  which  the 
vortex-ring  state  occurs.  Momentum  theory  and  actuator  disc  concepts 
are  employed;  despite  the  simplicity  of  this  approach,  the  results  are 
in  good  agreement  with  experiment. 

The  chapter  is  organized  as  follows:  Section  2  presents  an 
analytical  method  of  calculating  the  vortex- ring  state  boundary  for 
vertical  descent  using  simple  momentum  theory  modified  to  include 
certain  viscous  effects.  Section  3  extends  the  method  to  the  case 
of  nonvertical  descent.  In  Section  4  the  experimental  data  and 
published  analyses  are  reviewed  and  compared  with  the  theory  of  this 
report.  Section  5  gives  conclusions  and  recommendations  for  further 
work. 

III. 2  FLOW  MODEL  AND  BQUA TICKS  FOR  VERTICAL  DESCENT 

The  analysis  considers  both  vertical  and  inclined  descent,  but 
the  theory  is  most  easily  understood  by  considering  the  vertical  descent 
case  first. 


Consider  an  actuator  disc  in  a  uniform  stream  of  air 
rising  with  velocity  V.  This,  of  course,  corresponds  to  a  vertical 
descent  at  velocity  V  on  an  actual  helicopter.  The  flow  model  used  is 
illustrated  in  Figure  6.  It  consists  of  a  slipstream  with  uniform  flow 
at  any  cross-section,  surrounded  by  a  tube  of  vorticity.  This  tube  is 
modeled  by  a  series  6f  vortex  cores.  Thus,  near  the  rotor  outside  the 
tube  the  stream  velocity  is  V  upward,  and  inside  the  tube  the  velocity 
is  (v  -  V)  downward,  where  vis  the  induced  velocity  at  the  actuator 
disc.  The  rate  of  descent  of  the  centers  of  the  vortex  cores  is  the 
mean  of  these  velocities,  i.e.,  -  V^downward. 


The  key  assumption  of  the  analysis  is  that  the  vortex-ring 
state  will  occur  when  the  relative  velocity  of  the  vortex  cores  normal 
to  the  disc  falls  to  zero.  That  is,  when  the  rate  of  descent  is 
increased  to  the  point  where  vortex  cores  no  longer  move  away  from  the 
actuator  disc,  unsteady  flow  occurs.  The  critical  rate  of  descent 
Vcrit  'k^15  occurs  is  given  by 


V 

vcrit 


v 


2 


(3) 
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Figure  6.  Flow  Model  for  Vertical  Descent. 
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From  simple  momentum  theory  (e.g.,  Reference  5)> 


v  = 


(*) 


where 


Vt* irust/2pnR 


2 


o  A  Air  Density 

R  A  Rotor  Radius 


Combining  (5 )  and  (4 )  yields  the  following  formula  for  the  rate  of 
vertical  descent  at  which  the  vortex-ring  state  commences . 


V  = 
crit 


=  0.707  vh 


b) 


This  formula,  although  derived  by  the  simplest  possible  methods,  agrees 
well  with  experiments,  as  will  be  shown. 

An  alternative  formulation  of  Eq.  b)  can  be  obtained  by 
putting  J  =  v/nD,  =  T/ pn^D  >  which  gives 


^crit 


(6) 
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A  further  alternative  is  obtained  by  using  C^,  =  T/  pnfl^R^  and  X  =  v/ftR, 
giving  the  critical  rate  of  vertical  descent  as 


*crit 


(7) 


Tip  Lo§§  Effects 

It  is  necessary  to  refine  the  theory  slightly  to  allow  for 
nonuniform  flow  when  considering  high  disc  loading  rotors  and/or 
inclined  descents.  This  is  done  by  introducing  a  tip  loss  factor  B. 
Blade  elements  outside  of  radius  BR  are  assumed  to  produce  no  thrust. 
Thus,  near  the  rotor,  instead  of  a  uniform  induced  velocity  v  over  an 
area  nR,  we  assume  a  uniform  induced  velocity  v  over  an  area  nR2B  . 

This  leads  to  a  slight  increase  in  V  ..  for  vertical  descents,  but  the 
effect  is  small  for  typical  helicoptlr  disc  loadings.  The  tip  loss 
effect  for  vertical  descents  is  easily  obtained  as  a  special  case  of 
the  more  general  formulation  of  the  theory  for  inclined  descents, 
derived  below. 

III. 3  GENERAL  THEORY  FOR  INCLINED  DESCENT 

Figure  7  shows  the  assumed  flow  model  for  inclined  descent 
along  a  flight  path  inclined  at  an  angle  a  to  the  horizontal  (for 
vertical  descent  a  =  90  degrees).  The  velocity  of  the  vortex  cores 
normal  to  the  actuator  disc  is  assumed  to  be  the  mean  of  the  normal 
velocities  inside  and  outside  the  slipstream,  i.e.,  (v/2)  -  V  sina. 
Marked  unsteadiness  is  predicted  to  occur  when  the  rate  of  descent 
becomes  sufficiently  high  to  cause  the  velocity  to  fall  to  zero.  This 
condition  gives  the  following  general  formula : 


Vcrit  '  2  sina  ^ 

To  express  Eq.  (8)  in  a  more  convenient  form,  it  is  necessary 
to  rewrite  it  in  terms  of  horizontal  velocity  (V  cosa),  rate  of  descent 
(Vsina),  and  the  tip  loss  factor  B.  This  is  done  below  by  manipulating 
some  equations  of  momentum  theory. 
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Figure  7.  Flow  Model  for  Inclined  Descent  With  No  Parasite 
Drag,  Showing  Velocity  Components  Normal  to  the 
Disc. 
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From  momentum  theory  (e.g.,  Reference  5)  the  thrust  of  an 
actuator  disc  descending  at  an  angle  a  to  the  horizontal  can  be 
expressed  as 


T  =  2pjrR2B2v|Vj 


(9) 


2  2  2 
where  V  =  v  -  2Vv  sina  +  V 
r 


(10) 


Eqs.  (9)  and  (10)  differ  from  familiar  equations  of 
momentum  theory  only  in  that  the  rotor  radius  R  has  been  replaced  by 
an  effective  radius  BR.  Substituting  Eq.  ( 10)  in  Eq.  (9),  squaring, 
and  using  the  definition  of  vh  from  Eq.  (4)  gives 


(2*pR2  vh2)2  =  (2npR2  v)2  (v2  -  2V  v  sina  +  V2) 


(ID 


Expanding, 


v^  =  B^v^  -  2BV  (sina)  bV  + 


(12) 


1  - 


2B  Vsina 
Bv 


2  2 

B  V 

_2  2 
B  v 


(1 9) 


/ 

/Bv  \ 
\Vh  ) 


1  - 


2  B  Vsina 
Bv 


(5) 


(1*0 


This  is  similar  in  form  to  the  standard  quartic  equation  of  momentum 
theory  usually  derived  for  B  =  1j  i.e.,  quoting  from  Reference  5, 


1 _ 

2  Vsina  +  /v\2 
v  \v/ 


(15) 


Eq.  (19)  is  conveniently  graphed  with  v/vh  as  ordinate  and  Vsina/ 
as  abscissa,  for  a  fixed  a  forming  a  nondimens ional  graph  of  induced 
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velocity  versus  rate  of  descent  for  given  angle  of  descent.  (Similar 
graphs  are  given  in  References  6  and  7  for  B=1 . )  The  graphs  can  be 
formed  from  Eq.  (l4)  for  a  fixed  B  using  Bv/v,  as  ordinate  and 
B  Vsina/v  as  abscissa.  Figure  8  illustrates  n  these  graphs  for 
several  values  of  a,  from  5  to  90  degrees. 

To  establish  the  boundary  of  the  vortex-ring  state,  for 
a  given  B  and  a,  simply  include  in  Figure  8  the  critical  condition  of 
Eq.  (8).  Expressed  in  terms  of  the  variables  of  Figure  8,  Eq.  (8) 
becomes 


2  B  Vsina 


(16) 


This  describes  a  straight  line  of  slope  2,  as  shown  on  Figure  8.  The 
intersection  of  this  line  with  the  remaining  graphs  denotes  the 
critical  value  of  rate  of  descent  for  a  given  B  and  a.  It  is  convenient 
to  plot  the  results  in  terms  of  nondimens ional  horizontal  velocity, 
Vcosa/vh,and  nondimens  ional  rate  of  descent,  Vsina/v^,  for  a  given  B  as 
shown  on  Figures  9,  10,  and  11.  In  these  figures,  Eq.  (l6)  corresponds  to 
the  line  marked  "  zero  parasite  drag"  for  which  case  a  =  -  7,  as 
explained  below. 


Effect  of  Parasite  Drag 

The  parasite  drag  of  the  rotor  and  of  other  components  has 
not  been  included  up  to  this  point.  This  drag  causes  the  boundary  of 
the  vortex-ring  state  to  occur  at  a  steeper  angle  of  descent  for  a 
given  airspeed.  The  effects  of  parasite  drag  can  be  included  in  the 
analysis  as  indicated  below. 

Figure  12  shows  the  forces  acting  on  a  helicopter  in  a 
steady  descent;  the  aerodynamic  forces  are  assumed  to  consist  only  of 
drag  acting  parallel  to  the  flight  path  and  rotor  thrust  normal  to 
the  actuator  disc. 

From  Figure  12,  resolving  parallel  and  normal  to  the  flight 
path,  noting  that  y  is  negative  for  descent, 

T  sina  +  D  +  W  sin  7=0 

T  cosa  -  W  cos  7=0 


(17) 

(18) 
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RATE  OF  DESCENT  PARAMETER,  B  VsiNoc/i/h 


Figure  8.  Induced  Velocity  Versus  Rate  of  Descent  From 
Momentum  Theory. 
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INDUCED  VELOCITY  PARAMETER,  Bi//i/h 


HORIZONTAL  VELOCITY  PARAMETER,  VCOS(-7)/i/h 


“/I/  (X  -)NISA  *y313WVdVd  1N30S30  JO  31Vd 
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Figure  9.  Predicted  Upper  Boundary  of  Vortex-Ring  State  With 
Zero  Tip  Loss  Factor,  B  =  1 .0  . 


W-) NISA  *a313IAIVaVd  1N30S3Q  30  31Vd 
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HORIZONTAL  VELOCITY  PARAMETER  ,  VCOS  7/ 


NISA  <a313INVaVd  1N3DS3CI  30  31Vd 


2k 


Figure  1 1 .  Predicted  Upper  Boundary  of  Vortex-Ring  State 
With  Tip  Loss  Factor  B  =0.8. 


ACTUATOR  DISC 


Figure  1 2.  Balance  of  Forces  in  Descending  Flight. 
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These  equations  yield  the  required  relation  between  y  and  a  as: 


tan  a 


-(pAO 

cos  y 


tan  y 


09) 


The  drag/weight  ratio  is  conveniently  expressed  in  terms  of  a  "drag 
area"  Af  defined  as 


Af  =  D/^pV2  (20) 

2  2 

Thus,  using  the  definition  of  (v,  =  T/2p«R  ),  and  assuming  thrust 

at  hover  equal  to  weight,  Eq.  (19)  becomes 


tan  a  =  - 


-  tan  y 


(21) 


To  modify  Figures  9#  10,  and  11  to  show  the  effect  of  profile 
drag,  we  require  the  relationship  between  y  and  Af/UnR d  for  a  given  V 

and  a.  This  is  gbtained  by  manipulating  Eq.  (21)  to  give  y  in  terms  of 
a,  V,  and  Af/4rtR  ,  as  follows.  Multiply  (21)  by  C0S7  and  cosa : 


sin  a  cos  y  - 


sin  (a  +  y)  - 


7  = 


(22) 

(23) 


(2*0 
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Eq.  (24)  has  been  used  to  show  the  effect  of  parasite  drag  in 
Figures  9,  10,  and  11.  For  a  fixed  value  of  the  drag  parameter 

p 

Af/4nR  =  0.05,  the  flight  path  angle  at  any  V  changes  from  a  to  y 

as  determined  from  Eq.  (22).  For  zero  parasite  drag,  a  =  -  y. 

Note  that  the  inclusion  of  parasite  drag  yields  the  characteristic 
shape  of  the  vortex-ring  boundary,  first  rising  and  then  falling, 
as  horizontal  speed  is  increased. 

2 

The  drag  parameter  Af/4jtR  was  deliberately  chosen  to  be 

rather  high,  so  that  most  helicopters  display  vortex-ring  boundaries 
falling  between  the  "zero  drag"  and  "with  drag"  curves  of  Figures 
9,  10,  and  11.  To  appreciate  the  physical  significance  of  the  parameter 
Af/4nR2,  note  that  a  value  of  0.05  corresponds  to  a  parasite  drag  area 

of  20  percent  of  the  disc  area.  This  is  much  larger  than  would  be 
expected  in  level  flight  for  most  helicopters.  However,  in  descent, 
the  fuselage  may  be  operating  at  a  large  angle  of  attack,  causing  a 
corresponding  increase  in  drag  area.  When  comparing  Figures  9,  10 
and  11  with  exp  rimental  results  on  wind  tunnel  models,  note  that  the 
parasite  drag  of  the  fuselage,  tail,  and  other  components  will  be 
disproportionately  large  at  low  Reynolds  numbers. 

To  predict  vortex-ring  boundaries  of  specific  configurations 
using  Figures  9,  10  and  11  requires  the  value  of  the  tip  loss  factor,  B. 
Well-known  formulas  exist  for  B  in  hover J  e.g.,  Payne  (Reference  8 ) 
quotes  a  formula  due  to  Sis singh, 


Payne  notes  that  Sissingh's  formula  is  in  reasonable  agreement  with 
tests  on  untwisted,  untapered  blades  but  yields  values  of  B  which  are 
too  small  for  uniformly  loaded  blades.  For  the  latter,  Payne  suggests 
an  alternative  formula,  which  at  hover  reduces  to 


B  =  1 


(26) 
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For  high-disc-loading  VTOL  rotors,  both  the  above  formulas 
appear  to  give  values  of  B  which  are  too  small.  An  alternative  semi- 
empirical  expression  is  derived  in  Chapter  IV  for  rotors  of  ideal 
twist  and  is  given  below.  This  is  more  complicated  than  Eqs.  (29)  or 
(26) but  gives  good  accuracy  for  vertical  descent  as  well  as  hover. 


B  =  1  -  £  ^  T  (9t  +  V/BnR)  (27) 


This  gives  a  cubic  equation  for  B,  for  a  given  rotor  geometry,  pitch 
setting,  and  rate  of  descent.  Empirically  it  is  suggested  that  Eq.  (27) 
can  be  used  for  inclined  as  well  as  vertical  descent,  replacing  V  by 
Vsina. 

Ill  .4  COMBARISCW  WITH  EXPBRBflSET 

Although  investigations  of  the  behavior  of  rotors  in  the 
vortex-ring  state  have  been  made  since  the  1920's,  only  the  more  recent 
references  contain  data  useful  for  comparison  with  the  theory  of  this 
report.  Most  of  the  earlier  references  (e.g.,  Reference  9 )  are  confined 
to  measurement  and  empirical  prediction  of  the  mean  thrust  and  through- 
flow  velocity.  Reference  10  summarizes  various  presentations  of  parameters 
such  as  l/f  and  l/F.  These  parameters  are  derived  from  measurements  of 
averaged  quantities,  so  these  data  are  unsuitable  for  determining  the 
velocity  and  angle  of  descent  at  which  noticeable  unsteadiness  occurs. 
References  11,  12,  19,  1^,  and  19  do  not  include  any  unsteady  force 
measurements,  but  they  do  present  some  flow  visualizations,  which  quali¬ 
tatively  indicate  a  "region  of  roughness."  These  data  are  discussed 
later.  For  purposes  of  comparison  with  theory,  it  is  unfortunate  that  the 
available  test  data  are  for  blades  with  non- ideal  twist;  e.g.,  References 
11,  12,  19,  14,  and  19  employ  an  untwisted  rotor.  The  resulting  induced 
velocity  distribution  is  at  variance  with  the  uniform  distribution 
assumed  in  the  present  theory.  However,  even  full-scale  helicopters 
do  not  have  ideal  twist,  so  this  compromise  must  be  accepted  in 
order  to  compare  theory  and  experiment. 

Vertical  Pea cent 

Azuma  (Reference  16)  presents  data  on  thrust  and  inflow 
fluctuations  for  a  9-bladed  rotor,  with  a  diameter  of  1100  mm,  a  solidity 
of  0.0979,  an  NACA  0012  section  and  -8°  twist  from  root  to  tin,  free  to 
flap  and  lag.  The  rate- of-des cent  parameter  was  expressed  as  Vjl/vfi, 
because  the  nominal  rotational  speed  of  1,000  rpm  was  not  maintained 
precisely. 
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Lacking  data  on  v/ftR,  we  cannot  make  allowance  for  the  tip 
loss  factor.  However,  Eq.  (5),  vcrit  =  0.707  \^>  gives  good  agreement 

with  the  observed  onset  of  significant  unsteady  thrust  and  inflow. 

This  is  demonstrated  in_Figure  13.  Eq.  (5)  corresponds  to  r.ra.s. 
thrust  fluctuations  AT/T  ranging  from  kfo  to  1  k’jo  depending  on  the  blade 
pitch  setting.  (The  thrust  fluctuation,  AT,  is  defined  as  T  -  T, 
where  T  is  the  mean  thrustj  thus  the  r.m.s.  peak  to  peak  fluctuation 
is  8je  to  2&Jo  of  the  mean  thrust.)  Note  that  the  graph  marked  0  =  8 

degrees  corresponds  to  zero  blade  angle  at  the  tip.  As  predicted  by 
Eq.  (3)and  Figures  9,  10,  and  11,  the  critical  rate  of  descent  normali¬ 
zed  with  respect  to  v.  and  fL  increases  as  the  tip  blade  angle  is 
increased  (by  increasing  0  ).  The  experimental  data  of  Reference  16 
indicate  that  the  mean  thrust  in  descent  variesj  for  0  =10  degrees  at 
V  =  0.7  v,fl/nh,  the  mean  thrust  falls  to  8o/o  of  the  holering  thrust. 

To  convert  the  data  to  the  constant  mean  thrust  assumed  by  the  theory 
(and  appropriate  to  steady  descents),  the  value  of  V  ..  must  be 
increased  above  that  observed  experimentally.  Withou£1data  on  CT, 
this  cannot  be  done  precisely,  but  the  trend  will  certainly  be  to  move 
the  peaks  of  the  thrust  fluctuation  curves  to  the  right,  thus  increas¬ 
ing  agreement  with  the  predicted  boundary  for  the  onset  of  unsteady 
flow. 


Washizu  (Reference  7)  describes  tests  on  a  similar  rotor 
in  vertical  and  inclined  descent.  Figure  14  shows  envelopes  of  CL, 
versus  the  rate  of  descent  parameter  X,  for  vertical  descent.  For 
two  of  the  tests,  the  blade  pitch  setting  is  very  low  and  the  tips  are 
actually  lifting  downward  in  hover.  This  invalidates  comparisons  with 
our  theory,  which  assumes  uniform  inflow  except  at  the  tips, 
corresponding  to  the  rather  high  pitch  setting  appropriate  to  maximum 
static  thrust/power.  For  this  reason,  the  data  taken  at  Q,rf9  =  4.5 
degrees  are  also  unsuitable  for  comparison.  However,  using  Eq.  (7) 
and  B  =  0.95>  the  case  with  0^  =8.0  degrees  gives  ^-crit  =  0.033. 

This  corresponds  to  a  mean- to- peak  thrust  fluctuation  of  12 Jo  of  the 
mean  thrust. 


Inolined  Pea  pent 

Washizu  et  al  presents  experimental  data  on  thrust  and 
induced  velocity  fluctuations  in  the  vortex-ring  state  for  single¬ 
rotor  configurations  (Reference  7)  and  tandem  configurations  (Reference 
17).  Contours  of  percentage  thrust  fluctuation  are  shown  on  the  hori¬ 
zontal  velocity/rate-of-descent  plane  j  however,  the  data  points  are  not 
indicated.  Hence  results  for  low  pitch  angles  cannot  be  removed  to 
enable  the  remaining  data  to  be  compared  with  the  theory.  Washizu 
et  al  gives  more  complete  data  on  induced  velocity  fluctuations , 
averaged  over  the  disc  based  on  measurements  of  thrust  and  torque 
fluctuations,  and  these  can  be  compared  with  theory.  Figures  1 3  and  16 
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Figure  13  •  Measured  Thrust  Fluctuation  for  a  Model 
Helicopter  Rotor  in  Vertical  Descent. 
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THRUST  COEFFICIENT 


NONDIMENSIONAL  RATE  OF  DESCENT,  X 


Figure  14.  Measured  Thrust  Coefficient  for  a  Model 
Helicopter  Rotor  in  Vertical  Descent. 
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illustrate  the  data  of  Reference  7  for  rotor  with  9  _  =  7«5  degrees 

and  a  =  50  and  70  degrees.  Note  that  the  lower  boundary  of  the  induced 
velocity  is  close  to  the  momentum  theory  solutions  and  that  the  onset 
of  marked  unsteadiness  corresponds  to  the  critical  condition  predicted 
by  Figure  8.  This  critical  condition  corresponds  to  the  shaded  lines 
on  Figures  15  and  1 6. 

Yaggy  and  Mort  (Reference  18)  describe  tests  on  a 
flapping  propeller  with  the  shaft  axis  inclined  at  180,  165,  150,  135, 
and  120  degrees  to  the  freestream.  This  corresponds  to  descent  angles 
of  90,  75,  60 ,  45,  and  30  degrees.*  Contours  of  the  thrust  fluctua¬ 
tion  AT_(measured  from  the  mean  thrust  to  the  peak)_divided  by  the  mean 
thrust  T  are  shown  for  various  mean  disc  loadings  T/nR  ,  in  Figure  7 
of  Reference  18.  Figure  17  illustrates  the  results  for  T/rtR 2  =  16 
lb/ft2.  The  vortex- ring_boundary  has  also  been  calculated  assuming 
B  =  0.86  and  a  mean  v^,  A  (T/2p7tR2)  V2,  The  boundary  corresponds 

to  a  AT/f  of  between  1 5$  and  25 fj>.  A  more  accurate  boundary  can  be 
calculated  if  v  is  based  on  the  lower  value  of  thrust  T  -  AT.  This  is 
in  accord  with  Figures  15  and  16,  which  indicate  that  the  minimum 
induced  velocity  during  the  vortex-ring  state  is  approximately  equal 
to  the  velocity  predicted  by  momentum  theory.  To  follow  this  sugges¬ 
tion,  substitute  the  critical  condition  v  =  2Vsina  in  Eqs.  (9)  and 
(10).  This  gives 


T  -  AT  =  UnR^Cv  sin  a)2  /sin  a 


(28) 


The  critical  rate  of  descent  can  be  calculated  itera¬ 
tively  by  first  putting  AT  =  0.  This  gives  the  solid- line_boundary 
on  Figure  17.  Since  in  this  case  we  have  knowledge  of  AT/T,  we  can 
refine  the  solution  by  reducing  the  critical  rate  of  descent  in  the 
ratio  (l  -  AT/t)1/2.  The  resulting  boundary  is  also  shown  on 
Figure  17  as  a  broken  line.  Without  prior  knowledge  of  AT/T  (or  an 
estimate),  it  would  be  possible  to  calculate  only  the  solid-line 
boundary. 


Flight  Tests 

Published  data  on  flight  characteristics  of  helicopters 
in  steep  descents  are  insufficiently  complete  to  permit  mapping  the 
boundaries  of  the  vortex-ring  region.  However,  some  related  data 
points  can  be  obtained  as  listed  below. 


The  descent  angles  of  75,  60,  ^5,  and  30  degrees  are  necessarily 
approximate,  because  Ref,  18  does  not  give  the  inclination  of  the 
thrust  from  the  shaft  axis. 
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Figure  15.  Fluctuations  of  Mean  Induced  Velocity,  Deduced 

From  Thrust  and  Torque  Measurements,  a  =  50  Degrees. 
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Figure  16.  Fluctuations  of  Mean  Induced  Velocity,  Deduced  From 
Thrust  and  Torque  Measurements,  a  =  JO  Degrees. 
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Jteftrtftet  19  not**  that  for  •  WtaiUnd  Whirlwind  helicopter, 
iht  vortex- ring  «tete  mm  wall  established  at  1*00  ft/mln  vortical 
deaoonti  oorroa  pond  Inc  to  V/vh  *  0.6t.  Toata  on  an  early  helicopter, 
the  Sikorsky  IMS,  described  Xn  Reference  abowed  onaot  of  roughness 
in  vortical  daaoent  at  $00  f.p.e.  and  "shuddering"  coanonoinc  at 
900  f.p.e,  eorrocpondlnc  to  O.TOTv^.  In  the  discussion  of  Reference  «?i, 

Honnott  aUtca  that  for  the  IMP,  the  eritioal  rate  of  doaoent  decreases 
as  forward  speed  is  increased,  almost  reaching  level  flight  at  Co  e.p.h. 


Die  above  results  are  in  aeoord  with  the  theory,  but  it 
oust  be  roe-bored  that  factors  other  than  the  behavior  of  the  ro*or 
Itself  in  the  vortex-ring  state  nay  influence  helicopter  flying  dualities 
in  desoents.  Die  unsteady  flow  nay  induce  largs  pitching  moments  on  the 
ruse  lege.  Bsferenoe  22  discusses  this  phononenon,  noting  that  the  else 
and  cross-section  of  the  fuselage  tail  cone  have  an  important  effect. 

For  example,  the  IMb  had  particularly  bad  characteristics  in  the  vortex¬ 
ring  state  because  of  its  thick  slab-sided  rear  fuselage.  Quite  apart 
from  considerations  of  unsteady  flow,  changes  in  derivatives  nay  affect 
flying  qualitias.  Ro  fa  ranee  19  notes  that  the  magnitude  of  the  plunge 
damping  derivative,  Z  ,  decreases  markedly  between  hover  and  the  vortex- 
ring  state.  Dlls  treXd  is  confirmed  by  the  data  of  Rsfarenca  IB.  if 
I Z  j  become^  sufficiently  small,  the  helicopter' ■  flying  qualitlee  will 
do  degraded. 


D>t  T^ftT  Boundary  of  tbs  Vortex- Rins  State 

After  entering  the  vortex-ring  state,  if  the  rate  of  descent 
is  increased,  eventually  steady  flow  is  re-established.  The  horizontal 
speed  and  rate  of  descent  at  which  this  occurs  define  the  lower  boundary 
of  the  vortex-ring  state.  The  practical  importance  of  the  vortex-ring 
state  depends,  in  part,  on  the  distance  between  the  lower  and  upper 
boundaries  on  the  V  cos  7,  V  sin  7  hodograph  plane.  If  large  areas  of 
the  plane  are  associated  with  unsteady  flow,  the  operational  utility 
of  the  helicopter  will  be  correspondingly  reduced. 

A  critical  condition  for  the  lower  boundary  will  now  be 
derived.  As  before,  we  postulate  that  the  unsteady  flow  is  associated 
with  a  breakdown  in  the  protective  sheath  of  vorticity  surrounding  the 
slipstream.  However,  for  the  lower  boundary  the  "pileup"of  vorticity 
occurs  some  distance  above  the  rotor,  as  shown  in  the  flow  model  of 
Figure  18. 


See  Ref.  3^  for  an  explanation  of  the  influence  of  Z  on  handling 
qualities.  w 
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Figure  18.  Flow  Model  for  Lower  Boundary  of  Vortex-Ring 
State  in  Vertical  Descent. 
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Using  mcmatm  theory  and  ttoa  flow  nodal  of  flgura  the 
orltloal  eondltion  la  darlvad  u  follows.  Ins  Ida  the  •llpttrwam  the 
induuad  rtloclty  la  tv,  whore  the  value  of  k  varies  fro*  1 .0  at  the 
rotor  disc  to  2.0  at  Infinity.  The  rate  of  asoent  of  the  vortloee  le 
Vsli*»(kv/2)«  Zn  determlnlns  the  upper  boundary,  conditions  at  the 
disc  are  relevant,  and  so  a  value  of  k  *  1 .0  ms  used,  giving  the 
critical  condition  as  Bq.  (}),  or  nore  gene  rail/,  Bq.  ( \t>) .  for  the 
lower  boundary  the  cholos  of  k  Is  not  so  clear-cut.  At  Infinity  k  »  2, 
but  this  value  Is  Inappropriate  because  an  aeewulatlcn  of  vortlelty 
at  Infinity  would  not  oauee  unsteady  flow  to  be  reflected  back  to  the 
rotor.  Therefore,  the  value  of  k  wet  be  less  than  2.0.  Tentatively, 
l.b  <  k  <  1.6  Is  suggested.  Cross-plotting  lines  of  flops  2/k  on 
figure  P  with  k  •  FTel  and  k  ■  1  &  gives  a  lower  boundary  as  indicated 
on  figure  19* 


It  Is  interesting  to  note  that  the  "region  of  roughness " 
reported  by  Drees  in  Referenoe  H  on  the  basis  of  flow  visualisation 
agrees  approxieately  with  the  predicted  boundaries  at  low  speeds. 

The  lower  boundary  is  not  as  sharply  defined  as  the  upper 
boundary,  in  that  the  development  of  unsteady  flow  is  less  sudden 
(see,  e.g.,  Figure  13).  One  would  expect  this,  because  the  ac emulation 
of  vorticity  and  consequent  breakdown  of  the  slipstream  are  occurring 
farther  away  from  the  rotor.  One  would  also  expect  meuentun  theory  to 
become  less  accurate  far  predicting  the  lower  boundary  at  higher  forward 
speed,  because  distortion  of  the  slipstream,  rolling  up  of  vortices, 
and  other  viscous  phenomena  not  included  in  the  theory  become  more 
Important  as  the  critical  region  moves  away  from  the  rotor  with  increase 
of  k. 


Tsndem  Rotors 

Reference  17  describes  experiments  on  tandem  rotors  in 
vertical  and  inclined  descent.  Data  an  fluctuating  inflow  (similar 
to  Figures  15  and  16)  are  presented,  normalized  with  respect  to 

alternative  parameters  vh  and  where  vh  =  T/2pnR2,  =  T/2pSe, 

e  e 

where  T  =  thrust  per  rotor,  end  S  is  one-half  the  projected  area  of 

the  tandem  rotors.  The  overlap  bf tween  the  rotors  is  small,  so  for 

the  configuration  of  Reference  17, S  should  be  defined  in  this 

6 

way  for  momentum  theory  calculations;  Reference  6,  p.  312,  suggests  a 
larger  area.  This  is  the  area  formed  by  two  semicircles  of  radius 
equal  to  the  rotor  radius,  with  centers  spaced  a  distance  equal  to  the 
hub  separation  of  the  actual  configuration,  plus  the  area  enclosed  by 
straight  lines  joining  the  semicircles. 
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The  experlnental  data  of  Reference  *7  Indicate  that  severo 
Induced  floe  fluctuations  occur  at  a  •  lightly  lover  value  of  — 

vh 

than  that  for  the  single- rotor  configuration.  However ,  the  perlodlolty 
noted  In  Reference 7  1>  hardly  observable.  Because  of  the  Mutual 
lnduotlon  effect  of  tandea  rotors,  it  la  expected  that  the  theory 
proposed  here  will  require  Modification  for  such  configurations* 


nx  *9  tuwmnr  /mo  nxco-CMPATiom 

i<— t* 

1 .  ft*  bowdarlee  of  the  vortex-ring  stoto  oar  ba  pre- 
dlotad  by  mmmtm  theory  aa  daaerlbad  in  this  chapter.  Tba  predicted 
uppar  bowdary  oorreeponde  to  neeeured  nean-to-peek  thruat  fluctuations 
of  approKiaataly  of  the  mean  throat. 

2.  Zb  vertloal  da# oast,  tba  uppar  boundary  la  approximately 
0*707  tinee  tba  naan  induoad  veloolty  at  bovor.  For  inclined  descent* 
tba  uppar  boundary  la  vary  aanaltlva  to  paraalta  drag,  and  for  hellooptere 
with  large  paraalta  dreg, the  vortex-ring  atata  nay  vanlab  for  daaoant 
an&Ies  a  bailout  r  than  SO  degrees. 


1  *  Znprovad  ballooptar  ataap  daaoant  capability 
requires  a  nora  oonplata  undo ra tending  of  tba  vortex-rlng 
atata.  fba  "uppar”  (k  ■  1)  boundary  la  of  principal  interact  for 
ballooptara  ;  but  for  8T0L  aircraft  with  reverse-pitch  propa liars  daaoant 
angle  la  United  by  the  "loner"  (k  •  1.9)  boundary,  as  shown  in  Figure 
20*  Further  knowledge  of  tba  vortex- ring  state  Is  thus  of  lnportanoe 
in  obtaining  steeps  slow  flight  of  props ller-drivsn  8T0L  aircraft  as  well 
as  helicopters. 


2.  lbs  present  theory  is  United  by  the  assumptions  of 
nomen  tun  theory,  e.g.,  the  rather  arbitrary  Olauerb-Kussner  hypothesis 
regarding  the  "area  of  capture"  and  the  assumption  of  uniform  induced 
velocity  across  the  disc.  Greater  accuracy  could  be  obtained  by  using 
more  refined  theories.  A  further  area  for  exploration  by  such  theories 
Is  the  nature  of  the  flow  within  the  vortex-ring  state.  For  example,  it 
appears  that  the  accumulation  of  vorticity  either  at  or  above  the  disc 
increases  the  induced  velocity  at  the  disc,  thus  removing  the  accumulated 
vorticity  and  returning  the  flow  to  that  predicted  by  momentum  theory. 

In  turn,  this  causes  the  vorticity  to  re- accumulate,  causing  a  cyclic 
oscillation  characteristic  of  the  vortex-ring  state.  Mathematical 
formulation  of  this  physical  model  probably  lies  beyond  the  grasp  of 
momentum  theory. 


3 .  Although  the  upper  boundary  is 
slight  unsteadiness  does  occur  at  smaller  rates 


fairly  sharply  defined, 
of  descent,  and  even 
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Figure  20.  Effect  of  Reverse- Thrust  Propellers 
on  Descent  Angle  of  a  STOL  Aircraft. 
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•t  hover  according  to  the  data  of  Reference  A  nor*  sophisticated 
theory  is  required  to  analyte  this  phenomenon. 

k.  Hors  experimental  date  on  the  vortex- ring  stats  ere 
required*  The  present  analysis  oould  serve  as  a  quids  to  the  test 
conditions  of  primary  interest. 

5*  the  implications  of  the  vortex- ring  state  on  headline 
qualities  and  operational  maneuvers  require  laves  ties  ti  on.  Suitable 
automatic  control  of  collective  pitch  may  smooth  the  thrust  fluctuations 
sufficiently  to  permit  comfortable,  controlled  flight  at  all  desoent 
speeds  and  angles. 
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CHAPTIK  IV.  THMET  Of  AM  IfOLATCP  DC  SO*  DING  MOTOR 


Standard  msthods  for  calculating  bcliooptcr  stability 
derivatives  are  baaed  on  momenta*  and  blade-element  theories  (e.g., 
Mofbresoes  2k,  25) .  At  will  bo  shown,  thoso  procedures  do  not 
give  good  agreement  with  experimental  data  on  rotoro  In  daenendlng 
flight*  This  ehaptor  therefore  proton t«  an  enalyele  of  rotor*  In 
vo rt leal  doaoont  using  conventional  momentum  and  blade- element 
assumptions,  plus  sons  additional  rofinoaonts.  These  lnoiudo: 

(i)  non  uni  font  downvasb  assoelatod  with  rotor  tip 

vortloos 

(11)  ttao  loss  of  lifting  effectiveness  of  tho  tips 

duo  to  thrse-dlmcnslonal  flow  (l.o.,  "spillage") 

With  thoso  factors  included,  tho  theory  agrees  well  with  experiment. 

The  basic  flow  model  is  shown  in  figure  6*  The  tip 
vortices  are  modeled  by  vortex  cores  with  boundary  conditions 
matched  inside  and  outside  the  slipstream.  This,  of  course.  Is  a 
highly  simplified  model  of  the  actual  flow,  neglecting  swirl,  non- 
uniform  downvasb  due  to  non-ideal  blade  twist,  unsteady  effects,  etc. 
However,  this  model  gives  results  which  are  sufficiently  accurate 
for  stability  and  control  calculations,  predicting  derivatives  to 
within  +  20J&  accuracy. 

The  chapter  is  organized  as  follows.  In  Section  2,  mass 
flow,  momentum,  and  energy  relationships  are  derived  for  the  assumed 
flow  model,  with  the  ratio  of  the  parallel  sides  of  the  trapezoidal 
distribution  left  variable,  These  basic  relationships  relate 
conditions  at  the  rotor  disc  to  those  in  the  ultimate  wake.  Because 
of  the  nonuni  form  flow  at  the  outer  edge  of  the  slipstream,  the  rela¬ 
tionships  are  more  complicated  than  for  the  standard  model  where 
the  flow  is  assumed  to  be  uniform  across  the  slipstream.  Section  } 
expresses  rotor  thrust  in  terms  of  blade-element  angles  of  attack, 
taking  into  consideration  the  trapezoidal  induced  velocity  profile. 

In  Section  4  the  vortex  core  model  boundary  conditions  are  established 

in  order  to  fix  the  geometry  of  the  trapezoidal  distribution  both  at 
the  disc  and  in  the  ultimate  wake. 
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Motion  *j  dismisses  tbo  loot  of  lifting  effectiveness  of 
bind#  eleasats  no  or  tbo  tips  duo  to  three-dlaraslonal  flow.  This 
it  usually  aooounted  for  by  a  "tip  loss  factor".  l.i.,  do  fining  an 
offbotlvo  rotor  radius  BR  where  R  •  1  -  yzcl/b.  ibis  faotor  gives 
reasonable  aoouraoy  at  hover,  but  it  is  not  apparent  bow  B  should  bo 
Modified  for  nonboverlng  oonditiono*  Motion  5  thoroPoro  pro  son  ts 
an  alternative  derivation  of  B  btsod  on  oons  ldo rations  of  tbs  loading 
at  tbs  rotor  tips  ratbsr  than  tbs  average  loading  par  blade.  this 
pa ral to  B  to  be  oaloulatod  In  a  logloal  fashion  for  nonbovorlng 
oonditiono.  it  does,  however,  load  to  a  oyollo  variation  of  B 
which  eausoo  oons  additional  oonplioation  in  oaleulating  derivatives. 
Tbo  proooduro  adoptad  for  Incorporating  a  oyollo  B  in  tbo  rotor 
binds  tins  history  oaloulatlono  lo  too  dstailsd  for  dloouoolon 
in  tbs  aain  taxt  of  tbs  report  and  Is  tbaraforo  dloouosod  In 
Appendix  II,  In  conjunction  with  tbo  dsooription  of  tbs  subroutine 
8WKBP  which  Integra  too  tbs  blade-alsaant  loads  radially  and  with 
rosnoot  to  tins. 

Beotian  6  shows  bow  tbs  tip  loss  factor  B  can  be  oonblnsd 
with  tbs  pnsoedlng  analysis  to  ooaputa  thrust  versus  rats  of 
descent  for  a  given  rotor.  Bastion  7  oospores  tbsoratioal 
predictions  with  wind- tunnel  data  on  a  descending  rotor.  Because 
no  satisfactory  data  were  available  for  inclined  desoent,  it  was 
not  possible  to  validate  tbs  assuMd  oyollo  ocaponent  of  tip 
losses,  but  good  agreement  is  bbtainsd  with  data  fOr  vertical 
desoent  for  which  only  noncyclio  tip  losses  are  predicted. 


IV. 2.  BASIC  RELATIONSHIPS  FOR  MASS  FLOW,  MOMENTUM,  AND  POWER 

Axial  Velocity  Distribution 

A  trapezoidal  velocity  distribution  such  as  that  shown  in 
Figure  6  can  be  described  as  follows: 

|  v  -  V  0  £  x  £  1-6 

v(x)  =  v  ^  _  x)_v  1-8  *  x  &  1 
l  8 


(29) 


This  distribution  can  be  nondimensionalized  by  dividing 
through  by  QB.  The  result  is 


a(x)  = 


X  -  A 
\  (l-x)-A 


0  £  x  £  1-6 
1-8  <  x  <  1 


(30) 


The  above  equations  apply  near  the  rotor  disc.  Similar  expressions 
hold  in  the  wake,  except  that  a  is  used  in  lieu  of  a.  Of  course, 

X  and  8  must  be  used  in  the  aw  expression  corre sponding  to  Eq.  (30) .  In 
general,  Y  *  X  and  8^  f  8,  butWA  is  the  same  everywhere  and  requires 
no  special  subscript. 

Man  Plow 

The  mass  flow  through  the  rotor  disc  can  be  obtained  by 
integrating  the  trapezoidal  velocity  distribution  over  the  entire 
disc.  The  result  is 


m  = 


m  = 


1 2ns  p  v(| )  ds 
o 

ittR20R  ^2^*  xa(x)dx  j 


(31) 


(32) 


Mf  “ 


i 


2  /  xa(x)dx 


(33) 


In  the  distant  wake,  Eq.  (52)  is  applicable  with 
substituted  for  R2  in  the  area  factor  (note:  the  R  in  the  fiR 
factor  remains  nonsubscripted  because  this  term  was  used  to  trans¬ 
form  v  and  V  to  their  nondiraensional  forms  and  A  ) . 


m  =  pnR  2ftR  l  2#  xa  (x)dx 

w  L  Jn  w  J 


OM 


Dividing  ( 5 4)  by  (32)  and  rearranging, 


A  Rw 

e  =  T 


if 

xa(x)dx 

V  /  xaw(x)d, 


(55) 


MflMPttt  fiflBittmfclflBI 

The  momentum  flux  of  the  air  that  is  eventually  encountered 
by  the  rotor  while  it  moves  in  the  distant  atmosphere  is 


F 


o 


=  -  mV 


06) 


where  positive  F  is  downward. 

In  the  distant  wake,  the  momentum  flux  is  calculated  by 
integrating  the  trapezoidal  distribution: 


or 


/*Rw  2/  s 

W  2™pvwhr)ds 

o  '  w' 

Nondimensionalizing  the  integral, 

bl  'xa*,x) 

y  xaw2(x) 


Fw  = 


F  =  TTR2c(flR)2e2 


] 

] 


(57) 


(58) 

09) 
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Clearly,  the  change  In  momentum  flux  must  equal  the  rotor  thrust. 


T  =  +  mV  +  TTR2p(nR)2e2 


\?f  xotw2(x)dx  ] 


(40) 


Nondimensionalizing  with  the  help  of  Eq.  (33), 


CT  =  Mf  A  +  e 


2[A2 


(x)dxj 


(41) 


An  alternative  form  to  (4l)  is 


CT  =  x[Aa(x)  +  e2aw2(x)]dx  (42) 

o 

Energy  (Power)  Considerations 

The  energy  flux  (power)  in  the  air  before  it  reaches  the 
rotor  (i.e.,  when  this  air  is  in  the  distant  atmosphere)  is  given 
by 


P  =  1/2  mV2 

o 

The  power  added  to  the  air  by  the  rotor  is 


(43) 


-  /EvW4T  =/''W(et) 


dx 


(44) 


Eq.  (44)  can  be  verified  from  either  force  and 
momentum  considerations  or  blade  element  considerations. 


The  assumed  blade  element  has  span  ds  and  produces  lift  dL. 
The  shaft  torque  produced  by  this  element  is 


dQ  =  sdL  sinct^ 

v  (x) 

For  small  a^,  x  sina^  x  tana^  x  ^  ‘  • 


Then  the  torque  is 


dQ  =  s 


v(x)dL 

fts 


v(x)dL 

a 


Multiplying  through  by  fl  and  integrating, 


Qn  =  PT 


-/*<*>“  -  f}ix)  Sr 


dx 


Eq.  (44)  can  be  easily  nondimensionalized.  Putting 


L  =  T: 


/•I  Id  C 

PR  =  TTR  p(flR)^y  a(x)^_i_Jdx  (45) 

o 

The  power  in  the  distant  wake  is  assumed  to  be  completely  embodied 
in  kinetic  energy  fluxj  i.e.,  the  ambient  pressure  in  the  distant 
wake  has  returned  to  atmospheric  pressure. 

In  the  distant  wake,  currents  flow  in  radial,  tangential, 
and  axial  directions.  The  power  could  be  broken  into  three  corre¬ 
sponding  portions  denoted  by  the  subscripts  r,t,  and  a. 


P 

rw 


+  P 


tw 


+  P 


aw 


(46) 


The  radial  flows  occur  because  of  vortices  (and  vortex 
cores)  in  the  wake.  Tangential  flows  occur  due  to  "swirl",  which 
is  caused  by  momentum  imported  to  the  slipstream  by  rotor  shaft 
torque.  The  axial  power  is  associated  with  the  assumed  trapezoidal 
velocity  distribution. 

Assume  for  the  present  that  P  +  P.«  P  .  The  axial  power 
comes  from  integration  of  the  trapezoidal  distribution: 
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1 


/•  w 

Pw  =/  (1/2)  pans  vHr)  ds 


(**7) 


In  the  usual  manner,  the  right-hand  side  of  (47)  is  non- 
dimensionalized . 


or 


Pw  =  (l/2)PTTRw2(nR)5  [^2^  xaw3(x)dxj 

o 

Pw  =  1/2  TTR2p(flR)5e2  xaw3(x)dx] 


The  following  balance  of  power  must  exist: 


m 


p  =  pD  +  p 
w  R  o 


Substituting  (4.3) ,  (45)  and  (48)  into  (49), 


0*9) 


(1 


/2)TTR2p(nR)5e2  |~2  /*  xaw5(x)dxl  =  TtR2p(nR)^f  a(x)[-~jdx  +  (l/2)mV2, 
•*b  W  “ 1  o 

/I  1  j  V 

xa  *(x)dx  -J  a(x)(-^— T j  dx  +  A2/  x  a(x)dx 

*n  •A* 


This  equation  can  be  refined  further  after  the  blade 
element  analysis,  which  follows,  provides  an  expression  for 
dC^/ds.  Looking  ahead  to  Eq.  (5$), 


dC 


dx 


T  =  %  x[e  -  a(x)3 


Hence, 


1  1  i 

xa^ (x)dx  =  ~ — y  x[0t  -  a(x)]a(x)dx  +  A2^*  xa(x)dx 


or 


50 


(50) 


Up  to  this  point,  we  have  discussed  only  the  mass  flow, 
momentum,  and  energy  in  the  slipstream,  using  the  assumed  trapezoi¬ 
dal  distribution  of  induced  velocity.  Eqs.  (33),  (34),  (33),  (42) 
and  (50)  summarize  the  principal  results.  The  next  section  develops 
the  blade-element  expressions  for  the  rotor. 


IV. 3. 


biade-eummt  combideratioms 


If  xil-8,  the  axial  velocity  is  constant  and  equal  to 
v-V.  In  this  region  the  thrust  per  blade  is 

dO^  =  1/2  p(fis)2ac^0  -  jds  Oss  <R(l-B)  (50 


Assume  the  twist  distribution* 


e  =  —  8, 

s  t 


Then 


p  I  ® 

“b  -  1/2  Pfns)^  ao|-£ns 


[6tnR  (v-V)l 
L  ns  ns  J 


ds 


(52) 


=  l/2  p(fis)ac|fiR0t  +  V  -  vj  ds 


(55) 


Comparing  this  expression  to  the  Kutta-Joukowski  expression  with 
T  =  L  ,  one  sees  that 


dT^ 

ds 


=  pvr 


(5*0 


r  is  constant  for  Os;x  si-6  and  is  given  by  the  expression 


T  =  1 


/2  ac  £nR  0^.  +  V  -  v] 


0  sx  <  1-5 


(55) 


Since  T-*0  at  the  tip,  the  T  strength  given  by  (55)  must 
shed  as  training  vortices  over  the  region  1-8<x<1. 


$ 

This  distribution  is  usually  called  "ideal 
uniform  inflow  (v  constant  over  the  disc), 
along  the  blade  bound  vortex  system,  i.e., 
for  ideal  twist. 


twist",  since  it  causes 
and  constant  circulation 
in  Eq.  55  r  is  constant 
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From  this  system,  it  seems  reasonable  to  assume  that  the  r 
strength  fed  to  the  vortex  core  is  given  by  (55)  •  Nondimens ionalizing 
(27),  assuming  constant  chord  c, 


r 


ca 

2b 


TTR2fi 


0,  +  A 


(56) 


Continuing  with  the  blade-element  theory,  from  Eq.  (55) 
dT^  =  l/2  p(fts)  ac  -  v(x)j  ds  (57) 

where  v(x)  has  been  substituted  for  v-V  in  Eq.  (55).  Note  that 
Eq.  (57)  holds  over  the  entire  blade,  unlike  Eq.  (55). 


■(x)  J  ds 

-  a(x)Jdx 

=  1/2  p{nR)STTR2a  jig-  (-£-)*[et  -  »•(*)]  a*  (58) 

For  all  b  blades,  the  thrust  expression  becomes 


dT  =  l/2  TTR2p(fiR)2  aa  x|g^  -  a(x)J  dx 

For  the  case  when  the  blade  chord  is  constant,  this  thrust 
expression  can  be  integrated  and  nondimensionalized  to  yield  the 
equation 


Note  that  the  integral  in  (59)  also  appears  in  (55)  •  Thus, 
an  alternative  form  of  (59)  is 
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ur.4. 


VQRUX  OOCT  MODEL 


The  trapezoidal  induced  velocity  distribution  shown  in 
Figure  6  is  composed  of  a  uniform  distribution  plus  two  idealized 
vortex  cores.  These  cores  have  a  tangential  velocity  V.  which 
increases  linearly  with  distance  from  the  center  of  the  vortex  core 
until  the  core  radius  R6/2  is  reached,  beyond  which  point  the  vortex 
is  assumed  to  induce  no  velocity  either  inside  or  outside  the  slip¬ 
stream.  Of  course,  this  is  a  very  simplified  representation  of  the 
actual  flow,  and  to  obtain  accurate  results  with  this  mode],  some 
consideration  must  be  given  to  the  characteristics  of  real  viscous 
vortices.  This  forms  the  main  topic  of  this  section. 

A  viscous  vortex  model  is  given  in  Reference  27  which  has 
the  tangential  velocity  profile 

I  -r2/4vet\ 

This  velocity  distribution  for  some  time,  t,  is  shown  in 

Figure  21 . 


For  r<r  ,  the  model  tends  toward  the  simple  inviscid 
vortex  model 


r 


2rrr 


For  r<r  ,  the  velocity  profile  tends  toward  that  of  a 
rigid  body  spinniSg  at  a  rate  o>,  where 


CJO 


V 

max 


r 

c 


(62) 


It  was  shown  in  Chapter  III  that  the  boundary  condition  required 
to  fit  the  peripheral  velocities  of  a  vortex  core  to  the  slipstream 
velocities  of  a  descending  rotor  is 


v  =  2Vt 

where  is  the  tangential  velocity  of  the  core. 


(63) 
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Figure  21 .  Vortex  Model. 


I 
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Assume  that  V  ^  =  V  as  depicted  by  Figure  21 .  Then 

Eq.  (63)  can  be  used  to  determine  rc,  which  will  later  be 

related  to  the  coordinate  6  on  the  trapezoidal  velocity  distributions. 

To  determine  r  ,  one  must  first  differentiate  Eq.  (61)  with 
respect  to  r:  c 


d 


V(r.t)  „ 

a 


) 


V _ occurs  when  this  derivative  vanishes,  leaving  the  result 


For  simplicity,  define 


e 


Then  (6*0  is  written 

1  +  2  £  *  eC 

The  approximate  numerical  solution  to  this  equation  is 

C  *  5A  -  t, 

as  can  easily  be  determined  by  direct  substitution. 


(64) 


(65) 


(66) 


(67) 
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Combining  (65)  and  (67), 


4v  t  =  — r 

e  i 


(68) 


8 


and  substituting  (68)  into  (6l)  to  get  V  , 

max 


max  2ttt 


(- 


ro2  <. 


Using  (66)  to  substitute  for  e  , 


V 


max  2tt  r 


(’"  1+’20 


)-+A"\ 

e 


(69) 


max 


*  "*e  (’  ^  0 


(70) 


Using  the  value  5/4  for  C  as  given  by  (67),  Eq.  (70) 
becomes 


V 

max 


Tiro- 

c 


(71) 


Assuming  that  V  =  V  ,  Eqs.  (63)  and  (71)  can  be 
combined  to  get  *  max 


or 


v  = 


2I_ 

7ttt 

c 


(IP-) 
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Eq.  (72),  when  applied  to  the  cores  near  the  rotor, 

can  be  manipulated  to  get  (since  6  =  2r  /K  ) 

0 


ior 

7TTRv 


(73) 


The  value  of  r  applies  for  a  single  blade  and  is  given 
by  Eq.  (56) . 


the  T  strength  of  the  cores  is  constant.  Thus, 


*  _  ior 

w  "  7TTR  v 
w  w 

An  alternative  form  for  (7*0  would  be 


(7*0 


R  _  ior 

w  "  7nRvwe 


(75) 


At  this  point,  it  is  appropriate  to  combine  Eqs 
(56),  (73)  and  (75). 


ior  10  /  aa  \_rt>2  / _  »  ^ 

7ttRv  "  7rrR(fiR)X  \  2b  /  R  n  ^et  A  “  x) 


/  5oa  \  0t  +  A  '  X 

(76) 

l  7b  )  X 

1  5oa  \  0t  +  A  -  X 

(77) 

<■> 

£ 

Having  related  the  core  model  to  the  flow  inside  and  outside  the 
slipstream,  the  next  step  in  calculating  the  rotor  thrust  is  to 
allow  for  tip  losses  due  to  three-dimensional  flow,  or  "spillage". 
This  is  discussed  in  Section  5.1but  before  beginning  that  topic,  it 
is  convenient  to  summarize  the  major  equations  obtained  up  to  this 
point  and  to  restate  them  in  convenient  nondimensional  forms. 
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Equation  Summary 


a(x) 

Flow  model 

a  (x) 
w  ' 


X  -  A 
I  (1-x)  -  A 


A 

T  ('-*) 
w 


-  A 


OS.  x  S  1-6 
1-6  x  £  1 


Oixi  1-8 
w 

1-5  mi 
w 


(30) 


Continuity  e 


x  a(x)dx 


x  ay(x)dx 


(35) 


Momentum 


xa(x)dx  +  2e^x  ay2(x)dx 


(42) 


2/*1 

Energy  e  J  x 

0 

aw5(x)  dx  «  - 

"F”  y  x^et  a(x)dx 

0 

+  A* 

i 

7  xa(x)dx 

''o 

(50) 

Blade  element 

n  oa 

°t  *  ~r 

x  o(x)dx  ) 

(^9) 

Vortex  Core  model 

f,  -  1 ?ga 

\  VA"X 

(76) 

6  '  \  7b 

/  X 
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Vortex  Core  Model 


(77) 


To  simplify  the  equations,  introduce  the  notation 


R  n  2f 1  xa(x)dx  =  R  (A,X,6) 

1  Jo 

(78) 

R2  =  if'  xa2(x)dx  =  R2(A,X,6) 

(79) 

W,  4  ijf'  xaw(x)clx  -  W,(A,Xw,6w) 

W2  -  if  xaw2(x)dx  =  Ws(A,Xw,6w) 
o 

(80) 

(81) 

Wj  •  2/  xow}  (x)dx  -  W}  (A,  \w,  6^ 

O 

(82) 

Hie  equations  become 

(83) 

CT  .  AH,  +  e2  W2 

(84) 

e^j  =  -?r  (et8,  -  Rg)  +  A®R, 

(85) 

cT  •  ^  (et  -  H,) 

(86) 

B  .  p2S\  9t  +  A  -  X 

5  \  7b  )  X 

(87) 
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(88) 


The  problem  we  have  to  solve  can  now  be  stated  as 

Given-  A, 9^,0, a, b>  Find:  e,CL,X,X  ,5,6 

w  ^ 
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IV.5. 


UPLQflS  FACTOR 

The  analysis  presented  to  this  point  has  not  included  a 
blade  tip  loss  factor.  The  purpose  of  the  following  presentation 
is  to  derive  a  method  for  inciting  such  a  factor. 

Eq.  (53)  expresses  blade  element  thrust  for  a  single 

blade  as: 


pac 

2 


ns 


praet  +  v-vj 


ds 


This  can  be  written  in  a  somewhat  more  general  form  as 


«b  ■  -¥■  (fiR)2x(et  +  A  -  X)ds 


(89) 


The  blade  loading  is  defined  by 


£|2-  (flR)2x(et  +  A  -  X) 


(90) 


Eq.  (89)  can  be  integrated  to  get  the  thrust  produced 
by  one  blade*  For  b  constant- chord  blades,  the  thrust  is 


T  -  l/4  bpac(nR)2R(0t  +  A  -  X) 

-  a/4  TTC2p(nR)2(et  +  A  -  X))  (91) 


* 

The  integration  assumes  uniform  inflow,  i.e.,  X  is  not  a  function 
of  x.  This  is  equivalent  to  replacing  the  trapezoidal  inflow 
distribution  by  the  classical  uniform  distribution.  We  do  this/ 
at  this  point  in  the  derivation  to  obtain  the  classical  tip  loss 
factor.  Subsequently  we  return  to  the  trapezoidal  distribution 
to  get  an  improved  tip  loss  factor. 
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The  thrust  coefficient  is  therefore  given  by  Eq.  (66)  as 


CT  *  Hr  (®t  +  A  ■ 


C_  can  be  expressed  in  terms  of  the  blade  loading  at  the 
tip,  by  com cining  this  equation  and  Eq.  (90)  and  setting  x  =  1  in 
Eq.  (90).  The  result  is 

ct  ■  S  — )  (98) 


In  Reference  27,  a  tip  loss  factor  for  hovering  rotors  is 
presented  in  terms  of  as  follows; 


B  =  1 


(93) 


The  "tip  losses"  occur  because  air  flows  from  the  bottom 
region  to  the  top  region  of  the  blades  at  the  tips.  The  flow  occurs 
because  of  lift  (pressure  gradient  between  the  upper  and  lower  blade 
surfaces).  Such  flow  "destroys"  the  lift  on  the  blade  elements  in  the 
imnediate  region  of  the  tips,  because  it  interferes  with  the  two- 
dimensional  lifting  mechanism  of  the  airfoil  section. 

Since  tip  losses  are  caused  by  flows  at  the  tips  due  to 
lift  at  the  tips,  it  seems  that  the  blade  distributed  loading  at 
the  tips  really  controls  the  quantitative  value  of  B,  not  the  overall 
rotor  thrust  coefficient  as  suggested  by  Eq.  (93)  •  However, 

Eq.  (93)  is  valid  at  hover  (assuming  the  important  factor  is 
distributed  loading  at  the  tip),  since  Cm  is  directly  related  to  T' 
by  Eq.  (64) .  A  more  general  form  of  (93;  can  be  expressed  in 
terms  of  T/,  and,  as  shown  below,  this  leads  to  a  simple  expression 
for  B  valla  both  at  hover  and  in  descent. 


Eq.  (94)  has  been  derived  for  rotor  systems  with 
uniform  induced  velocity.  Because  B  is  expressed  in  terms  of 
T.I  ti_,  however,  it  seems  reasonable  that  it  should  be  applicable 
to  rotors  operating  with  nonuniform  inflow,  inasmuch  as  the  tip 
loss  is  a  local  phenomenon  depending  primarily  on  blade  tip  operating 
conditions . 
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For  the  descending  rotor,  Eq.  (57)  is  applicable. 


dT^  =  l/2  p(ns)ac|fiR9^  -  v(x)Jds 

(fiR)2x(9t  -  a(x))  (95) 

At  the  tip,  becomes* 

VI  tip  -  («H)2(et  +  A)  (96) 

Substituting  Eq.  (96)  into  Eq.  (94)  gives  the 
major  result  of  this  section,  a  general  formula  for  B. 

B  -  1  -  5  yf-(et  +  A)  (97) 


Note  that  X  disappears  from  the  ^  expression  at  the  tip.  This  is 
because  a(x)  =  -  A  at  x  =  1 ,  i.e.,  at  the  tip.  The  disappearance 
of  X  is  a  consequence  of  the  assumption  of  a  trapezoidal  inflow 
distribution.  It  accounts  for  the  difference  between  the  tip  loss 
factor  derived  here  (Eq.  97)  and  the  classical  tip  loss  factor, 

Eq.  93.  Because  the  trapezoidal  inflow  model  gives  a  higher  tC' 
than  the  uniform  inflow  model,  at  a  given  C_  and  A,  the  tip  loss 
is  larger,  i.e.,  B  is  smaller. 
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I V.6. 


INOORPORATICH  OF  THE  TIP  LOSS  FACTOR 


All  the  equations  required  to  solve  for  C_  versus  rate  of 
descent  of  a  given  rotor  have  now  been  obtained.  At  this  point,  the 
problem  is  to  include  the  effects  of  B  on  Eqs.  (83)  through  (88). 

One  way  to  do  this  is  outlined  below. 


Solve  Eqs.  (83)  through  (88)  in  a  normal  way,  without 
considering  tip  losses.  Results  of  such  a  solution  process  can  be 
expressed  in  nondimensional  form,  using  the  actual  rotor  radius  R  as 
the  characteristic  length. 

To  include  tip  losses,  say  that  the  effective  rotor  radius, 
R^  is  given  by 

Re  4  RB  .  R  [l  -  1  (8t+A)j  (98) 


The  solutions  to  (83)  through  (86)  can  be  considered  valid 
for  a  rotor  with  radius  R  and  no  tip  losses  (the  no  tip  loss 
assumption  being  already  embodied  in  Eq.  (83)  through  (88)). 

Thus,  when  the  nondimensional  data  representing  the  solutions  to 
Eqs.  (83)  through  (88)  are  dimensicnalized,  an  effective  radius  Rfi 

must  be  used  to  account  for  tip  losses. 

Obviously,  it  is  desirable  to  present  data  nondimens ionali zed 
by  a  constant  characteristic  length  (e.g.,  R)  rather  than  a  variable 
characteristic  length  (e.g.,  R  ).  To  accomplish  this,  the  non- 
dimensional  data  expressed  in  terms  of  R  must  be  restated  in  terms 
of  R.  As  an  example  of  this  process,  consider  a  graph  of  C  Tvs  A, 
predicted  for  a  rotor  of  given  geometry  by  solving  Eqs.  (83)  through 
(88) .  For  a  given  A,  the  actual  thrust  T  is  given  by 

T  *  TTRe2p(fiRe)2CT  =  prtl\kCT  (99) 
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Define  C_  as  the  desired  thrust  coefficient,  i.e.,  the  thrust  non- 
dimensionalized  with  respect  to  R.  Then 


T  =  pTTfiVc  * 


(100) 


Combining  (99)  and  (100)  to  eliminate  T 


:T*  ■  (-T-J 


Eqs.  (101)  and  (98)  can  be  combined  to  yield 


CT*  ■  ('  -  5  Pf  K  +  A))  °5 


(101) 


produce  A : 


The  characteristic  length,  R  ,  has  also  been  used  to 

© 


A 


(102) 


This  process  is  applicable  to  the  value  for  A  in  the  original  non- 
dimensionalized  data  C^,  =  C^(A) . 


The  question  arises,  does  the  nondimensionalizing  process 
(102)  apply  to  Eq.  (98)?  [if  R  is  included  in  the  A  in  Eq.  (98), 
then  Eq.  (98)  becomes  *n  implicit  expression  in  R  .]  It  is  expected 
that  the  difference  in  results  as  to  whether  the  A  in  (98)  is 
nandimensionalized  by  R  or  R  will  be  extremely  small.  Assuming 
this  for  the  time  being,  substitute  A*  into  (98)  and  (101), 

*  A  JL  (103) 

OR 
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However,  the  effect  of  B  should  be  included  in  determining 
the  A  in  the  original  data.  Combining  (102)  and  (103)  to  eliminate  V, 


A*  =  A  (tt)  =  (’  ■  ~rif¥~ lv"^))A  (,o4> 


which  can  be  rewritten  as 


(105) 


Summarizing,  the  major  results  of  this  section  are  the  following 
equations  : 


B  = 

1  *  $  y  t  (et + a*  > 

(97) 

A  = 

A* 

B 

(103) 

r  *- 
CT  “ 

b4ct 

(101) 

An  example  hand  calculation  of  C  *  versus  the  rate-of- 
descent  parameter  is  done  below,  to  show  the  method.  Note  that  to 
obtain  the  graph  of  Cj*  vs  J,it  is  necessary  to  assume  an  initial 
value  of  J*,  and  later  to  determine  what  value  of  J  this  implies. 

Thus  the  method  doesn't  predict  C  *  for  a  given  single  value  of  J,  but 
instead  predicts  the  entire  C^*  Versus  J*  graph. 
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The  derivations  given  in  this  chapter  employ  standard 
notation  used  in  helicopter  analysis#  A  different  notation  is 
traditional  in  propeller  analysis.  For  example,  J  is  used  to 
denote  nondimensional  axial  propeller  velocity.  J  is  related  to 
A  hy 


J  =  ttA 

The  propeller  thrust  coefficient  is  denoted  here  as  CT  and  is 

P 

related  to  the  ‘'helicopter"  thrust  coefficient  hy 


CT  =  64  tt3  0 
P 

In  the  comparison  of  our  predicted  versus  measured,  given  in 
the  next  section,  the  "propeller"  notation  is  used,  because 
the  experimental  data  are  presented  in  this  notation. 

The  results  of  the  calculation  are  shown  in  Figure  22 
and  are  compared  with  the  CT  versus  J  graph  that  would  be  obtained 
assuming  zero  tip  loss  factor.  The  difference  in  magnitude  of 
CT  and  slope  dCT/dJ  at  a  given  is  considerable. 
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TABUS  I.  EXAMPIE  CALCUIATION  OF  CT  VERSUS  J 

T 

0  =  .186  p  *  1.0  6  =  .25  a  =  5.73  b  =  3 

».^sm - T-r 

mum 

1 

9 

H 

PI 

S3 

f  x(3) 

(5) 

Vw 

(5)/b 

B 

i 

B2 

BU 

B 

-.40 

-.30 

— .  20 

-.10 

0 

•05 

.10 

.15 

.20 

•25 

-.127 

-.0955 

-.06366 

-.03183 

0 

.01592 

.03183 

.04775 

.06366 

.07958 

.0658 

.0828 

.0995 

.1160 

.1550 

.1420 

.1500 

•1590 

.1662 

.1750 

.256 

.289 

.315 

.341 

.365 

•577 

.387 

•398 

.410 

.420 

.0858 

.0965 

.1050 

.114 

.122 

.126 

•129 

.133 

•  137 

.140 

.9142 

•9035 

.8950 

.886 

.878 

.874 

.871 

.867 

.8^3 

.860 

-.437 

-.332 

-.223 

-.113 

0 

•0571 

.1145 

.17^0 

.2310 

.2910 

.123 

.152 

.182 

.210 

.236 

.249 

.260 

.271 

.280 

.287 

.838 

.818 

.800 

.785 

.769 

.765 

.760 

.751 

.746 

•7^0 

•  700 

.669 

.640 

.617 

•590 

.585 

•  579 
.565 

•  557 
.548 

.086 

.102 

.116 

.129 

.137 

.146 

.150 

.153 

.155 

.156 
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IV.7  CCMPARISQW  WITH  fXFKBJMSMT 

Hardly  any  data  on  thrust  of  helicopter  rotors  In  vertical 
descent  has  been  published.  Reference  19  describes  flight  tests  on  a 
Westland  Whirlwind  helicopter  In  vertical  descent*  but  unfortunately 
the  data  obtained  are  not  sufficiently  complete  to  permit  correlation 
with  our  theory.  We  therefore  have  to  employ  data  on  VTOL  propellers . 
Reference  18 presents  wind-tunnel  data  an  VTOL  propellers  at  shaft  angles 
of  attack  representing  vertical  and  inclined  descents,  and  Reference  28 
gives  data  on  the  same  propellers  in  climb  conditions  (though  not  all  the 
blade  settings  were  duplicated) .  By  combining  these  references  the 
experimental  data  graphed  on  Figure  23  were  obtained  (solid  lines) .  The 
dash-dot  lines  on  Figure  23  show  the  results  obtained  by  solving  Eqs . 

(78)- (89)  for  the  propeller  and  then  modifying  the  results  (as  specified 
by  Eqs.  (97)- ( 101)  to  include  the  tip  loss  effects. 

It  will  be  seen  that  the  theory  agrees  well  with  experiment 
for  climb,  hover,  and  rates  of  descent  outside  the  unsteady  region  of 
the  vortex-ring  state. 

Figures  2k  through  27  are  presented  to  show  the  comparison  be¬ 
tween  experiment  and  theory  for  the  following  different  theoretical 
approaches : 

(1)  Vortex  core  wake  model  with  variable  tip  loss  factor 
(Eq.  (97) )  - 

(2)  Variable  tip  loss  factor  with  no  vortex  cores  (i  .e . , 
wake  represented  by  uniform  downwash  over  rotor  disc 
of  radius  BR) . 

(3)  Conventional  analytic  model;  i.e.,  constant  tip  loss 
factor  (B  =  1  in  this  case)  and  uniform  downwash. 

For  all  three  of  these  cases,  the  solid  lines  represent  the  experimental 
data. 


Figures  2k  through  27  indicate  a  striking  improvement  in  the 
accuracy  of  theoretical  results  when  the  descending  rotor  models 
described  in  this  chapter  are  applied. 

Another  important  conclusion  that  can  be  reached,  by  inspec¬ 
tion  of  the  figures,  is  that  correlation  between  experiment  and  theory 
is  good  if  the  tip  losses  are  accounted  for  while  retaining  the  uniform 
inflow  model.  This  does  not  mean  that  the  vortex  cores  are  unimportant, 
since  this  core  model  was  an  essential  part  of  the  derivation  leading 
to  the  tip  loss  factor  (Eq.  (97) ) .  However,  once  the  existence  of 
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Figure  2k.  C  vs  J  for  6-Degree  Blade  Angle. 
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Figure  25.  CT  vs  J  for  9-Degree  Blade  Angle. 
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Figure  26.  CT  vs  J  for  12-Degree  Blade  Angle. 
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Figure  27.  C_  vs  J  for  15-Degree  Blade  Angle. 
P 
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the  viscous  "boundary  layer"  between  the  slipstream  tube  and  the 
surrounding  air  mass  is  recognized ,  the  simple  tip  loss  factor  can 
be  derived.  This  situation  is  quite  analogous  to  Frandtl's  fixed- 
wing  boundary  layer  model.  The  presence  of  the  viscous  region  must 
be  acknowleged  to  assess  correctly  the  physical  flow  mechanisms  which 
produce  the  wing  loads.  The  results  for  lift  and  induced  drag  derived 
from  inviscid  theory  (which  does  not  include  the  boundary  layer), 
however,  are  still  accurate  if  the  boundary  layer  has  not  separated, 
catusing  stall.  For  the  descending  rotor,  the  viscous  region  is 
analogous  to  the  boundary  layer.  The  nonuniform  downwash  indirectly 
causes  a  thrust  loss  by  inducing  higher  blade  loadings  near  the  tips. 

This  thrust  loss  greatly  exceeds  that  due  to  the  momentum  loss  associated 
with  the  nonuniform  downwash. 

IV. 8.  Summary 

The  performance  of  a  vertically  descending  rotor  can  be 
calculated  using  a  "vortex-core"  wake  system  and  a  suitable  blade  tip 
loss  factor.  The  results  obtained  by  analyzing  this  flow  model  show 
good  correlation  with  available  experimental  data  over  the  region 
from  fast  climb  to  descent  rates  outside  the  vortex-ring  buffet 
boundaries .  Although  the  vortex  core  model  is  required  to  derive  the 
descending  tip  loss  factor,  the  wake  inflow  distributions  can  be  assumed 
uniform  (retaining  the  tip  loss  factor  only  and  assuming  negligible 
vortex  core  radius)  without  seriously  degrading  the  accuracy  of  numerical 
results . 


The  application  of  the  tip  loss  factor  to  derivative 
calculations  is  discussed  in  Volume  II.  The  necessity  to 
consider  inclined  as  well  as  vertical  descent  requires  a  cycliely 
varying  B.  Experimental  data  are  not  sufficiently  complete  to  check 
the  accuracy  of  derivatives  calculated  using  this  approach;  however, 
it  follows  logically  from  the  assumption  that  conditions  at  the  tip 
determine  the  tip  losses. 


CHAPTER  V. 
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V.l.  IRTRCDUCnCil, 

This  chapter  describes  hew  stability  derivatives  were 
calculated  for  the  S-58  and  AH-56A.  (As  explained  below,  published 
derivatives  were  used,  for  the  YHC-1A.)  The  chapter  is  organized 
as  follows: 


Firstly,  the  MOSTAB  modular  stability  derivative  program 
is  described.  MOSTAB  was  specially  evolved  for  this  study;  it 
calculates  derivatives  for  any  type  of  helicopter,  with  a  wide 
variety  of  possible  flow  models.  The  description  in  this  chapter 
is  very  general;  full  details  of  the  equations  used  in  MOSTAB  and 
a  listing  are  given  in  Volume  II. 

Secondly,  the  forms  of  the  input  data  required  for  MOSTAB 
are  summarized,  and  the  specific  data  for  the  S-58  and  AH-56A  are 
presented.  These  data  include  three-view  drawings,  estimated  aero¬ 
dynamic  coefficients  for  fuselages,  rotors,  etc.,  and  the  appropriate 
reference  areas. 


Thirdly,  the  flight  conditions  for  which  derivatives  were 
calculated  are  tabulated.  (The  actual  derivatives  for  each  flight 
condition  are  listed  in  Volumes  III  and  IV,  with  the  transfer 
functions . ) 


Finally,  the  chapter  concludes  with  a  brief  discussion 
of  the  accuracy  of  the  MOSTAB  derivatives. 


V.2.  THE  NEED  FOR  MOSTAB 

During  the  course  of  the  study  described  here,  it  became 
apparent  that  it  was  difficult  to  obtain  reliable  stability  derivatives 
for  single-rotor,  tandem-rotor,  and  compound  helicopters  in  descending 
flight.  In  some  cases  derivatives  were  available  for  level  flight, 
but  the  programs  used  to  calculate  these  derivatives  contained  "  built 
in"  assumptions  regarding  flow  models  which  were  inappropriate  for 
descending  flight.  It  was  not  feasible  to  modify  these  programs  to 
incorporate  alternative  flow  models,  such  as  those  described  in 
Chapter  IV.  Furthermore,  only  limited  experimental  data  on  derivatives 
in  descending  flight  were  available. 
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Accordingly,  it  was  decided  to  start  afresh,  by  writing 
a  "universal”  program  to  compute  trim  and  stability  derivatives 
for  any  type  of  helicopter,  with  little  or  no  restriction  as  to 
the  flow  models  that  had  to  be  used.  Further  requirements  were 
that  the  stability  derivatives  should  be  printed-out  in  all  axis 
systems  of  interest,  in  dimensional  form  and  also  divided  by  the 
appropriate  inertias,  and  that  the  contribution  of  each  component 
to  the  complete  derivative  should  be  printed-out  so  that  the 
reasons  for  unexpected  values  of  derivatives  could  be  traced. 

The  result  of  these  requirements  is  the  MOSTAB  program  described 
below. 


The  following  description  of  MOSTAB  is  the  bare  minimum 
required  to  follow  the  subsequent  discussion.  For  more  details 
of  MOSTAB,  see  Volume  II. 

v.3.  PRINCIPAL  imam  or  MOBTAB 

MOSTAB  currently  calculates  performance,  trim,  and  stability 
derivatives  for  any  type  of  subsonic  aircraft.  The  program  can  be 
extended  for  use  in  hybrid  simulation  and  may  also  be  used  to 
calculate  limiting  conditions  such  as  'g'  boundaries.  The  unique 
feature  of  MOSTAB  which  gives  such  versatility  is  its  modular 
construction,  explained  below. 

Figure  28  shows  a  helicopter  separated  into  elements.  The 
physical  characteristics  of  each  element  are  specified  in  the  input 
data.  Separate  sections  of  MOSTAB  handle  each  element  type.  The 
influence  of  each  element  on  the  aircraft  is  summed  to  calculate 
dynamic  characteristics  of  the  assembled  flight  vehicle. 


When  an  element  (e.g.,  wing,  propeller)  produces  a  force  (lift, 
drag,  thrust),  the  air  in  the  vicinity  of  the  aircraft  is  set  in 
motion.  The  induced  velocities  (sometimes  called  downwash  and 
sidewash)  affect  other  elements  by  changing  their  local  airspeed 
and  angle  of  attack.  Interference  velocities  are  thus  very 
important  and  must  be  calculated  accurately. 

MOSTAB  includes  all  of  these  "aerodynamic  coupling" 
(interference  velocity)  effects  in  a  single  subroutine. 

Any  desired  model  for  interference  velocities  can  be  included  in 
this  subroutine, which  performs  all  of  the  necessary  interference 
velocity  calculations.  No  fundamental  program  changes  are  required 
to  modify  interference  velocity  models  as  required  to  suit  particular 
vehicles;  only  the  interference  subroutine  need  be  changed. 
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ASSEMBLED  AIRCRAFT 
IS  SEPARATED  INTO 
ITS  COMPONENTS 


HORIZONTAL 

STABILIZER 


FUSELAGE 


INTERFERENCE  VELOCITY 
CALCULATION 

•  Receives  Loads  Produced 
By  Each  Component 

•  Generates  Interference 
Velocities  Between  Components 


/ 


TAIL  ROTOR 


Figure  28.  Modular  Representation  Used  by  MOSTAB. 
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Either  rigid  rotors  (propellers)  or  flexible  bladed  rotors 
(helicopter  rotors  -  hinged  or  hingeless)  can  be  studied  by  MOSTAB. 

The  modular  nature  of  the  program  allows  one  to  represent  an  air¬ 
craft  with  as  many  rotor s,  aerodynamic  bodies  (fuselages,  nacelles, 
etc.), and  lifting  surfaces  (wings,  empennage  surfaces,  canard 
surfaces,  etc.)  as  desired. 

Because  aerodynamic  rotors  are  the  most  complicated  element 
types,  a  brief  description  of  the  rotor  analyses  incorporated  in 
MOSTAB  is  presented  here.  "  Rotors"  are  either  of  the  flexible 
bladed  type  (e.g.,  helicopter  rotors  with  hinged  blades,  or  the  so- 
called  "  rigid"  rotors,  whose  blades  deflect  structurally  to  a 
substantial  degree)  or  of  the  truly  rigid  type  (propellers)  •  The 
same  subroutines  compute  flexible  and  rigid  rotor  characteristics 
in  MOSTAB.  Since  the  rigid  rotor  (propeller)  blades  do  not  deflect, 
they  provide  no  particular  analytic  difficulty.  MOSTAB  numerically 
integrates  computed  inertial  and  aerodynamic  loads  radially  and 
azimuthally  to  determine  the  characteristics  of  rigid  rotors  (propellers) . 
The  blade  motion  of  f?.exible  rotors  must  be  determined  before  loads 
from  such  rotors  can  be  determined.  The  "normal  mode"  method  is 
used  in  MOSTAB  to  represent  the  flexible  blades.  The  normal  mode  shape 
and  frequency  for  the  first  flapping  mode  of  the  flexible  blade  are 
inputted  to  MOSTAB.  (Other  modes  can  be  added  easily  to  MOSTAB  if 
required,  but  the  first  mode  adequately  represents  flexible  bladed 
rotors  when  determining  stability  derivatives.)  Generalized  mass 
and  force  are  generated  in  the  program,  and  the  flapping  motion  of 
the  blades  is  then  computed  numerically.  Blade  motion  histories 
for  all  flexible  blade  rotors  are  printed  out  for  the  trim  flight 
condition. 

Other  Versions  of  MOSTAB 

As  the  MOSTAB  program  undergoes  more  development,  various 
versions  will  appear,  each  useful  for  analyzing  aircraft  under 
certain  conditions.  The  first  version,  MOSTAB-A,  is  considered 
obsolete  and  unusable.  This  report  describes  the  MOSTAB-B  version, 
which  will  always  have  applicability  if  its  basic  assumptions  and 
aerodynamic  models  are  appropriate  for  the  vehicle  and  flight 
regimes  being  studied.  More  advanced  versions  (-C,  -D,  etc.)  will 
subsequently  appear,  each  with  more  refined  aerodynamic  models,  or 
with  special  effects  added  for  specific  vehicle  studies  (e.g., 
tilt-wing,  ducted- fan,  etc.,  configurations). 
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As  of  June  1970,  MOSTAB-C  had  been  completed,  together 
with  a  user's  manual.  The  principal  improvements  of  MOSTAB-C 
over  MOSTAB-B  are: 

(1)  Airfoil  data  can  be  automatically  corrected  for 
compressibility  effects  by  means  of  a  subroutine 
which  reads  local  Mach  number.  The  subroutine 
applies  to  both  stalled  and  unstalled  flow  regimes. 
This  compressibility  correction  is  important  at 
moderate  and  high  speeds.  MOSTAB  was  originally 
developed  for  approach  conditions;  hence,  com¬ 
pressibility  is  not  explicitly  included  in  MOSTAB-B 
although  an  overall  correction  can  be  made  by 
modifying  the  two-dimensional  lift-curve  slope  in 
the  wing  and  rotor  input  data. 

(2)  MOSTAB-C  has  an  option  to  print  out  forces  and 
moments  acting  on  each  element  due  to  overall 
vehicle  velocity  and  the  induced  velocity  of  the 
element  itself,  but  excluding  induced  velocities 
due  to  other  elements.  This  is  important  for 
correlating  "  isolated  rotor"  tests  and  calculations 
with  data  determined  for  complete  aircraft. 

(3)  Rotor  integration  subroutines  have  been  refined  by 
removing  a  number  of  small- angle  assumptions. 

(4)  Derivatives  can  be  called  in  "stability  axes"  for 
any  flight  condition  including  hover  (this  is  done 
by  extending  the  definition  of  stability  axes  to 
include  hover). 

(5)  A  number  of  improvements  and  extensions  to  the 
data  input  and  print- out  formats  have  been  made. 

v.it  MOBttB  mamma  verity  BtraRewimrc 

The  interference  velocities  between  the  various  modules 
are  calculated  by  an  extension  of  the  classical  Glauert  expression 
which  applies  to  fixed  wings  as  well  as  rotors.  The  Glauert 
expression  is 


v 


(102) 


82 


where  T 


v 

P 

A 


Resultant  aerodynamic  force  on  the  module  (e.g.,  wing 
or  rotor) 


The  induced  velocity  at  a  specified  location  (e.g.,  down- 
wash  at  the  rotor) . 

Air  density 


(  =  jtR  for  a  rotor  of  radius  R, 
for  an  elliptic  wing  of  span  b) 


"  Characteristic  Area” 
jt  b2 


VT  =  Resultant  velocity  at  the  module 

k  =  A  factor  which  relates  the  downwash  at  the  desired 

location  to  the  downwash  at  the  module.  Thus  at  some 
distance  away  from  a  rotor,  when  the  slipstream  has 
fully  contracted,  k  =  2,  whereas  in  the  plane  of  the 
rotor  k  =  1. 


MOSTAB  extends  the  above  Glauert  expression  to  a  multidimensional 
form  so  that  all  linear  and  angular  components  of  the  "  wash"  can  be 
related  to  the  forces  and  moments  acting  on  the  module.  Thus  such  effects 
as  swirl  due  to  rotor  torque,  and  nonuniform  downwash  can  be  included. 

The  generalized  form  of  Eq.  (102)  consists  of  two  matrix  equations: 


(103) 


(104) 


where  v.  =  A  6-component  vector,  the  elements  of  vhlch  are 

the  linear  and  angular  induced  flow  velocities  at 
module  i 


VT  J  =  Absolute  magnitude  of  the  airspeed  at  module  i 

A.  1  =  A  6  x  6  matrix  which  is  input  to  MOSTAB.  Its 
.  1  J  elements  have  units  of  (i/area)  and  represent 
the  reciprocals  of  the  characteristic  areas  of 
module  i. 
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=  A  6-component  vector,  the  elements  of  which  are 
the  X,  Y,  Z  forces  and  L,  M,  N  moments  acting 
on  module  i. 

=  A  6 -component  vector,  the  elements  of  which  are 
the  linear  and  angular  induced  velocities  at 
the  location  of  module  j  induced  by  the  forces 
and  moments  acting  on  module  i  (e.g.,  i  could  be 
the  main  rotor  and  ,1  the  horizontal  tail) 

A  6  x  6  matrix  of  interference  velocity  coupling 
constants 

In  the  complete  MOSTAB  program,  the  order  of  vectors  is 
increased  from  6  to  6n  where  N  is  the  number  of  modules  (e.g.,  N  ■  5 
for  a  helicopter  with  a  fuselage,  a  main  rotor,  a  tail  rotor,  a 
horizontal  tail,  and  a  vertical  tail).  Thus  a  wide  variety  of 
interference  effects  can  be  modeled  with  ease. 

V.5*  HQSIAB  HTPOTDA^ 

The  numerical  values  of  the  derivatives  calculated  by  NOSTAB 
for  the  S-58  and  AH-56A  are  given  in  Volumes  III  end  IV.  To  supplement 
these  derivatives,  this  section  presents  data  which  permits  the 
basic  assumptions  underlying  the  derivatives  to  be  checked.  That 
is,  relevant  dimensions  and  aerodynamic  coefficients  for  each  "module" 
(wing,  fuselage,  etc.)  are  tabulated,  together  with  the  assumptions 
regarding  interference  velocities  between  the  modules.  Because  of 
the  low  speeds  of  interest  to  the  present  study,  compressibility 
corrections  are  not  included.  Hence,  the  aerodynamic  coefficients 
for  each  module  are  the  same  for  all  flight  conditions. 

To  understand  the  input  data,  it  is  first  necessary  to 
review  the  axis  systems  used  in  MOSTAB. 

Figure  29  presents  a  three-view  drawing  of  the  Hughes  0H-6A 
helicopter  with  various  coordinate  systems  attached.  Hiese  axes 
will  now  be  discussed  under  individual  headings. 

(1)  Overall  vehicle  reference  axes  -  Axes  x,  y,  z  of  Figure 
29  we  the  "  overall  vehicle  reference  axes "  or  "  overall 
vehicle  axes".  These  coordinates  are  defined  in  the  main 
text  and  are  the  most  basic  axes  used  in  MOSTAB.  The 
positions  of  all  vehicle  elements  are  specified  in  the 
basic  MOSTAB  input  data  deck  in  vehicle  coordinates,  and 
the  first  set  of  stability  and  control  derivatives  is 
output  in  vehicle  axes. 
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(2)  Individual  element  axes  -  basic  element  axes  -  All 
coordinate  systems  shovn  by  Figure  29,  with  the 
exception  of  the  overall  vehicle  axes,  are  "  individual 
element  axes"  or  "  individual  element  basic  axes" 

These  axes  all  lie  parallel  to  overall  vehicle  reference 
axes,  but  with  origins  located  at  the  "  reference  points" 
of  their  respective  aircraft  components,  as  shewn  in  the 
figure.  The  FyA,  FyI  and  Fc  matrices  printed  as  part  of 

MOSTAB  output,  and  the  columns  OH,  WE,  WO,  VIO,  VAO  and  FO 
printed  as  the  trim-iteration  solution  all  refer  to  these 
basic  element  axes,  as  discussed  in  detail  in  Volume  II. 

(3)  Tndiwiduml  element  local  yy  ■  These  coordinate  systems 
are  not  shewn  on  Figure  29,  but  they  are  defined  in 
detail  in  Volume  II.  The  "individual  element  local 
axes"  have  their  origins  at  the  respective  element 
reference  points  (as  do  the  individual  element  basic  axes 
discussed  in  (b)  above),  but  they  are  rotated  from  the 
basic  axes  to  lie  in  a  convenient  position  for  analysis 
of  the  particular  aircraft  component.  Far  example,  the 
main  rotor  local  axes  are  rotated  from  the  basic  axes 
(coordinates  x  ,  y  ,  z  in  Figure  29)  so  that  their  z 
axis  lies  along  the  rotor  shaft  centerline.  Also  from 
Figure  29,  one  sees  that  the  tail  rotor  basic  axes  must 
be  rotated  approximately  90  degrees  abort  the  x_  axis 
to  become  the  tail  rotor  local  axes  (with  the  z  lixis 
along  the  tail  rotor  shaft). 

(U)  Stability  mi  -  These  coordinates  are  translated  and 
rotated  from  "  overall  vehicle  axes"  so  that  their 
origin  lies  at  the  aircraft's  center  of  gravity  and 
so  that  the  v  and  w  airspeed  components  in  stability 
axes,  at  trim,  vanish.  This  is  the  classical  definition 
of  stability  axes.  In  NOSTAB,  stability  axes  are  defined 
with  .their  z  axis  in  the  xz  plane  of  the  overall  vehicle 
axes.  For  hover,  stability  axes  are  arbitrarily  defined 
such  that  the  x  and  y  axes  are  horizontal  with  the 
stability  x  axis  parallel  to  the  projection  of  the  over¬ 
all  vehicle  x  axis  on  the  horizontal  plane. 

V.6.  NUMERICAL  DATA  FOB  S-58,  AH- 56 A,  and  THC-1A 

Figures  JO,  31,  and  32  show  3-view  drawings  of  the  S-58, 

AH- 56a  and  YHC-1A.  The  input  data  for  the  S-58  and  AH-5&A  are 
summarized  in  Table  II.  All  symbols  are  defined  in  Volume  II. 
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Figure  30.  Three-View  Drawing  of  Sikoraky  S-58. 
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TABLE  II  -  MOSTAS  INPUT  DATA 

AIRCRAFT  LOADING  CONFIGURATION  AND  INERTIA  (POUND»FOOT«SLUG) 


MOSTAB 

DEFINITION 

SIKORSKY 

LOCKHEED 

SYMBOL 

S-5A 

AH-56A 

GROSS  1 

HEIGHT 

• 1 560*05 

•) 187*05 

XCG 

CENTER  OF  GRAVITY 

.3300-00 

•  0000 

YGG 

RELATIVE  TO  OVERALL 
VEHICLE  COORDINATES 

.0000 

•  0000 

ZCG 

•  8750 

•  0000 

IXX 

• 1 150*05 

• 5895*04 

I XY 

•  0000 

•  0000 

IXZ 

MOMENTS  OF  INERTIA  WITH 

•  0000 

•  0000 

IYX 

RESPECT  TO  THE  CENTER 

.0000 

•  0000 

IYY 

OF  GRAVITY  WITH  THE  EXE 

. 5000*05 

•2750*05 

IYZ 

AXIS  ALIGNED  WITH  O.V.C. 

•  0000 

•  0000 

IZX 

(O.v.c.  =  OVERALL  VEHICIE 

•  0000 

•  0000 

IZY 

COORDINATES) 

•  0000 

•  0000 

IZ7 

.5500*05 

.2309*05 

CHARACTERISTIC  AREA  MATRIX  ELEMENTS  A(1  •l)»i<2t2)«A(3«3) 
ALL  OTHER  ELEMENTS  ARE  ZERO# 


SIKORSKY 

S-58 

LOCKHFFD 

AH-54A 

FUSELAGE 

.2770-01 

•1350-01 

HORIZONTAL  TAIL 

.4140-01 

•1080-01 

VERTICAL  TAIL 

•3000-01 

.2500-01 

MAIN  ROTOR 

•4100-03 

•4800-03 

TAIL  ROTOR 

•1400-01 

•1270-01 

WING 

.1750-02 

PROPELLER 

•1270-01 

< 


i 
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TAflLE  II  -  MOSTAB  INPUT  DATA  (COISITTNurm 
fuselage  CHARACTERISTICS  <RADIAN*EOOT. SECOND) 


MOSTAB 

SYMBOL 

definition 

SIKORSKY 

LOCKHEED 

S-58 

AH-S6A 

PSI  8 

•  0000 

THETA  B 

PHI  B 

ANGULAR  DISPLACEMENTS 

•  0000 

relative  to  overall 

VEHICLE  COORDINATES 

•  0000 

•  0000 

W 

•  0000 

•  0000 

A 

LINEAR  DISPLACEMENT 

.0000 

•  0000 

Y 

7 

RELATIVE  TO  OVERALL 
VEHICLE  COORDINATES 

•  0000 

•  0000 

c. 

•  0000 

•  0000 

AB 

REFERENCE  AREA 

.5400*02 

•1500*03 

CO 

drag  coefe.  alpha*o. 

•4600*00 

•1000-00 

Cl 

DRAG  COEEE/ALPHA 

.0000 

•3000-03 

C2 

drag  coeee/alphax alpha 

•  0000 

.6000-03 

CYO 

SI DEFORCE  COEFE, BETA-0# 

•  0000 

•  0000 

CY1 

-SI DEFORCE  COEEE/BETA 

•1400*01 

•1000*01 

CZO 

-Z-EORCE  COEFE. ALPHA-0. 

•  0000 

•2800-01 

CZ1 

-Z-EORCE  COEEE/ALPHA 

•7000-00 

•4000-00 

LB 

REFERENCE  LENGTH 

•4400*02 

•5260*02 

CMO 

PITCH  COEFE. ALPHA-0. 

•  0000 

•2000-02 

CM1 

PITCH  COEEE/ALPHA 

•3400-00 

•2500-00 

CNO 

TAM  COEFE  ALPHA-0. 

•  0000 

•  0000 

CN1 

TAM  COEEE/BETA 

-.6600—00 

-•6000-00 
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TABLE  II  -  MOSTAB  INPUT  DATA  (CONTINUED) 


HORIZONTAL 

TAIL  CHARACTERISTICS  (RADI ANfFOOT* 

SECOND) 

MOSTAB 

DEFINITION 

SIKORSKY 

LOCKHEED 

SYMBOL 

S-58 

AH-56A 

PSI  L 

ANGULAR  DISPLACEMENTS 

•  0000 

•  0000 

THETA  L 

RELATIVE  TO  OVERALL 
VEHICLE  COORDINATES 

•  0000 

.5240-01 

PHI  L 

•  0000 

•  0000 

X 

LINEAR  DISPLACEMENT 

-•2830*02 

-•2770*02 

Y 

RELATIVE  TO  OVERALL 
VEHICLE  COORDINATES 

•  0000 

•  0000 

7 

-.2000*01 

•  0000 

AW 

SECTION  LIFT  CURVE  SLOPP 

.8000*01 

•6000*01 

SW 

REFERENCE  AREA 

.1238*02 

.3450*02 

CHORD 

REFERENCE  CHORD 

.2220*01 

•3180*01 

COO 

SECTION  CD  AT  ALPHA*0« 

•  0000 

•  0000 

CD1 

SECTION  CD/ALPHA 

•  0000 

•  0000 

CD? 

SECTION  CO/ALPHAXALPHA 

•  0000 

•  0000 

AWCLO 

ALPHA  AT  CO-MINIMUM 

•  0000 

•  0000 

CAP  GAMMA 

DIHEORAL  ANGLE 

•  0000 

•  0000 

BW 

REFERENCE  SPAN 

•5550*01 

•1083*02 

LAMOA  w 

TAPER  RATIO 

•  0000 

.5700-00 

CMO 

CM  AT  ALPHA»0. 

•  0000 

•  0000 

CMA 

CM/ALPHA 

•  0000 

•  0000 
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TABLE  II  -  MOSTAB  INPUT  DATA  (CONTINUED) 

VERTICAL  TAIL  CHARACTERISTICS  (RAOT AN* FOOT* SECOND) 


MOSTAB 

DEFINITION 

SIKORSKY 

LOCKHEED 

SYMBOL 

S-5B 

AH-56A 

PSI  L 

ANGULAR  DISPLACEMENTS 

•  0000 

•  0000 

THETA  L 

RELATIVE  TO  OVERALL 
VEHICLE  COORDINATES 

•  6600 

•  0000 

PHI  L 

.1570*01 

.1570*01 

X 

LINEAR  DISPLACEMENT 

-.3100*02 

-•2666*02 

Y 

RELATIVE  TO  OVERALL 
VEHICLE  COORDINATES 

•  0000 

•  0000 

1 

-.3300*01 

.2770*01 

AW 

SECTION  LIFT  CURVE  SLOPE 

•6000*01 

•6000*01 

SW 

REFERENCE  AREA 

.2400*02 

•2460*02 

CHORD 

REFERENCE  CHORD 

•3BB0*01 

•4250*01 

CDO 

SECTION  CD  AT  ALPHA«0. 

•  0000 

•  0000 

CD1 

SECTION  CD/ALPHA 

•  0000 

•  0000 

CD2 

SECTION  CO/ALPHA/ALPHA 

•  0000 

•  0000 

AWCLD 

ALPHA  AT  CD-MINIMUM 

•  0000 

•  0000 

CAP  GAMMA 

DIHEORAL  ANGLE 

•  0000 

•  0000 

BW 

REFERENCE  SPAN 

.6200*01 

•5800*01 

LAMDA  W 

TAPER  RATIO 

•  0000 

•  0000 

CMO 

CM  AT  ALPHAsO. 

•  0000 

•  0000 

CMA 

CM/ALPHA 

•  0000 

•  0000 
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TABLE  I!  -  MOSTAB  INPUT  DATA  (CONTINUED) 
VINO  CHARACTERISTICS  (RADIAN# FOOT* SECOND) 


MOSTAB 

DEFINITION 

SIKORSKY 

LOCKHEED 

SYMBOL 

S-5R 

AH-56A 

PSI  L 

ANGULAR  DISPLACEMENTS 

•  0000 

THETA  L 

RELATIVE  TO  OVERALL 
VEHICLE  COORDINATES 

•1660-00 

PHI  L 

•  0000 

X 

LINEAR  DISPLACEMENT 

-.6500-00 

Y 

RELATIVE  TO  OVERALL 
VEHICLE  COORDINATES 

•  0000 

Z 

•3400*01 

AW 

SECTION  LIFT  CURVE  SLOPE 

•6000*01 

sv 

REFERENCE  AREA 

.1950*03 

CHORD 

REFERENCE  CHORD 

•7?5ft*01 

CDO 

SECTION  CD  AT  ALPHA*0# 

.8000-0? 

COl 

SECTION  CD/ALPHA 

-.7000-02 

CD2 

SECTION  CD/ALPHA/ALPHA 

•1470-00 

AWCLO 

ALPHA  AT  CD-MINIMUM 

.1750-01 

CAP  GAMMA 

DIHEDRAL  ANGLE 

•  0000 

RW 

REFERENCE  SPAN 

•2700*02 

LAMDA  W 

TAPER  RATIO 

•5100-00 

CMO 

CM  AT  ALPHA»0. 

-.6100-01 

CMA 

CM/ALPHA 

•  0000 
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TABLE  II  -  MOSTAB  INPUT  DATA  (CONTINUED! 

MAIN  ROTOR  CHARACTERISTICS  (RADI AN. FOOT. SECONO.SLUG! 


MOSTAR 

definition 

SIKORSKY 

LOCKHEED 

SVMROL 

S-5R 

AH-56A 

omega 

ROTATIONAL  VELOCITY 

•2325*02 

•2570*02 

BO 

CONSTANT  TIP  LOSS  FACTOR 

.*700-00 

•9700-00 

RA 

VARIABLE  TIP  LOSS  FACTOR 

.1680-02 

•1760-02 

RB 

EXTRA  TIP  LOSS  FACTOR 

•  0000 

•  0000 

DELTA  0 

SECTION  CD  AT  ALPHA>0« 

•7500-02 

•7500-02 

OFLTA  1 

CO/ALPHA 

•  0000 

•  0000 

DELTA  2 

co/alpha/alpha 

.3200-00 

•3200-00 

DELTA  3 

PITCH/FLAP  COUPLING 

•  0000 

•  0000 

THETA  1 

TOTAL  BLADE  TVIST 

-.13*5-00 

-•8730-01 

R 

BLADE  RADIUS 

.2800*02 

.2570*02 

SMALL  B 

NUMBER  OF  BLARES 

.6000*01 

•4000*01 

SMALL  A 

SECTION  LIFT  CURVE  SLOPF 

.6000*01 

•6000*01 

PSI  R 

ANGULAR  DISPLACEMENTS 

.0000 

•  0000 

THETA  R 

RELATIVE  TO  OVERALL 
VEHICLE  COORDINATES. 

.3300-01 

•  0000 

PHI  R 

.0000 

•  0000 

X 

LINEAR  DISPLACEMENT 

•2700-00 

•  0000 

Y 

RELATIVE  TO  OVERALL 
VEHICLE  COORDINATES 

•  0000 

•  0000 

z 

-.8200*01 

-.4150*01 

p 

NAT. FLAP. FREDUENCY/OMEGA 

•1020*01 

•1150*01 

ROE 

ESTIMATED  RETA#PSI*0. 

•  0000 

-.5000-02 

RDOE 

ESTIMATED  BETADOT. PSI*0. 

•  0000 

•1000-02 

PR 

BETA  PERTURBATION 

•5000-02 

.5000-02 

PRO 

BETADOT  PERTURBATION 

.1000-00 

.1000-00 

CHORO 

CONSTANT  CHORD  BLADES 

•1360*01 

.2320*01 
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TABLE  If  •  MOST  All  INPUT  DATA  (CONTINUED) 
MAIN  MOTOR  DISTRIBUTED  BLADE  PRORERT ICS- 


BLADE 

SIKORSKY  S-S8 

LOCKHEED  AH-56A 

STATION 

RAOIAL 

DISTRIBUTED 

RADTAL 

DISTRIBUTED 

NUMBER 

OISTANCE 

MASS 

OISTANCE 

MASS 

1 

•  0000 

•1360-00 

•0000*01 

•1000*01 

2 

•1000*01 

•1360-00 

.7700*02 

•3800-00 

3 

•8400*01 

•1360-00 

•1160*02 

•3800-00 

4 

•1340*02 

•1360-00 

•1540*02 

•2700-00 

S 

•1690*02 

•1360-00 

•1800*02 

•2700-00 

6 

•1960*02 

•1360-00 

•2050*02 

•2700-00 

7 

•2180*02 

•1360-00 

•2170*02 

•2700-00 

8 

•2370*02 

•1360-00 

•2310*02 

•2700-00 

9 

•2530*02 

•1360-00 

•2430*02 

•2700-00 

10 

•2670*02 

•1360-00 

.2570*02 

•2700-00 

11 

•2800*02 

•1360-00 

INITIAL 

FIRST 

INITIAL 

FIRST 

SHAPE 

FLAPPINO 

SHAPE 

FLAPPING 

MODESHAPE 

MODESHAPE 

1 

•  0000 

•  0000 

•  0000 

•  0000 

2 

•  0000 

•  0000 

-•1000-00 

-.3800*01 

3 

•  0000 

-•7400*01 

-.1500-00 

-.7180*01 

4 

•  0000 

-•1240*02 

-.1500-00 

-.1100*02 

5 

•  0000 

-•1590*02 

-.1500-00 

-.1350*02 

6 

•  0000 

-.1860*02 

-.1500-00 

-.1610*02 

7 

•  0000 

-•2080*02 

-.1500-00 

-.1730*02 

8 

•  0000 

-.2270*02 

-•1500-00 

-.1850*02 

9 

•  0000 

-•2430*02 

-.1500-00 

-.1980*02 

10 

•  0000 

-.2570*02 

-.1500-00 

-.2110*02 

11 

•  0000 

-•2700*02 

POINT  MASSES 

i- 

STATION 

RADIAL 

MASS 

RADIAL 

MASS 

NUMBER 

OISTANCE 

DISTANCE 

1 

•2800*02 

.3477-00 

•2570*01 

•7700*01 

2 

•8950*01 

•1000*01 

3 

,2570*02 

.5000-01 

INITIAL 

FIRST 

INITIAL 

FIRST 

SHAPE 

FLAPPING 

SHAPE 

FLAPPING 

MOOESHAPE 

MODESHAPE 

1 

•  0000 

-•2700*02 

-.7500-01 

-.7400-00 

2 

-.1500-00 

-•4880*01 

3 

-•1500-00 

-.2110*02 
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TABLE  I!  -  NOSTAB  INPUT  DATA  (CONTINUED) 

TAIL  ROTOR  CHARACTERISTICS  (PAD I AN* FOOT* SECOND* SLUG) 


MOSTAR 

DEFINITION 

SIKORSKY 

LOCKHEED 

SYMBOL 

S-SR 

AH-S6A 

omega 

ROTATIONAL  VELOCITY 

•1300*03 

•1300*03 

BO 

CONSTANT  TIP  LOSS  FACTOR 

•0700*00 

•9500-00 

RA 

VARIABLE  TIP  LOSS  FACTOR 

•  0000 

•  0000 

BR 

EXTRA  TIP  LOSS  FACTOR 

•  0000 

•  0000 

DELTA. 0 

SECTION  CD  AT  ALPHAbO. 

•1000*01 

•1000-01 

DELTA  1 

CO/ALPHA 

•  0000 

•  0000 

DELTA  2 

CD/ALPHA/ALPHA 

•3200-00 

•3200-00 

DELTA  3 

PITCH/FLAP  COUPLING 

-•1000*01 

•  0000 

THETA  1 

TOTAL  BLADE  TVIST 

•  0000 

•  0000 

9 

BLADE  RADIUS 

•4670*01 

•5000*01 

SMALL  B 

NUMRER  OF  BLADES 

•4000*01 

•V>oo*oi 

SHALL  A 

SECTION  LIFT  CURVE  SLOPF 

•6000*01 

•6000*01 

PSI  R 

ANGULAR  DISPLACEMENTS 

•  0000 

•  0000 

THETA  R 

RELATIVE  TO  OVERALL 
VEHICLE  COORDINATES 

.0000 

•  0000 

PHI  R 

.1570*01 

•1570*01 

X 

LINEAR  DISPLACEMENT 

-.3300*02 

-•2990*02 

Y 

RELATIVE  TO  OVERALL 
VEHICLE  COORDINATES 

-.1750*01 

-•6000*01 

7 

-.5750*01 

•  0000 

P 

NAT .FLAP. FREOUENCY/OMEGA 

.1000*0) 

•1000*01 

ROE 

ESTIMATED  RETA*PSI«0. 

•  0000 

.0000 

ROOF 

ESTIMATED  RETAD0T*PSI*0. 

.0000 

•  0000 

PR 

BETA  PERTURBATION 

,5000-02 

•5000-0? 

PRO 

BETADOT  PERTURBATION 

.1000-00 

.1000-00 

CHORD 

CONSTANT  CHORD  BLADES 

.6120-00 

.1170*01 
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TABLE  II  -  NOSTAB  INPUT  OATA  (CONTINUED) 
TAIL  ROTOR  DISTRIBUTED  blade  properties- 


BLADE 

SIKORSKY  S-58 

LOCKHEED  AH-56A 

STATION 

RAOUL 

DISTRIBUTED 

RADIAL 

DISTRIBUTED 

NUMBER 

DISTANCE 

NASS 

DISTANCE 

NASS 

1 

•  0000 

•2560*01 

•0000 

•1000-00 

2 

•1400*01 

•2560*01 

•1000*01 

•1000*00 

3 

•2230*01 

•2560-01 

•1500*01 

•1000*00 

4 

•2820*01 

•2560*01 

•2000*01 

•1000*00 

5 

•3270*01 

•2560*01 

•2500*01 

•1000-00 

6 

•3640*01 

•2560*01 

•3000*01 

•1000-00 

7 

•3950*01 

•2560*01 

•3000*01 

•1000-00 

S 

•4220*01 

•2560-01 

•4000*01 

•1000-00 

9 

•4460*01 

•2560-01 

•4500*01 

•1000-00 

10 

•4670*01 

•2560-01 

•5000*01 

•1000-00 

INITIAL 

rfRST 

INITIAL 

FIRST 

SHAPE 

FLAPPIN6 

SHAPE 

FLAPPING 

MODfSHAPE 

NODESHAPE 

1 

•  0000 

•  0000 

•  0000 

•  0000 

2 

•  0000 

•1400*01 

•  0000 

•1000*01 

3 

•  0000 

•2230*01 

•  0000 

•1500*01 

4 

•  0000 

•2820*01 

•  0000 

•2000*01 

5 

•  0000 

•3270*01 

•0000 

•2500*01 

6 

•  0000 

•3640*01 

•0000 

•3000*01 

7 

•  0000 

•3950*01 

•0000 

•3500*01 

8 

•  0000 

•4220*01 

•0000 

•4000*01 

9 

•  0000 

•4460*01 

•  0000 

•4500*01 

10 

•  0000 

•4670*01 

•  0000 

•5000*01 
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TABLE  11  •  MOST  All  INPUT  DATA  (CONTINUED! 

PROPELLER  CHARACTERISTICS  <RAOIAN*FOOTt SECOND. SLUG) 


MOSTAB 

DEFINITION 

SIKORSKY 

LOCKHEED 

SYMBOL 

S-5B 

AH-56A 

OMEGA 

ROTATIONAL  VELOCITY 

•1800*03 

BO 

CONSTANT  TIP  LOSS  FACTOR 

•*500-00 

BA 

VARIABLE  TIP  LOSS  FACTOR 

•0000 

«R 

EXTRA  TIP  LOSS  FACTOR 

•0000 

DELTA 

0 

SECTION  CD  AT  ALPHA«0« 

•8000-0? 

DELTA 

1 

CO/ALPHA 

•0000 

DELTA 

2 

CD/ALPHA/ALPHA 

•3200-00 

DELTA 

3 

pitch/flap  coupling 

•  0000 

theta 

1 

TOTAL  BLADE  TMIST 

•2600-00 

Q 

BLADE  RAOIUS 

•5000*01 

SMALL 

■ 

NUMBER  OF  BLADES 

•3000*01 

SMALL 

A 

SECTION  LIFT  CURVE  SLOPE 

•6000*01 

PS!  R 

ANGULAR  DISPLACEMENTS 

•  0000 

THETA 

R 

RELATIVE  TO  OVERALL 
VEHICLE  COORDINATES 

•1570*01 

PHI  R 

•  0000 

X 

LINEAR  DISPLACEMENT 

-•3129*0? 

Y 

RELATIVE  TO  OVERALL 
VEHICLE  COORDINATES 

•  0000 

Z 

•  0000 

P 

NAT. FLAP. FREOUENCY/OMEGA 

•  0000 

BOE 

ESTIMATED  BETA*PST»0. 

•  0000 

BOOE 

ESTIMATED  BETADOT*PSI«0. 

•  0000 

PB 

BETA  PERTURBATION 

•  0000 

PBD 

BETADOT  PERTURBATION 

•  0000 

CHORD 

CONSTANT  CHORD  BLADES 

•1000*01 
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TABLE  II  -  M0STA8  INPUT  DATA  (CONCLUDED) 
PROPELLED  DISTRIBUTED  BLADE  PROPERTTES- 


OISTRIBUTED  BLADE  PROPERTIES- 


BLADE 

SIKORSKY  S-5B 

LOCKHFTD  AM-54A 

STATION 

RADIAL 

DISTRIBUTED 

RADIAL 

DISTRIBUTED 

NUMBER 

DISTANCE 

MASS 

DISTANCE 

MASS 

1 

•  0000 

•1000*00 

2 

•1000*01 

•1000*00 

3 

•1500*01 

•1000*00 

A 

•?000*01 

•1000*00 

5 

•P500*01 

•1000*00 

6 

•3000*01 

•1000-00 

7 

•3000*01 

•1000-00 

B 

•4000*01 

•1000*00 

9 

•4500*01 

•1000*00 

10 

•5000*01 

•1000*00 

INITIAL 

EIRST 

INITIAL 

EIRST 

SHAPE 

FLAPPING 

SHAPE 

FLAPPING 

MOOFSMAPE 

MOOESHAPE 

1 

•  0000 

•  0000 

? 

•  0000 

•1000*01 

3 

•  0000 

.1500*01 

4 

•  0000 

•?000*01 

5 

•  0000 

•?500*01 

6 

•  0000 

.3000*01 

7 

•  0000 

•3500*01 

8 

•  0000 

•4000*01 

9 

•  0000 

•4500*01 

10 

•  0000 

.5000*01 
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TABIC  !!  -  HOST AS  INPUT  DATA  (CONTINUED) 
INDUCED  VELOCITY  INTERFERENCE  FACTORS- 
VELOCITIES  ARE  IN  OVERALL  VEHICLE  COORDTNATFS 
ALL  ELEMENTS  NOT  SHOWN  ARE  ASSUMED  ZERO 


TNOUCEO  VELOCITY 

COMPONENT 

INOUCFD  VELOCITY 

ORIGINATING 

FROM 

VELOCITY 

MULTIPLYING 

FACTORS 

LOADS  ON  THE 

AFFECTED  PY 

FOLLOWING 

THE  INDUCED 

SIKORSKY 

LOCKHFED 

COMPONENTS 

VELOCITY 

S-5G 

AH-S6A 

FUSELAGE 

(U) 

FUSELAGE 

(0) 

•  0000 

•1000*01 

(V) 

(V) 

•  0000 

•1000*01 

(W) 

(W) 

•  0000 

•1000*01 

(U) 

HORIZONTAL 

(U) 

•1000*01 

•oooo 

(U) 

VERTICAL 

(U) 

•1000*01 

•  0000 

(If) 

PROPELLER 

(U) 

•1000*01 

(U) 

TAIL  ROTOR 

(U) 

•1000*01 

(W) 

(W) 

•1S00*01 

HORIZONTAL 

(U) 

HORIZONTAL 

(II) 

•  0000 

•1000*01 

tail 

(V) 

TAIL 

(V) 

•  0000 

•1000*01 

(W) 

(W) 

•1000*01 

•1000*01 

WING 

(U) 

WING 

(U) 

•1000*01 

(V) 

(V) 

•1000*01 

(W) 

(W) 

•1000*01 

(W) 

MAIN  ROTOR 

(W) 

•1000*01 

VERTICAL 

(U) 

VERTICAL 

(If) 

•  0000 

•1000*01 

tail 

(V) 

TAR 

(V) 

•1000*01 

•1000*01 

(W) 

(W) 

•  0000 

•1000*01 

MAIN  ROTOR 

(U) 

MAIN  ROTOR 

(U> 

•1000*01 

.1000*01 

(V) 

(V) 

•1000*01 

•1000*01 

(W) 

(W) 

•1000*01 

•1000*01 

(W) 

WING 

(W) 

•1600*01 

(W) 

HORIZONTAL 

(W) 

•1600*01 

•  0000 

TAIL  ROTOR 

(U) 

TAIL  ROTOR 

(If) 

•  0000 

•1000*01 

(V) 

(V) 

•1000*01 

•1000*01 

(W) 

(W) 

•  0000 

•1000*01 

PROPELLER 

(U) 

PROPELLFR 

(U) 

•1000*01 

(V) 

cv> 

•1000*01 

(W) 

(W) 

•1000*01 
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For  the  YHC-1A  it  was  decided  not  to  employ  NOSTAB  to 
calculate  the  derivatives.  This  was  because  published  derivatives 
were  available  from  Ref.  29.  These  derivatives  were  calculated  by 
the  manufacturer  and  use  simple  down  wash  models,  not  including 
the  cyclic  variation  of  tip  loss  factor,  B.  Some  accuracy  is  lost 
because  of  this.  On  the  other  hand,  the  manufacturer  undoubtedly 
possesses  good  information  on  the  important  interference  effects 
between  front  and  rear  rotors.  This  would  have  to  be  estimated 
if  the  YHC-1A  derivatives  were  calculated  by  MOSTAB.  It  was  felt 
that  the  errors  incurred  in  such  estimates  would  exceed  those 
caused  by  oversimplified  tip  loss  models.  Thus,  the  manufacturer's 
derivatives  are  used. 

The  resulting  derivatives  are  presented  in  Volumes  I II  and 
IV  (together  with  the  transfer  functions)  for  the  flight  conditions 
tabulated  in  Table  III  below.  For  the  S-58  and  AH-56A,  the  maximum 
rate  of  descent  corresponds  approximately  to  the  upper  boundary 
of  the  vortex-ring  state. 


TABLE  III.  FLIGHT  CONDITIONS 


HELICOPTER 

RATE  OF  DESCENT 

(knots) 

(fi») 

AH-56A 

20,  40,  60,  100 

0,  9.2,  19.2, 

28.8,  at  each  speed 

Vertical  Descent 

0.  9.2,  19.2, 

28.8 

S-58 

20,  40,  60,  100 

0,  7.5,  15 

22.5  at  each  speed 

Vertical  Descent 

0.9,  9.2,  19.2, 

28.8 

60.1 

0,  25 

YHC-1A 

80.1 

0,  25 

Vertical  Descent 

0,  25 
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v.7  ooRREAnar  or  bboyatives  with  gaggE  MSI  MU 


Before  presenting  transfer  functions  based  on  the  published 
YHC-1A  derivatives  and  the  AH-5^A  and  S-58  derivatives  calculated 
by  MOSTAB,  it  is  highly  desirable  to  check  the  accuracy  of  these 
derivatives.  The  only  really  satisfactory  check  is  to  compare  cal¬ 
culated  derivatives  with  those  deduced  from  flight  tests.  Unfortunately, 
many  derivatives  are  hard  to  obtain  from  flight  tests,  and  only 
incomplete  flight  test  data  were  available  for  the  helicopters  of 
interest.  However,  an  incomplete  check  is  better  chan  none  at  all, 
so  this  section  sunmarizes  the  comparisons  that  were  made  between 
calculated  derivatives  and  the  available  flight  test  data. 

In  brief: 

(i)  For  the  YHC-1A,  the  results  of  steady  sideslip  flight 
tests  were  available.  These  substantiate  the  published 
derivative.  Y^,  L^,  Yy,  Iy,  H„,  and  Y6R,  L£>R,  N5R. 

(ii)  For  the  S-58  flight  test  data  shoving  time  histories  of 

responses  to  B^  c,  longitudinal  control  inputs  were  compared 

to  time  histories  predicted  from  the  calculated  derivatives. 
The  agreement  is  good  for  the  longitudinal  degrees  of 
freedom,  but  less  so  on  the  lateral  motions  induced  by 
longitudinal-lateral  cross-coupling.  Unfortunately,  it 
is  not  known  whether  any  lateral  control  Inputs  were 
Lpplied,  so  the  accuracy  of  the  cross-coupling  derivatives 
remains  uncertain.  However, the  longitudinal  derivatives 
appear  to  be  substantially  correct. 

(iii)  No  flight  test  data  were  available  for  the  AH-56A; 

wherever  possible,  manufacturer's  wind-tunnel  data  were 
used  as  inputs  to  MOSTAB,  and  derivatives  appear  to  be 
reasonable . 

In  general,  there  is  little  value  in  making  comparisons 
between  the  MOSTAB  derivatives  and  other  calculated  derivatives. 

This  is  because  the  other  derivatives  are  generally  calculated  by 
less  sophisticated  programs  which  do  not  include  all  the  interference 
velocities  between  components.  For  example,  Ref.  2b  presents  S-58 
derivatives  and  flight  test  data;  however,  the  derivatives  yield 
responses  which  do  not  agree  with  the  test  data.  Ibis  is  principally 
due  to  the  neglect  of  main  rotor  downwash  on  the  tail  and  tail  rotor, 
giving  inaccurate  Adequate  flight  test  data  are  difficult  to 
obtain,  but  they  constitute  the  only  valid  criteria  against  which 
calculated  derivatives  may  be  checked. 
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IHC-1A  flight  Tggtg 


Ref.  30  describes  steady  sideslip  tests  performed  on  a 
YHC-1A  in  level  flight  at  forward  speeds  of  45,  60,  and  80  knots. 
Derivatives  for  the  latter  two  speeds  are  presented  in  Volume 
III}  and  as  shown  in  Ref.  30,  the  steady- state  variation  of  lateral 
stick  position  (6.)  and  rudder  pedal  position  (&_)  with  sideslip 
angle  ( B)  can  be  predicted  from  the  following  equation,  using 
stability  axis  derivatives. 


(105) 


This  solution  applies  even  to  a  basically  unstable  aircraft  such 
as  the  YHC-1A,  provided  it  is  stabilised  by  feedbacks  from  yaw 
and  roll  rates,  these  feedbacks  being  supplied  either  by  the  pilot 
or  the  stability  augmenter  system.  Figures  33  and  34  (taken  from 
Ref.  31)  indicate  that  the  calculated  v,  &_,  and  6A  derivatives  are 
in  reasonable  agreement  with  flight  test.  Furthermore,  the  helicopter's 
behavior  is  fairly  linear  with  sideslip  angle  for  moderate  per¬ 
turbations. 


«-«  iu«ht  n«t« 

Seckel  (Ref.  25)  presents  geometric  and  inertial  data  on 
the  S-58,  wind-tunnel  data  on  components,  and  manufacturer's 
estimates  of  derivatives  for  the  helicopter  flying  at  5,000  feet 
with  a  gross  weight  of  11,600  pounds.  Flight  test  data  are  also 
shown,  i.e.,  time  histories  of  a  response  to  a  longitudinal  stick 
input  in  pitch  rate,  anfle  of  attack,  normal  acceleration,  and 
airspeed.  These  time  histories  sure  reproduced  from  the  original 
flight  test  report  (Ref.  32),  but  for  some  variables  the  reproduction 
is  not  sufficiently  precise  for  purposes  of  checking  derivatives. 

Ref.  32  should,  therefore,  be  used  in  preference  to  Ref.  25.  In 
addition,  Ref.  32  includes  the  time  histories  of  the  lateral  degrees 
of  freedom  (although  the  lateral  control  time  histories  are  not 
recorded ) . 
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PEDAL  POSITION,  PERCENT  OP  TOTAL  TRAVEL 


LATERAL  STICK  POSITION,  PERCENT  OF  TOTAL  TRAVEL 


FULL  RIOHT 


FULL  RIOHT 


FULL  RIGHT 


SIDESLIP  ANGLE,  DEGREES 

LEFT  RIOHT 


Figure  YHC-1A  Stick  Position  Versus  Sideslip  Angle. 


106 


The  time  history  response  following  a  step  B1  input 

from  trim  (-0.015  rad.)  at  74  knots  C.A.S.,  5000  ft.,  level 
flight,  was  selected  from  the  available  flight  test  histories 
of  Kef.  52.  The  appropriate  six- degree* of- freedom  stability 
derivatives  generated  from  MOSTAB  are  shown  in  Table  IV.  The 
derivatives  were  input  into  program  ZEPLIN,  which  solves  the 
six-degree-of-freedom, ,  small-  perturbation,  linear  equations  of 
motion  and  obtains  eigenvalues  (R  +  JI)  and  the  modulus  and 
angle  of  corresponding  residue  vectors  (L,A)  for  response  to 
a  unit  B^  step.  The  eigenvalues  and  residue  matrix  are  shown 

on  Table  V.  The  aircraft  response  was  then  obtained  using  the 
following  relationship  : 


e Rrt  (L  CosAj  Cos  I  t  -  (L  Sin  A)  Sin  I  t 
r  r  v  t  r 


n 


Response  *  B. 


i  E 


r  =  1 


006) 


The  comparison  of  the  MOSTAB-generated  and  experimental 
time  histories  for  the  selected  flight  condition  is  shown  on 
Figure  35*  For  longitudinal  motions,  there  is  excellent  agreement 
with  airspeed,  pitch  angle,  and  angle  of  attack,  but  poor  agreement 
with  pitch  rate.  This  discrepancy  is  attributed  to  malfunctioning 
pitch  rate  instrumentation,  because  the  measured  pitch  rate  does 
not  agree  with  the  differentiated  pitch  attitude. 

For  lateral  -  directional  motion  (Figure  36), the  flight  test 
results  show  an  initial  roll  and  yaw  to  the  right)  however,  MOSTAB 
derivatives  show  a  roll  and  yaw  to  the  left.  There  is  a  possibility 
that  the  pilot  applied  some  roll  control  as  he  pulled  back  the  stick, 
but  the  disagreement  may  be  due  to  too  large  a  MOSTAB-generated  value 
for  L  (-10,940  ft -lb  -sec  ),  causing  a  strong  left  rolling  moment 
as  thi  vehicle  pitches  up.  Further  studies  and  flight  test  comparisons 
will  be  required  to  produce  a  suitable  explanation  for  this  disagree¬ 
ment  between  time  histories. 


Summary 

The  MOSTAB  program  has  been  described  in  general,  avoiding 
technical  detail  which  is  given  in  Volume  II.  Flight  test  data  are 
limited,  but  they  generally  support  the  values  of  derivatives  used 
for  the  YHC-1A  and  S-58. 


107 


TABLE  IV  S-58  DERIVATIVES  FOR  FLIGHT  TEST  CONDITION 


SIKORSKY  S-58  73  KNOTS  5000  FT.  LEVEL  FLIGHT  FWD.  C.O. 

SEE  REFS.  25  AND  32 

VEHICLE 

HEIGHT  SCO  YCG  ZGG 
•11*7*05  .4400-00  .0000  .0000 

STABILITY  DERIVATIVE  MATRICES 

U  V  H  R  8  R 

X  -.1361*02  -.2267-00  -.1557*01  -.7005*03  .1244*04  -.1142*03 

Y  .1098*01  -.3349*0?  -.2662*01  -.1360*04  -.6969*03  .4766*03 

Z  -.7200*01  -.3777*01  -.2163*03  -.7174*03  -.3467*02  .9191*03 

L  .2137*02  -.2062*03  .1326*02  -.2033*05  -.1094*05  .2179*04 

M  .1106*03  .1790*02  .5892*02  .1070*05  -.2400*05  .1189*04 

N  -.7655*02  .3494*03  -.2985*03  -.4756*03  .1070*05  -.1969*05 

U  DOT  V  DOT  W  DOT  R  DOT  8  DOT  R  DOT 

X  .3690-02  -.8526-03  -.4382-01  -.2570*02  -.9115*01  -.2130*01 

Y  -.7242-02  .6685-01  .8636-01  .9947*01  -.2444*02  -.1500*01 

Z  .1053-03  -.1471-02  -.6709-03  -.6369-00  -.6005-00  -.2350-01 

L  -.1981-00  .5054-00  .2353*01  .2673*04  -.3945*03  .2080*03 

H  -.1607-00  -.5245-01  .1959*01  .4154*03  .2672*04  .3676*0? 

N  .8099-02  -.2248*01  -.1042*00  -.4841*02  -.5217*02  .7643*0? 

C(  1)  C<  2)  C(  3)  C<  4) 

X  -.3872*04  -.6976*03  .1218*05  -.1047*03 

Y  -.6522*03  .1258*05  .1009*04  .4204*04 

Z  -.9315*05  -.3025*02  .2802*05  -.1032*03 

L  .6394*03  .1882*06  .3210*04  .3163*05 

M  .9605*05  .5190*04  -.1940*06  -.3360*04 

N  .7851*05  .4805*03  .3463*05  -.1414*06 

THE  INERTIA  TENSOR  .5940*04  .0000  -.8800*03 

•0000  .2750*05  .0000 

-.8800*03  .0000  .2304*05 

TRIMMED  VELOCITIES  WITH  RESPECT  TO  OVERALL  VEHICLE  REFERENCE  AXES- 
U  V  W  R  0  R 

.1230*03  -.0000  -.6314*01  -.0000  .0000  -.0000 

TRIMMED  ITERATION  COLUMN  VECTOR.  TE- 
.2507-00  -.1755-01  .2118-01  .8137-01  -.5127-01  -.1546-01 

STABILITY  AXIS  SYSTEM  EULER  ANGLES-  THETA*  .1220-04  PHI  -.1544-01 
AIRCRAFT  INERTIAL  SPEED*  .1232*03 
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TMLF  V  S-56  RESIOUES  AND  EIGENVALUES  FOP  FLIGHT  TEST  CONDITION 


p 

I 

-.1178*00 
•  0000 

-•6261*01 
•  0000 

-•6549-00 
•  0000 

•4182-01 

•3269-00 

•4162-01 

-•3269-90 

-•3367-00 

•1316*01 

L  •1178*00 
A  .1800*03 
VECTORS 

•6281*01 

•1800*03 

•6549-00 

•1600*03 

•3316-00 

•8275*02 

•3316-00 

-•6275*02 

•1356*01 

•1043*03 

U  L 

A 

•4269*01 
•  0000 

•6361-01 
•  0000 

•7156*02 
•  0000 

•8357*03 

•1754*03 

•8357*03 

-•1754*03 

•2774*01 

-•9867*02 

V  L 

A 

•1443*02 
•  0000 

•6195-00 
•  0000 

•1719*04 
•  0000 

•3606*03 

•1592*03 

•3606*03 

-.1592*03 

•1377*03 

-.1795*03 

N  L 

A 

•8962-00 
•  0000 

•4391-00 
•  0000 

•1775*04 
•  0000 

•4901*03 

•1308*03 

.4901*01 

-.1308*03 

.2077*02 

-.1427*03 

P  L 

A 

•1140*00 

•1800*03 

•2586-00 
•  0000 

•1626*02 

•1800*03 

•4366*01 

-.2918*02 

•4366*01 

•2918*02 

.1251*01 

-.3267*01 

0  l 

A 

•5689-02 

•1800*03 

•1626-01 

•1800*03 

•1407*01 

•1800*03 

•2941*01 

.1572*03 

•2941*01 

-.1572*03 

•2551-00 

-.6254*02 

P  L  .2623-00 

A  .0000 
EIGENVALUES 

•1257-01 
•  0000 

•1492*02 
•  0000 

.2467*01 

-.1046*03 

,2467*01 

•1046*03 

•1322*01 

.1077*03 

P 

I 

-•3357-00 

-.1316*01 

-.1267*01 
•  0000 

•0000 
•  0000 

•  0000 
•  0000 

•  0000 
•  0000 

.0000 
•  0000 

L  .1358*01 

A  -.1043*03 
VECTORS 

•1267*01 

•1800*03 

•  0000 
•  0000 

•  0000 
•  0000 

.0000 

.0000 

•  0000 
•  0000 

U  L 

A 

•2774*01 

.9867*02 

.1487*03 

•1800*03 

•1502*03 

•1800*03 

•  0000 
•  0000 

.0000 
•  0000 

.0000 
•  0000 

V  L 

A 

.1377*03 

.1795*03 

.9137*03 

•1800*03 

•  0000 
•  0000 

•1084*05 
•  0000 

•  0000 
.0000 

.0000 
•  0000 

W  L 

A 

.2077*02 

.1427*03 

•1158*04 

.1800*03 

•  0000 
.0000 

•  0000 
•  0000 

.1160*03 
•  0000 

•  0000 
•  0000 

P  L 

A 

.1251*01 

.3267*01 

.5993*01 
•  0000 

•0000 
•  0000 

•0000 

.0000 

•  0000 
.0000 

.3745-03 
•  0000 

0  L 

A 

.2553-00 

.6254*02 

.6573*01 
•  0000 

•  0000 
•  0000 

.0000 
•  0000 

.0000 
•  0000 

.4741-00 

•  0000 

P  L 

A  - 

.1322*01 

••1077*03 

.1028*02 

•1800*03 

.0000 

.0000 

•  0000 
•  0000 

•  0000 
.0000 

.3070*02 

.1800*03 
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Longitudinal  Stick  Angle  of  Attack  •  Pitch  Angle :  Pitch  Sate  C.A.S 

Position,  In;  A  O,  Deg  A#,  Deg  q.  Deg  /Sec  — j  Knots 


Sikorsky  8-56  Helicopter 
Trim  Speed  Jk  Knots  C.A.S. ,  $000  Feet 
'  W  -  11, 9^*  Lb#  Forward  C.O.  130.5  In. 


Figure  35.  A  Comparison  of  Longitudinal  Flight  Test  Response 
Following  a  Step  Pitch  Control  Input  to  a 
Response  Based  on  Stability  Derivatives  Obtained 
From  MOSTAB. 
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Figure  36.  A  Comparison  of  Lateral  Flight  Test  Response 

Following  a  Step  Pitch  Control  Input  to  a  Response 
Based  on  Stability  Derivatives  Obtained  from  MOSTAB. 


Ill 


CHAPTER  VI.  TANDEM- ROTOR  HELICOPTER 

WEE'-LMF  wmm - 


In  this  sactlon  we  consider  the  transfer  functions  of  the  YHC-1A, 
which  Is  taken  to  represent  a  typical  tandem-rotor  helicopter.  The 
transfer  functions  are  presented  in  Volume  III.  Unlike  the  single¬ 
rotor  and  compound  configurations^  for  the  tandem  configuration  lateral- 
longitudinal  cross-coupling  effects  were  neglected)  thus,  separate 
lateral  and  longitudinal  transfer  functions  were  calculated.  The 
format  of  the  transfer  function  print-out  is  explained  below. 

The  format  is  basically  similar  to  that  of  Reference  33,  but  a 
number  of  detailed  refinements  have  been  added.  These  include  addi¬ 
tional  sets  of  numerators  referred  to  a  station  distant  from  the  c.g. 

The  coordinates  of  this  station  are  conveniently  chosen  as  the  location 
of  the  pilot,  to  relate  pilot  cues  to  the  motion  of  the  c.g.  of  the 
aircraft. 

VI.  2  LCBQITUDIKAL  DYNAMICS 

As  noted  in  Volume  III,  the  derivatives  are  inputted  in  body 
(i.e.,  waterline)  axes.  The  input  format  of  the  program  used  contains 
a  number  of  redundancies,  e.g.,  Mach  number,  which  are  arbitrarily  put 
to  zero. 

The  print-out,  in  Volume  III,  gives: 

(i)  Derivatives  in  stability  axes 

(ii)  Eigenvalues,  and  associated  information,  such  as 
the  coefficients  of  the  denominator  quartic 

(iii)  Numerators  for  the  pitch  attitude  control,  which  is 
here  designated  "  CYCLIC”  . 

Actually,  the  YHC-1A  control  is  predominantly  differential  col¬ 
lective,  though  some  cyclic  variation  is  included  as  shown  by  the  non¬ 
zero  X&  at  hover.  For  each  numerator  the  roots,  D.C.  gain  (i.e.,  Bode 
gain ),  and  root  locus  gain  are  presented.  In  addition,  various 
secondary  items  of  information  such  as  the  coefficients  of  the  numera¬ 
tor  polynomial  are  printed  out.  The  numerators  correspond  to  0,  u.  w, 

a  (as  measured  by  an  accelerometer,  a  *  w  -  U  q  +  g  sin  Y  /q  dt),  and 
z  zoo 


* 

The  equations  of  motion  are  the  standard  longitudinal  and  lateral 
sets.  These  can  be  obtained  from  the  coupled  equations  given  in 
Chapter  VII  by  neglecting  all  cross- coupling  terms. 
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horizontal  velocity  perturbations  x  {■  u  for  y0  ■  0).o  In  addition,  the 
print-out  includes  rate-of- climb  perturbations  (  l?  -  h0)  and  perturba¬ 
tions  in  velocities  parallel  and  normal  to  the  undisturbed  flight  path. 
Note  that  only  for  y0  «  0  is  the  normal  velocity  perturbation  equal  to 
the  rate -of- climb  perturbation. 

For  steep  approaches,  the  "normal  velocity"  perturbation  indicates 
the  rate  of  deviation  from  the  desired  flight  path.  The  "parallel 
velocity"  perturbation  indicates  the  rate  of  deviation  from  the  desired 
position  along  the  desired  flight  path.  Hence,  this  quantity  is  impor¬ 
tant  for  stationkeeping,  e.g.,  in  formation  flight  or  scheduling  of 
multiple  operations. 

Longitudinal  Eigenvalues 

For  ease  of  reference  in  the  subsequent  discussion,  the  eigen¬ 
values  are  summarized  in  Table  VI  below. 


TABLE  VI.  TANDEM- ROTOR  LONGITUDINAL  DENOMINATORS 

Rate  of 
Descent 
(fpm  ) 

Horiz.* 

Speed 

(knots) 

Case 

No. 

Eigenvalues 

(rad/sec) 

0 

!  0 

L5 

-.974,  -.55,  .099  +  .45  i 

1500 

0 

l6 

-.824,  -.25,  .109  +  .70  j 

0 

60 

L5 

-2.54,  -.214  +  .545 j,  .475 

1500 

60 

l4 

-2.55,  -.55  +  .29j,  .4o6 

0 

80 

LI 

-2.50,  -.212  +  ,55j,  .45 

1060 

80 

— 

-2.57,  -.51  +  .52j,  .400 

1500 

80 

L2 

-2.56,  -.55  +  .29j,  .4o6 

—  ■  ■  ■  _  - .  .  _ 1 

Before  discussing  the  changes  that  occur  between  level  flight  and 
steep  descent,  the  nature  of  the  modes  for  level  flight  will  be  reviewed. 


This  column  shows  horizontal  speed,  not  total  speed,  and  is  therefore 
zero  for  vertical  descent. 
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For  level  flight,  the  eigenvalues  are  distributed  similarly  to 
those  of  the  HUP-1  given  in  Reference  34  (p»  258).  The  principal 
difference  is  about  a  5036  reduction  in  the  magnitude  (i.e.,  fre¬ 
quency).  This  would  be  expected  from  the  increased  mass  and  size 
of  the  YHC-1A  (see  Reference  35)  •  It  appears,  therefore, that  the 
general  conclusions  of  Reference  3 4  regarding  the  physical  nature  of 
the  modes  should  still  be  valid.  That  is,  away  from  hover,  the  motion 
consists  of  three  modes,  one  of  which  is  aperiodic  and  appreciably 
larger  than  the  others.  This  mode  is  a  degenerate  short- period  mode 
with  time  constant  T  approximated  by  the  standard  formula 

spi 


-(Mq  +  Zy  +  UQ  tyft) 


•V 


(M  n  +  V'oV 


+  UoMw-MqZw 


(107) 


This  mode  contains  little  u  and  is  dominated  by  w  and  0  motions. 

The  remaining  modes  involve  u,  w,  and  0  and  can  be  traced  to  the 
effect  of  the  negative  Mu  on  the  short-period  roots.  The  root-locus 
presentation  of  page  119  of  Reference  34  can  be  used  to  predict  the 
eigenvalues,  with  the  condition  that  -Zw  +  (Zu  <  0.  In  summary, 

the  open- loop  longitudinal  modes  in  level  flight  are  the  conventional  modes 
for  a  tandem-rotor  helicopter,  as  predicted  in  Reference  34. 

In  this  report,  our  principal  concern  is  to  detect  and  study  changes 
in  dynamic  characteristics  from  the  level  flight  condition,  so  we  shall 
now  consider  the  descent  cases;  a  more  complete  discussion  of  level  flight 
dynamics  of  tandem-rotor  helicopters  is  given  in  Reference  34. 

At  60  and  80  knots,  descent  angle  has  little  effect  on  the  denomina¬ 
tor  roots,  as  shown  by  Table  VI.  The  only  significant  change  is  an  in¬ 
crease  in  phugoid  damping.  This  would  be  expected,  since  the  phugoid 
damping  increases  with  -Xu.  Xu  depends  on  the  trimmed  drag/lift  ratio, 
which  of  course  must  increase  to  balance  out  the  X- force  gravity  com¬ 
ponent  in  descent. 

For  vertical  descent,  the  changes  are  also  rather  small.  Hence, 
if  any  marked  changes  in  tandem-rotor  control  characteristics  occur, 
when  going  from  level  flight  to  descent  at  the  same  horizontal  speed, 
the  changes  must  be  associated  with  the  transfer  function  numerators. 

These  will  now  be  discussed. 


Numerators 

At  a  given  horizontal  speed,  the  0/Be  numerator  shows  little 
change  with  rate  of  descent.  For  ease  of  reference,  the  6/Be 
numerators  are  summarized  below.  We  employ  the  convention  that  a 
positive  time  constant  is  "stable";  i.e.,  the  zero  is  in  the  left- 
half  complex  plane* 


TABLE 

VII.  TAN 

DEM- ROTOR  6/&e 

liOii 

Rate  of 
Descent 
(fp») 

Horiz. 

Speed 

(knots) 

Case 

No. 

Root- Locus 
Gain 

1/Tq. 

(rad/sec) 

1/T8g 

(rad/.ec) 

0 

0 

L5 

.555 

+.019 

+.5789 

1500 

0 

L6 

.5^56 

+.0156 

+  .267 

0 

6o 

L5 

>5 

+.055 

+  .817 

1500 

6o 

l4 

M 

+.0512 

+  .912 

0 

80 

Ll 

.48 

+.045 

+.955  | 

1500 

80 

L2 

_ 

.50 

_ 

+.058 

_ 

+.99  j 

_ 1 

amaapsapaai 

■aMBi 

The  effect  of  descent  on  l/T0p  1/T02  is  slight.  As  noted  above, 
the  changes  in  the  denominator  roots  are  also  small.  Hence,  any  effect 
of  descent  on  pitch  control  must  be  associated  with  feedbacks  other 
than  0.  One  of  the  most  common  feedbacks  for  tandem-rotor  configura¬ 
tions  is  u—6e.  This  is  employed  to  correct  the  negative  Mu  charac¬ 
teristic  of  tandem-rotor  helicopters  at  forward  speed.  Accordingly, 
the  u/6e  transfer  functions  are  examined  below. 

Volume  III  presents  u  numerators  in  stability  axes,  and  also  in 
axes  parallel  to  the  stability  axes,  but  located  17  feet  farther  forward, 
at  the  cockpit.  The  latter  axis  system  is  appropriate  when  considering 
feedbacks  from  air  data,  since  such  data  are  measured  well  forward  on  the 
fuselage  to  minimize  rotor  downwash  effects.  For  inertial  systems  using 
stable  platforms  or  frequently  updated  corrections  obtained  from  ground 
tracking  equipment  applied  to  strapdown components,  "stability  axes”  with 
origin  at  the  c.g.  appear  to  be  more  appropriate.  Thus,  for  S.A.S.  design 
and  such  pilot  tasks  as  stationkeeping,  the  "cockpit"  axes  are  appro¬ 
priate,  whereas  for  inertial  guidance  .the  usual  "stability"  axes  are 
relevant.  As  shown  below,  the  two  se^s  of  numerators  show  marked  dif¬ 
ferences. 


115 


First,  consider  the  u  -*&e  stability  axis  numerators  tabulated 
below.  Root  locations  are  presented  here,  not  inverse  time  con¬ 
stants;  i.e.,  a  negative  real  part  of  the  root  location  indicates  a 
"stable"  root. 


TABI 

jE  VIII.  TA 

NDEM-R0T0I 

t  u/&§  NUMERATOR 

.8  IN  STABILITY  AXES 

Rate  of 

Horiz . 

Case 

Root- Locus 

Descent 

Speed 

No. 

Gain 

Root  Locations 

(fm) 

(knots) 

(rad/sec) 

0 

80 

LI 

.185 

-9.99,  -1.14,  6.84 

1500 

80 

L2 

.312 

-10.55,  -.656,  7.45 

0 

6o 

L5 

.189 

-10.2,  -.87,  7.07 

1500 

6o 

l4 

.37 

-9.08,  -.567,  6.85 

0 

0 

L5 

00 

r- 

H 

• 

-9.05,  -.37,  7.118 

1500 

0 

l6 

,13'1 

-1.045,  6.68,  14.59 

~  : i 

At  60  and  80  knots,  there  is  little  change  in  the  numerator  with 
descent  angle.  Hie  vertical  descent  case  shows  a  marked  change;  however, 
u/6e  referred  to  stability  axes  is  unimportant  for  control  in  vertical 
descent,  since  u  is  then  vertical,  and  is  more  logically  controlled  by 
collective,  rather  than  by  5  ,  which  is  the  pitch  attitude  control. 

C 

Thus  we  conclude  that,  for  inertial  feedbacks  referred  to  the  c.g. 
and  to  the  aircraft’s  desired  flight  path,  there  is  little  effect  of 
descent  on  the  vehicle  open-loop  or  closed-loop  dynamics.  By  contrast, 
for  feedbacks  using  air  data,  or  other  data  referred  to  the  pilot's 
location,  ahead  of  the  c.g.,  the  effect  of  descent  angle  is  considerable, 
as  indicated  overleaf,  in  Table  IX. 
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IAHLEIX.  lAHPai-ROTOR  u/6e  IPMPUT0R8  II  COCKPIT  AH8 


Rate  of 
Descent 
(fpm) 

Horiz. 
Speed 
;  (knots) 

Case 

No. 

Root-Locus 

Gain 

Root  Locations 
(rad/sec) 

0 

60 

LI 

.164 

-1.71,  3.22,  16.1 

1500 

80 

L 2 

.146 

-.796  +.639J,  105.9 

0 

6) 

L3 

.142 

-2.53,  -1.3^>  24.8 

1500 

60 

L4 

.126 

-.789  +.53J,  122.1 

0 

0 

L5 

.177 

-9.05,  -0.37,  7.11 

1500 

0 

l6 

.166 

-1.32,  -0.285  ,  50.36 

Some  significant  changes  occur  in  going  from  level  flight  to  descent. 
At  80  knots,  the  right-half  plane  zero  at  3.22  rad/sec  joins  with  the  left- 
half  plane  zero  at  1*71  and  forms  a  complex  pair  in  the  left-half  plane. 
The  third  zero  moves  into  the  distant  right  half- plane,  becoming  too  far 
out  to  be  significant.  Similar  changes  occur  at  60  knots.  These  changes 
are  potentially  significant,  because  there  are  no  near-cancellations 
between  numerator  and  denominator  in  the  transfer  functions,  listed  below 
for  ease  of  reference. 


u/5e  Transfer  Functions  at  80  Knots  Horizontal  Speed 


Level  Flight: 

u  /  .164  (s  ♦  1.71)(s-3.22)(s  -  16.1) 

%  '  cocKplt  '  (St  8.50JU+  .as+  .}}))  (s  -  .4?) 

1500  fpm  Descent: 

u  .l46(s  ♦  .796  1  .639  .1 ) ( s  —  105.9) 

/  cockpit  “  (s  ♦  2.36)(  s  ♦"»  i  .29  J)(«  -  -*06) 


117 


The  rather  small  l/Te1  (see  Table  VII )  limits  the  degree  to 
which  the  unstable  eigenvalue  can  be  stabilized  by  6  and  q  feedback. 
Hence  it  is  quite  possible  that  the  u  feedback  will  be  required  to 
achieve  stability,  even  if  6  and  q  fe<*d  backs  are  also  employed.  In 
these  circumstances,  the  changes  in  the  u/6e  transfer  function  should 
be  considered  carefully  when  designing  stability  augmentation  systems 
using  air  data  recorded  away  from  the  c.g. 

Let  us  now  consider  the  effects  of  descent  on  collective  pitch 
control.  The  primary  function  of  collective  pitch  is  to  control  the 
deviation  of  the  aircraft  c.g.  normal  to  the  desired  flight  path.  The 
appropriate  transfer  function  numerators  are  summarized  below. 


TABU  X.  TANDEM- ROTOR  COLLECTIVE  PITCH  NUMERATORS  FOR  VELOCITY 

PERTURBATIONS  OF  C.  G.  NORMAL  TO  DESIRED  FLIGHT  RATH 

Rate  of 
Descent 
(fjm) 

Speed 

(knots) 

Case 

No. 

Root- Locus 
Gain 

Root  Locations 
(rad/sec) 

0 

80 

U 

-9.51 

-1.87,  -.209,  >95 

1500 

80 

1*2 

-8-99 

-1.15,  -.778,  .2148 

0 

60 

L5 

-8.56 

-1.84,  -.217,  .548 

1500 

60 

1*4 

-7.85 

-.952  ♦  .5755  j,  .267 

0 

0 

L5 

-7.*»J 

-.955,  .105  *  .445j 

0 

1500 

1*6 

-1.466 

-.652,  .054  +  1.26J 

At  hover,  the  usual  cancellation  of  the  denominator  roots  occurs, 
leaving  the  transfer  function  as  fi/6c  ■  Z5  /[s  -  z^).  For  vertical 
descent,  the  above  transfer  function  becomes  x/lc  and  is  not  important 
for  practical  flight  control.  Chapter  VII  shows  that  for  single-rotor 
helicopters*  control  of  velocity  normal  to  the  unperturbed  flight  by 
collective  pitch  becomes  difficult  at  very  low  speeds  (i  20  knots)  in 
steep  descents.  Unfortunately,  derivatives  were  not  available  for  the 
YHC-1A  in  such  conditions,  so  it  is  not  known  whether  this  effect  (due 
to  an  uncancelled  right-half  plane  zero)  also  occurs  on  tandem- rotor 
helicopters . 
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The  importance  of  the  changes  in  the  above  numerators  for 
descending  flight  is  hard  to  assess  in  general,  because  the  ”  normal 
velocity”  collective  control  is  affected  by  feedbacks  to  cyclic  from 
u,  q,  and  0.  Such  feedbacks  will  probably  be  employed  for  inner-loop 
stabilization. 

It  is  not  feasible  to  discuss  all  the  tandem-rotor  longitudinal 
transfer  functions  that  are  tabulated  in  Volume  III.  However,  from  those 
that  have  been  examined  here,  it  would  appear  that  descent  angle  does 
not  cause  radical  changes  in  the  transfer  functions.  The  most  important 
exception  to  this  generalization  is  u/&  /  cockpit.  The  numerator  of  this 
transfer  function  exhibits  seme  significant  changes ,  which  should  be 
considered  when  designing  stability  augmenter  systems  that  employ  M 
augmentation  based  on  air  data  sensors  located  away  from  the  c.g. 


VI.  3  pysMTcs 

The  lateral  transfer  function  denominators  for  the  YHC-1A  are 
presented  in  Table  XI.  The  table  has  been  arranged  to  fit  the  following 
discussion  which  first  considers  the  variation  of  the  level  flig.it 
denominators  with  forward  speed,  and  then  analyzes  the  effects  of 
descent. 


XABI£  XI.  TANDEM- ROTOR  LATERAL  EIGENVAIAJES 

Rate  of 
Descent 
(fpa) 

Case 

No. 

Eigenvalues 

(rad/sec) 

0 

0 

L5 

.0594  +  .5113 

-.057 

-.975 

0 

6o 

L3 

.225  1  .6043 

CO 

$ 

• 

1 

-1.29 

0 

8o 

U 

.281  +  .6423 

-.0768 

-1.36  ' 

1500 

0 

L6 

.052  ♦  .7073 

-.052 

-.925 

1500 

6o 

l4 

.134 

.542 

-.048 

-1.77 

1500 

80 

L2 

-.140 

.772 

.0023 

-1.832 
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The  hover  behavior  is  standard  and  closely  analogous  to  the 
longitudinal  case  illustrated  in  Table  VI .  Note  the  close  proximity 

of  the  frequencies  of  the  longitudinal  hovering  phugoid  (to  =  o.  U5 
rad/sec)  and  the  lateral  hovering  dutch  roll  («*>  =0.51  ra8/sec). 

As  stated  in  Chapter  V,  for  the  tandem  configuration,  longitudinal- 
lateral  cross-coupling  was  arbitrarily  neglected.  Uncbubtedly,  some 
coupling  exists,  and  the  near -coincidence  of  to  and  to  may  make  this 
coupling  important  for  the  YHC-1A.  The  quest  ion1, will  oRly  be  settled 
by  more  extensive  measurements  of  derivatives  and  flight  test  responses 
to  control  inputs. 

The  principal  effect  of  descent  at  the  higher  speeds,  is  to 
change  the  unstable  dutch  roll  oscillatory  roots  into  an  aperiodic 
pair.  At  60  knots,  both  these  roots  are  unstable }  but  as  speed  is 
increased  one  becomes  stable  while  the  other  root  moves  farther 
out  into  the  right-half  plane.  For  the  YHC-1A,  this  degeneration  of 
the  dutch  roll  roots  into  an  aperiodic  pair  is  a  consequence  of  the 
low  (actually  negative)  Nv« 


N  remains  negative  for  all  the  level  and  descending  forward 
flight  conditions  considered  here,  when  measured  in  stability  axes 
(which  is  the  appropriate  axis  system  for  defining  directional  sta¬ 
bility)  .  It  is  difficult  to  design  a  tandem-rotor  configuration  with 
the  c.g.  far  removed  from  the  mid-point  of  the  rotor  axes;  hence, 
negative  N  is  characteristic  of  these  configurations.  It  is 
usually  "fYxed"  by  feeding  back  p  to  6_  to  achieve  static  stability, 
plus  additional  feedbacks  of  p  to  d.,  and  possibly  r  to  6R  command* 
to  move  all  the  roots  into  the  right-half  plane  and  to  achieve  good 
turn  coordination. 

Feedback  of  sideslip  alone  is  usually  insufficient  to  ensure 
stability.  This  is  illustrated  in  Figure  57  for  the  80-knot  level- 
flight  case.  Note  the  awkwardly-placed  zero  of  the  sideslip 
numerator.  Because  of  the  fundamental  importance  of  directional 
stability  for  good  handling  qualities,  particular  interest  attaches 
to  the  variation  of  the  transfer  functions  with  descent  angle. 

Numerators  of  these  transfer  functions  are  tabulated  below. 
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Figure  37 •  Hoot  Locus  for  v  — *  5  Feedback  on 
Tandem- Rotor  Helicopter. 
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TABLE  XII .  TANDEM-ROTOR  v/&R  NUMERATORS 

Rate  of 
Descent 
(fpm) 

Horiz . 
Speed 
(knots ) 

Case 

No. 

Root 

-Locus 

Gain 

Root  Locations 

(rad/sec) 

0 

80 

LI 

.099 

-  .92,  -.00145, 

256.9 

1500 

80 

L2 

.07 

-  .41,  -.020, 

422.9 

!  0 

60 

L3 

.106 

-  .86,  -.0025, 

163.5 

1500 

60 

l4 

.073 

-  .151, -.0239, 

524.9 

0 

0 

L5 

.146 

-  .0226, .79  +  6.04j 

1500 

0 

L6 

.143 

-1.33,  -.055, 

77.96 

In  these  numerators,  v  is  referred  to  stability  axes .  For 
analyzing  stability  augmenter  systems,  it  is  perhaps  more  meaningful 
to  refer  v  to  a  location  near  the  cockpit,  where  the  sideslip 
sensing  ports  are  located.  However,  as  shown  in  Volume  III  , 
there  is  relatively  little  difference  between  the  v/d  transfer 
functions  referred  to  either  location  at  60  and  80  knots.  The  v/d 
transfer  function  is  not  particularly  important  at  zero  forward 
speed,  since  v  is  then  controlled  directly  through  lateral  cyclic 
pitch. 


From  Table  XII,  the  effect  of  the  descent  on  the  v/6  numerators 
is  seen  to  be  relatively  minor  at  60  and  80  knots.  The  low-frequency 
zero  limits  the  possible  improvement  of  the  unstable  spiral  root  at 
80  knots,  suggesting  that  additional  feedbacks  (e.g.,  {5  — *  5^)  are 
required  to  stabilize  this  root.  The  degenerate  dutch  roil  modes 
may  also  require  excessive  gain  for  stabilization  through  v— 
feedback  alone.  (The  low  accuracy  to  which  v  can  be  measured  constrains 
practical  v  -+  6  feedbacks  to  rather  low  gains. )  3he  net  effect  of 
descent  upon  v  d  feedback  control  seems  to  be  adverse,  and  it 
appears  that  a  system  designed  for  level  flight  would  not  necessarily 
be  suitable  for  steep  approaches. 
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From  the  above  discussion,  v  -»  0  feedback  appears  to  be  manda- 
tory  for  tandem-rotor  helicopters.  Suen  feedback  represents  an  "inner 
loop"  which  will  affect  the  other  loop  closures  such  as  #  —  $  .  In 
level  flight,  the  */\  numerator  exhibits  some  undesirable  non¬ 
minimum  phase  characteristics*  as  noted  in  Reference  34.  These  become 
appreciably  worsened  in  descent,  as  shown  in  Table  XIII. 


TABLE  XIII 

TANDEM  ROTOR  4/&A  NUMERATORS 

1  ^"1 

Rate  of 
Descent 
(fpm) 

Horiz. 

Speed 

(knots) 

Case 

No. 

Root 

-Locus 

Gain 

Root  Locations 
(rad/sec) 

0 

80 

LI 

.495 

-  .70, 

.543 

1 

i 

i 

1500 

80 

L2 

.48 

-1.09, 

.911 

i 

1 

0 

60 

L3 

.51 

-  .496, 

.368 

! 

1 

1500 

60 

L4 

.45 

-  .856, 

.70 

! 

i 

0 

0 

L5 

.55 

• 

. 

ro 

s 

* 

a 

.057 

1500 

0 

l6 

.52 

-  .055  ♦ 

•02j 

i 

=H 

The  poor  handling  qualities  associated  with  the  unaugmented 
p/ 6  closure  in  level  flight  noted  in  Reference  34  are  likely  to 
become  even  more  degraded  in  descent.  For  example,  at  80  knots, 
the  #/0.  transfer  fraction  has  one  right-half  plane  zero  and  two 
right-half  plane  real  poles,  thus  presenting  a  situation  which  is 
quite  difficult  to  stabilize. 

An  important  feedback  for  human-pilot  control  is  y^  0., 
where  y  is  the  lateral  displacement  of  the  helicopter  measured 
either  at  the  c.g.  or  at  the  cockpit. 


*  A  nonminimum  phase  transfer  function  is  one  with  one  or  more 
right-half  plane  zeros. 
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The  print-out  in  Volume  II  includes  both  alternatives,  however 
at  the  higher  speeds  the  differences  between  the  transfer  functions 
for  y/^/cockpit  811(1  y/*  /  eg.  81,0  relatively  minor.  For  vertical 

descent  the  differences  become  more  pronounced,  as  shown  by  the 
following  comparison. 


TABLE  XIV 

.  EFFECT  OF  VERTICAL  DESCSIT  ON  TANDB4- ROTOR 
y/b.  NUMERATORS  MEASURED  AT  THE  CHtTER  OF 

GRAVITY  AMD  AT  THE  COCKPIT 

Flight 

Condition 

KCMi 

Location 
of  Origin 

Root 

-Locus 

Gain 

Root  Location 
(rad/sec) 

L5 

c.g. 

.998 

-.082  ♦  4.22j,  -.037,  -.0177 

Hover 

L5 

cockpit 

3.26 

-.108  ♦  2.33J,  -.053,  -.0123 

Vertical 

Descent 

L6 

c.g. 

.98 

-.083  ♦  U.05j,  -.051,  0 

Vertical 

Descent 

L5 

cockpit 

-.0116 

-26.7,  -.083,  0,  52.8 

A  complete  assessment  of  the  significance  of  this  difference  in 
y/6.  is  beyond  the  scope  of  this  report.  It  would  involve  loop 
closures  using  assumed  models  for  the  pilot;  e.g.  the  "series  closure" 
model  used  in  Reference  36  and  in  the  PAPER  PILOT  program  described  in 
Reference  37*  The  validity  of  these  closures  for  conditions  other 
than  hover  is  subject  to  question.  However,  if  the  series  closure  is 
employed,  the  y/d./fVV,k_it  and  #/6  transfer  functions  would  be  relevant 
for  the  inner  ana'  outer1ioops .  As  shown  in  Tables  XIII  and  XIV,  marked 
changes  in  these  transfer  functions  occur  between  hover  and  vertical 
descent.  It  is  therefore  to  be  expected  that  correspondingly  large 
changes  in  handling  qualities  will  appear  when  going  from  hover  to 
vertical  descent. 

8ummary 

In  this  chapter,  we  have  reviewed  the  major  effects  of  descent 
on  the  transfer  functions  of  a  typical  tandem-rotor  helicopter.  It 
was  assumed  that  longitudinal  and  lateral  motions  were  uncoupled. 
Longitudinally,  the  effects  of  descent  (at  normal  approach  speeds  of 
60  to  60  knots)  are  generally  minor.  Some  increase  in  phugoid  damping 
occurs,  but  the  important  Q/be  pitch  attitude  control  numerator 
hardly  changes.  Similarly,  little  change  was  observed  in  the  u /bc 
numerator  referred  to  stability  axes. 
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However,  narked  changes  occur  in  the  u/^  numerator  referred  to 
body  axes  with  origin  at  the  cockpit.  These  changes  are  potentially 
significant  for  stability  augmenter  systems  which  employ  air  data  sensors 
to  augment  ^ 

Laterally,  the  effects  of  descent  are  more  noticeable.  At 
forward  speed  the  unstable  dutch  roll  degenerates  into  an  aperiodic 
pair,  having  at  least  one  unstable  root.  This  unstable  root  is  usually 

stabilised  by  v  bp  feedback j  however,  this  feedback  has  little 
effect  on  the  spiral  mode  (which  nay  be  unstable)  because  of  a 
M close- in"  zero  in  the  v/6g  transfer  function.  This  awkwardly 
placed  zero  is  retained  in  descending  flight.  Thus,  additional  feed¬ 
backs  (e.g.,  i  6a,  r  5p)  are  still  required  to  achieve  good  handling 
qualities.  The  nonniniaua  phase  characteristics  of  the  0/6* 
numerator  worsen  in  descent,  thus  increasing  the  need  for  still 
further  feedbacks,  such  as  r  bp.  Because  of  the  large  msnber 
of  feedbacks  required*  it  is  hard  to  fora  an  assessment  of  the 
effects  of  descent  on  closed-loop  characteristics  without  per¬ 
forming  detailed  calculations  on  *ious  alternative  loop  closures. 
However,  the  changes  in  the  transfer  functions  are  sufficiently 
large  and  numerous  to  indicate  the  necessity  of  performing  these 
calculations  for  descent  in  addition  to  level  flight. 
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C  HPH1  VU.  SINGLE-ROTOR  HELICOPTER  DYNAMICS 


VII.1.  INTRODUCTION 


Thia  chapter  diacusaea  the  open- loop  and  closed- loop 
dynamlca  of  a  typical  single-rotor  helicopter,  the  Sikorsky  S-58. 

The  baaic  data  are  preaented  in  Volume  III,  which  tabulates  deriva¬ 
tives  and  eigenvalues. 

The  flight  conditions  preaented  in  Volume  III  inclule 
level  flight  and  rates  of  descent  of  7*5  fpe,  15  fps,  and  22.5  fps. 
At  loir  speeds  the  latter  condition  is  virtually  at  the  boundary 
of  the  vortex-ring  state.  For  each  rate  of  descent^  true  airspeeds 
of  0,  20,  1*0,  60  and  100  knots  are  considered. 

All  lateral-longitudinal  cross-coupling  derivatives  were 
included,  as  shown  by  the  equations  of  motion  presented  in  Figure  58. 
To  illustrate  the  importance  of  cross-coupling,  one  case  was  re-run 
with  all  cross-coupling  derivatives  set  to  zero.  This  was  done 
for  level  flight  at  100  knots. 

The  cyclically  varying  tip  loss  factor  was  included  in  all 
cases.  To  demonstrate  the  effect  of  tip  losses  on  open- loop 
dynamics,  two  cases  were  re-run  with  tip  losses  removed.  These 
cases  were  22.5  fps  rate  of  descent  with  40  knots  airspeed,  and 
vertical  descent  at  22.5  fps. 

Eigenvectors  foar  u,  v,  w,  p,  q,  r,  in  stability  axes, 
residues  for  unit  impulse  inputs  in  A1g,  B1#,  0Q,  6qTR,  and 

transfer  functions  relating  u,  v,  v,  p,  q,  r,  to  A1g,  B1g, 

0 1  ,  6  were  also  calculated  for  all  the  above  cases.  To 
save  I  toe,  these  data  are  not  presented  in  full  in  Volume  III, 
but  selected  values  are  given  in  this  chapter  where  appropriate. 

The  main  points  discussed  here  are: 

(i)  The  effects  of  forward  speed,  rate  of  descent,  and 
tip  losses  on  the  eigenvalues 

(il)  Cross-coupling  effects  cu  open- loop  dynamics  as 
demonstrated  by  the  eigenvectors 

(Hi)  Cross-coupling  effects  on  c lea ed- loop  dynamics  and 
huean  pilot  control  as  demonstrated  by  the  transfer 
functions  and  residues 
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The  volume  of  data  generated  for  the  S-58  is  large,  and  it  is 
not  easy  to  condense  it  to  its  essentials.  Only  the  "highlights" 
can  be  indicated  here.  Qhe  reader  will  gain  worthwhile  insight  into 
helicopter  dynamics  by  examining  Volume  III  in  detail. 
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where 


[B]  |x|  +  [c]{x|  =  [F] 

T 

{x}  =  [  uwqOvp0r] 


Figure  38.  Aircraft  Equations  of  Motion  in  Stability  Axes 
for  Small  Perturbations  from  Straight-Line 
Flight  with  Zero  Bank  Angle. 
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VII.  2.  KTflBWVATJTBa 


The  eigenvalues  are  summarized  in  Table  XV  and  are  discussed 

below. 


The  effects  of  cross -coupling  on  most  eigenvalues  are  not 
large.  For  example,  at  100  knots  in  level  flight,  only  two  modes 
are  significantly  changed  when  the  cross-coupling  derivatives  are 
suppressed.  This  is  illustrated  by  the  excerpt  from  Table  XI 
reproduced  below. 

S-58,  100  Knots  Level  Flight 


Cross-Coupling  Eigenvalues 


Yes 

-.46  +  1.63J, 

-4.96, 

-.12, 

-.47, 

-1.82, 

.12  +  .33d 

No 

-.48  +  1.73d, 

-5.23, 

-.11, 

-.32, 

-1.79, 

.22  +  .34d 

As  will  be  shown  later,  the  fact  that  most  eigenvalues  are 
essentially  unchanged  by  the  cross-coupling  derivatives  does  not  mean 
that  cross-coupling  effects  are  unimportant. 

For  level  flight,  the  eigenvalues  change  only  slightly  as 
forward  speed  is  increased.  The  unstable  phugoid  mode 
and  the  degenerate  short-period  modes  characteristic  of  hover  remain 
remarkably  unchanged  with  forward  speed.  This  is  contrary  to  the 
usual  trend  noted  in  Reference  34,  where  the  phugoid  becomes  stable 
and  the  l/T  ,,  l/T  q  modes  merge  to  form  the  conventional  short- 
period.  sp1  sp^ 


The  roll  subsidence  root  remains  roughly  equal  to  Lp  over  the 
speed  range  0  to  100  knots.  The  dutch  roll  mode  is  stable  even  at 
hover  due  to  the  N  of  the  tail  rotor  as  explained  in  Reference 
38,  page  15).  The  frequency  increases  from  cd  =  0.73  to  co  =  I.I63 
rad/sec,  and  the  damping  ratio  rises  slightly  from  0.098  to  0.131. 

In  descent,  at  the  higher  speeds  there  is  not  much  change 
in  the  eigenvalues.  By  analogy  to  fixed-wing  aircraft,  in  which  the 
phugoid  damping  ratio  is  roughly  proportional  to(D/L),one  would 
expect  that  o>p  would  increase  in  descent.  This  does  occur,  but  the 
effect  is  not  large.  Thus  at  100  knots  in  level  flight, £  =  -.172, 
whereas  at  the  same  horizontal  speed  and  22.5  fps  rate  ofpdescent, 

£  =  -.14.  At  lower  speeds  a  similar  effect  occurs.  For  20  knots 

fo?ward  speed  in  level  flight, £  =  -.111;  with  the  same  horizontal 

speed  and  22.5  fps  rate  of  desclnt,  £  =  -.057* 

P 

The  largest  effect  of  descent  on  the  eigenvalues  occurs  at  very  low 
forward  speeds.  In  addition  to  the  increased  phugoid  damping  noted 
above,  the  degenerate  short-period  roots  coalesce  to  form  a  heavily 
damped  oscillatory  pair. 
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TABLE  XV.  SIKORSKY  S-^8  EIGENVALUES  WITH  SIX  DEGREES  OF  FREEDOM 
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In  summary ,  it  appears  that  the  variation  of  the  eigenvalues 
both  with  horizontal  speed  and  rate  of  descent  is  less  than  one  might 
have  expected.  In  part,  this  is  due  to  the  small  and  rather  low 
aspect  ratio  tail  surfaces  of  the  S-58«  These  are  extremely  ineffective, 
in  that  remains  positive  throughout  the  speed  range  0-100  knots.  A 
more  effective  tail  would  tend  to  produce  "airplane"  type  modes,  i.e., 
oscillatory  short-period  and  stable  phugoid,  as  speed  is  increased. 

VII .3*  CROSS-COUPLING  EFFECTS 


As  noted  above,  only  minor  changes  in  the  eigenvalues  occurred 
when  the  cross-coupling  derivatives  were  set  to  zero.  This  trend  is 
shown  on  Table  XV  for  100  knots  level  flight,  but  has  also  been  verified 
by  some  calculations  (not  published  here)  for  descending  flight  and 
lower  speeds  with  and  without  the  coupling  terms.  From  this,  it  might 
be  inferred  that  cross -coupling  effects  are  minor.  Nothing  could  be 
further  from  the  truth!  Examination  of  the  print-outs  (Volume  III) 
reveals  strong  cross- coupling  at  all  flight  conditions. 

In  this  section  we  demonstrate  lateral- longitudinal  cross¬ 
coupling  by  three  complementary  approaches 

1.  Calculation  of  transfer  functions  with  and  without  the 
cross -coupling  derivatives 

2.  Calculation  of  the  residues  of  the  partial  fraction 
expansion  of  the  Laplace  transform  of  the  time  response 
to  a  given  input,  such  as  a  step  or  impulse  in  Blg,  or  0q. 

3.  Study  of  the  eigenvectors  (the  ratios  of  which  equal  the 
ratios  of  the  residues  in  each  degree  of  freedom  for  a 
given  eigenvalue). 

We  show  that  the  cross -coupling  between  longitudinal  and 
lateral  motions  is  strong  for  many  of  the  modes.  Furthermore,  the 
cross -coupling  occurs  at  all  forward  speeds  and  rates  of  descent 
examined  here.  In  this  respect  the  S-58  is  believed  to  be  repre¬ 
sentative  of  single-rotor  helicopters  in  general.  Little  change 
occurs  in  the  lateral- longitudinal  cross -coupling  for  each  mode  as 
descent  angle  is  increased  at  a  given  forward  speed.  Hence,  the 
explanation  of  the  difficulties  encountered  in  flight  path  control  in 
steep  descent  must  be  elsewhere.  In  the  final  section  of  this 
chapter,  the  culprit  is  identified  as  a  nonminimum  phase  y/oo 
transfer  function.  This  transfer  function  is  crucial  in  flight  path 
control  at  low  speeds,  and  it  is  shown  that  in  descent  a  nonminimum 
phase  zero  appears  which  causes  a  marked  deterioration  in  the  achievable 
level  of  accuracy  with  which  the  desired  approach  path  can  be  followed. 

#  " 

A  nonminimum  phase  transfer  function  has  one  or  more  right-half 
plane  zeros.  The  terms  'nonrainimum  phase  zero'  and  'right-half 
plane  zero'  are  synonymous. 
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Effect  of  Cron-Coupling  Derivative!  on  the  Trtngftr  Fungtlong 


For  brevity,  transfer  functions  will  not  be  written  in  the 

usual  form;  e.g. , 

’  r  =  n 

Output (s)  (Root  locus  gain)  r  =  1 

Input (s)  -  i  =  m 

n 

i  =  1 


(s  +  z  ) 
r 


(s  +  Pi) 


Instead,  we 

Output 

Input 


adopt  the  equivalent  notation 

_  r  =  n  /  \ 

(Root  locus  gain)  r  =  1  zr 


n 


i  =  m 
i  =  1 


(-Pi) 


For  example,  the  transfer  function 


w(s 

6 

o 


154.4  (s  +  0.802)  (s  -  0.32)  _ _ 

(s  +  0.312) (s  -  0.035^  +  0.34j)(s  -  0.0354  -  0.34j ) 


is  written  as 

w  _  1^4.4  (-0.802)  (  +  O.32) 

0o  “  (-0.312)  (+  0.0354  +  o.34j) 

In  this  abbreviated  style,  right-half -plane  poles  and  zeros 
are  indicated  by  numbers  with  positive  real  parts. 

As  noted  in  Table  XV,  the  100- knot,  level- flight  case  was  run 
twice.  Firstly,  all  the  derivatives  calculated  from  MOSTAB  were 
included.  Secondly,  all  derivatives  coupling  longitudinal  forces 
or  moments  with  lateral  perturbations  and  vice  versa  were  dropped. 

In  addition,  all  lateral  forces  and  moments  produced  by  0  and  Blg 
were  put  equal  to  zero,  as  were  all  longitudinal  forces  and  moments 
produced  by  0qTR  andA^. 

The  results  are  compared  in  Table  XVI. 
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TABLE  XVI.  EFFECT  OF  CROSS- COUPLING  DERIVATIVES; 

SIKORSKY  S-58,  100  KNOTS,  LEVEL  FLIGHT 


Cross  Pitch  Rate/BL  Transfer  Function 

Coupling  'Ls 

No  8. 43 (-0.657) (-0.029)  (0) 

(-0.33) (-1.79) (0.219  +  0.336 j ) 

Yes  8,48(0.017)  (-O.061) (-0.109)  (-5.25)  (-0.46  -  1.693 ) (-0.66) 

(0.117) (-4.96) (-0.4o6) (0.12  +  0.33d) (-0.46  +  1.633 ) (-1.82) 


With  approximate  pole-zero  cancellations  removed,  this  becomes: 


Yes  8.48(0.017)  (-0.061)  (-0.66) 

(-0.466) (-1.82) (0.12  +  0.33d) 


Cross  Roll  Rate/A^  t  Transfer  Function 

Coupling  ■L# 

No  -59.0(0.0000017)  (-0.60  +  1.643) 

(-0.11) (-5.23) (-0.476  +  1.73d; 

Yes  -58.9(0.000033) (-0.58  +  1.69.1) (0.234  ±  0.328,1) (-1.86) (-0.31) 

(-0.117) (-4. 96)  (-0.466) (0.12  +  0.333K-0.46  +  1.633 )(-l. 82) 

With  approximate  pole-zero  cancellations  removed,  this  becomes: 
-58.9(0)(-0.577  +  1.69j  )(-:0. 311) 

(-0.117) (-4.96) (-0.46  +  1.63d)  (-0.466) 


Yes 


It  is  not  easy  to  assess  the  relative  magnitude  of  the  pitch 
versus  the  roll  motions  produced  by  B,  by  comparing  the  transfer 
functions  as  listed  in  Tables  XVI  ana  XVII.  A  more  direct  comparison  is 
obtained  by  expanding  the  q/B,  g  and  p/B^  transfer  functions  into 
partial  fractions  and  comparing  the  magnitude  of  each  coefficient  in 
the  expans iorv  i.e.,  each  residue,  as  shown  in  the  lower  half  of 
Table  XVII.  Table  XVIII  shows  such  a  comparison.  For  compactness,  some 
numbers  have  been  rounded  off  in  Table  XVIII,  and  only  the  residue 
corresponding  to  the  upper-half  plane  root  of  each  complex  pair  is 
shown. 
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TABLE  XVIII.  CROSS-COUPLING  EFFECTS  ON  B.  IMPUISE  INPUT,* 

SIKORSKY  S-p8  AT  100  KNOTS,  §EA  LEVEL 

Level  Flight 

- 1 

7.65  Deg.  Descent  j 

Root 

Residues 

Root 

Resj 

Ldues 

s,  rad/sec 

P 

<1 

s,  rad/sec 

P 

q. 

-  .117 

.061 

.0033 

-  .0926 

-  .0062 

-  .797 

-4.96 

-6.99 

.7245 

-5.25 

-5.59 

.469 

-  .466 

4.91 

-  .622 

-  .421 

6.l6 

-  .797 

.121 

1.82 

1.2 

.091 

2.39 

1.18 

+  .33d 

/  24° 

Z-1520 

+  -  31J 

A  196 

t-i47° 

-  .46 

1.63 

.264 

-  .541 

3.14 

.511 

+  1.63J 

/-  97° 

L  24° 

+  1.54J 

1.  107° 

^.29°  i 

-1.824 

-  .955 

-6.942 

-I.865 

-2.37 

-6.851  ; 

Table  XVII  shows  p  and  q  residues  for  inputs  at  two  flight 
conditions.  To  form  a  complete  picture  of  cros s -coupling, u,  v,  w,  and 
r  must  also  be  considered.  However,  Table  XVII  suffices  to  demonstrate 
that,  in  level  flight  at  100  knots,  Blg  induces  considerable  roll 
(compared  to  pitch)  and  that  this  crosl -coupling  is  not  much  changed 
by  descent. 

To  understand  Table  XVIII,  first  compare  the  relative  magnitudes  of 
the  residues  for  each  mode.  The  relative  importance  of  p  to  q  in  each 
mode  depends  upon  the  ratio  of  the  residues.  For  example.  Table  XVIII 
shows  that  the  rapidly  decaying  s  =  4.96  mode  is  predominantly  roll. 

On  the  otner  hand,  the  unstable  oscillatory  mode  displays  almost  as 
much  pitch  as  roll,  and  the  s  =  -1.824  mode  has  about  seven  times  as 
much  pitch  rate  as  roll  rate. 

The  above  remarks  indicate  the  relative  magnitude  of  pitch  and 
roll  in  each  mode,  but  this  information,  by  itself,  is  insufficient 
to  predict  the  cross -coupling  induced  by  the  pitch  control,  B^g.  To 
determine  that,  we  must  scrutinize  the  absolute  magnitude  of  each  of 
the  q  residues.  If  the  residues  of  the  "pitch"  modes  (s  =  0.12  + 

0*33j>  s  =  -1.324)  are  very  much  larger  than  the  residues  of  the  other 
"coupled"  or  "roll"  modes,  then  the  response  will  be  dominated  by  pitch, 
and  Bls  will  induce  little  cross-coupling.  Examining  Table  XVIII  leads 
to  the  conclusion  that  although  the  residues  of  the  "pitch"  modes  are 
somewhat  larger  than  those  of  the  other  modes,  the  difference  in 
magnitude  is  not  so  overwhelming  that  we  can  say  that  the  roll  induced 
by  B]_s  is  small.  We  therefore  conclude  that  B,  induces  appreciable 
roll  in  level  flight  at  100  knots.  S 
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The  descent  situation  is  summarized  by  the  right-hand  side  of 
Table  XVIII.  Overall,  there  is  very  little  change  in  the  residues  between 
level  flight  and  descent  at  100  knots.  If  this  result  is  representative, 
it  implies  that  the  difficulties  in  control  in  steep  approaches  due  to 
cross-coupling  are  no  worse  than  those  encountered  in  level  flight. 

Hence,  we  cannot  ascribe  the  difficulties  experienced  in  following 
steep  glide  slopes  to  changes  in  cross-coupling  between  level  flight 
and  descent. 

To  verify  that  the  above  result  is  representative,  we  have 
calculated  the  residues  for  u,  v,  w,  p,  q,  r,  to  Alg,  Bls,  0q, 

0  for  a  wide  range  of  flight  conditions. 

O  J.K 

Table  XIX  illustrates  a  typical  comparison  between  level  flight 
and  two  descending  cases,  all  at  40  knots  airspeed.  This  confirms  the 
conclusions  of  Table  XVIII that  descent  angle  has  only  a  minor  influence 
on  cross -coupling.  From  examination  of  the  residues  for  the  entire 
range  of  flight  conditions  presented  here  (0  to  100  knots  airspeed, 
and  zero  to  22.5  fps  rate  of  descent),  the  effect  of  descent  angle  is 
quite  small.  We  therefore  conclude  that  variations  in  cross - 
coupling  between  level  flight  and  descent  are  not  a  major  source  of 
difficulty  in  controlling  steep  approach  paths  on  single-rotor 
helicopters . 


Before  leaving  the  subject  of  cross -coupling,  it  should  be  noted 
that,  although  the  effect  of  descent  rate  is  small,  cross-coupling  varies 
significantly  with  airspeed.  This  can  be  seen  by  examining  the  transfer 
functions  presented  in  Volume  III  .  It  is  usual  to  reduce  the 
effects  of  cross- coupling  by  suitable  phasing  of  A1  and  B..  to  Blg, 
through  offset  of  the  servos.  This  correction  can ^be  exact? 
only  at  a  single  flight  condition. 
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TABLE  XIX.  EFFECT  OF  DESCENT  OK  CROSS-COUPLING 

SIKORSKY  S-58  AT  40  KNOTS  T.A.S.  SEA  IZVEL 
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-1.33  -1.23  -1.198  -6.045  -7.10  -7.11  2.45  2.64  4.89 


VII. 4.  HCOHI  OOraOL 

The  effect  of  descent  angle  on  the  transfer  function  relating 
flight  path  angle,  y ,  to  collective  pitch  is  most  significant.  The 
difference  between  level  flight  and  steep  descent  cases  for  this 
transfer  function  is  such  that  control  at  low  speeds  is  seriously 
degraded  in  descent.  As  explained  below,  the  effect  is  caused  by  a 
right- half  plane  zero.  It  appears  that  this  zero  is  the  primary  cause 
of  the  control  difficulties  experienced  in  steep  descent. 

It  is  first  necessary  to  define  the  quantities  of  importance  in 
precise  flight  path  control  for  slov*  steep  approaches. 

The  key  quantity  determining  how  accurately  a  given  approach 
flight  path  is  followed  is,  in  general,  not  height,  nor  flight  path 
angle,  but  the  displacement  of  the  helicopter  from  the  desired  approach 
path.  This  quantity  is  illustrated  in  Figure  39*  For  brevity  we 
shall  call  it  "height  error",  denoted  by  the  symbol  n. 

This  section: 

(1)  defines  h  in  terms  of  w  and  q  (for  which  transfer  functions 
have  been  calculated) 

(2)  presents  transfer  functions  relating  h  to  collective 
pitch 

/  A 

(3)  demonstrates  that,  at  low  speeds,  control  of  h  by  collec¬ 
tive  pitch  is  relatively  easy  in  level  flight  but 
difficult  in  steep  descents 

(4)  traces  the  reason  £or  the  difficulty  to  a  nonminimum 
phase  zero  in  the  h/0Q  transfer  function 

(5)  investigates  methods  for  alleviating  the  difficulty 
through  typical  stability  augmenter  systems 

Definition  of  Height  Error,  ft 

Consider  Figure  39*  This  indicates  the  distance  of  the  helicopter 
c.g.  from  the  desired  flight  path.  Call  this  distance  the  "height  error" 


Using  stability  axis  quantities , 

I  -  u0/*dt-w 

Note  that  this  expression  neglects  products  of  small-perturbation 
quantities,  and  also  products  of  small- perturb  at  ion  quantities  with 
the  small  trim  bank  angle  which  is  typically  of  the  order  of  1  degree. 


138 


Figure  j59« 


Definition  of  Height  Error, 


A 

h. 
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The  flight  path  angle  perturbation,  with  the  same  assumptions, 

becomes 


7  = 


A 

h 


For  vertical  and  almost- vertical  descents  the  above  expressions  are 
not  appropriate.  Vertical  descents  are  usually  made  close  to  the 
ground  and  the  key  parameter  is  the  height  above  the  ground.  Using 
stability  axes,  in  vertical  descent  this  parameter  is  given  by 


h  =  J  (U  +  u)dt 


A 

Transfer  Functions  for  h  to  Collective  Pitch 

Volume  III  does  not  give  the  h  numerators  directly,  but 
these  can  be  constructed  from  the  q/s  and  w  transfer  functions.  The 
results  are  displayed  in  Table  XX. 

As  will  be  shown,  the  significant  feature  to  note  in  Table  XX  is 
the  right-half  plane  zero.  At  20  knots  (Cases  6,  11,  and  20),  this 
zero  increases  with  descent  angle.  This  trend  is  reversed  at  the  higher 
speeds,  but  this  reversal  is  not  relevant  to  height- control  because 

above  about  4o  knots  B,  becomes  more  effective  than  0  in  controlling 

is  a  r  ^ 

height.  However  at  low  speeds  h  =  J  (U  qdt  -  wdt)  is  dominated  by  the 
w  terra  and  cannot  effectively  be  controlled  by  cyclic  pitch.  The 
h/0Q  transfer  function  is  thus  of  prime  importance, at  low  speeds. 

H>ight  Error  Rsapon«e  to  Collective  Input 

Figure  4o  shows  the  time  histories  of  the  responses  to  a  unit 
collective  pitch  input  for  the  S-58  at  20  knots  (l)  in  level  flight 
and  (2)  in  a  steep  descent  (yQ  =  -  4l  .8  deg).  Note  that,  in  the 
descending  case, the  height  error  goes  the  "wrong  way"  after  about 
8  seconds,  never  returning  to  the  "right"  direction.  By  contrast, 
for  level  flight  during  the  first  50  seconds  the  height  error  moves 
only  a  little  way  in  the  wrong  direction,  between  IT  and  22  seconds  after 
the  step  input. 

Both  responses  are  oscillatory.  This  is  due  to  the  unstable 
phugoidmode.  As  shown  later,  this  oscillatory  mode  can  be  stabilized 
by  standard  feedback  of  pitch  rate  to  longitudinal  cyclic.  Considering 
the  remaining  modes  in  Figure  40,  there  is  a  fundamental 
similarity  be twe^n  the  responses,  in  that  both  ultimately  go  in  the 
wrong  direction.  This  is  characteristic  of  the  step  response  of  a 
system  with  one  nonminiraum  phase  zero.  It  can  be  proved  (from  the  Laplace 
transform  initial  and  final  value  theorems)  that  for  such  a  system, 


* 

This  is  shown  by  the  partial  fractions  listed  on  page  145. 


TABLE  XX.  SIKORSKY  S-58  NUMERATORS  FOR  COLLECTIVE  CONTROL  OF  NORMAL  DISTANCE  FROM  UN 
FLIGHT  PATH  h/,  9  • 
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Figure  Uo.  Sikorsky  S-58:  20  Knots,  S.L.  Time  History  of  Height  Error  (Perpendicular 
Distance  From  Unperturbed  Flight  Path)  Due  to  +  0.1  Radian  Step  Input 
Increase  in  Collective  Pitch. 


the  sign  of  the  initial  slope  of  the  step  response  is  different  from 
the  sign  of  the  final  value  of  the  response. 


The  difference  between  the  responses  of  Figure  40  lies  in  the 
rapidity  with  which  they  move  in  the  wrong  direction,  and  it  is  this 
factor  which  causes  difficulty  of  control.  To  demonstrate  this,  form 
the  -unit  step  response  time  history  by  expanding  the  Laplace  transform 
of  i/s  times  the  transfer  function,  and  taking  the  inverse  transform 
of  each  partial  fraction.  This  yields,  per  radian  of  increased  collective, 

Level  Flight: 

h  =  -  670  +  377. 2t  -13. 6t2  +  869e"0,585t  +465e“°*197t 
-673e0,06t  (cos  0.33t  -0.28  sin  0.33t) 

+0.033e“^*76t  +21.7e"1,1^t  -28e“0ji|6t(cos  0.73t  -2.4  sin  0.73t) 
Steep  Descent 

h  =  -  1,298  +7. it  -9.94t2  +568e“0,385t  +1 ,496e“°* l8t 
-8l2e*^^  (cos  0o4t  -1.5  sin  0.34t) 

-0.031  e“^,98t  -6.9e“1,07t  +53e"0, 196t(cos  o.65t  +1.72  sino.65t) 

Figures  4l  and  42  show  how  the  above  components  sum  to  form  the  complete 
response.  The  n/0o  transfer  function  has  two  poles  at  the  origin  (see 
Figure  43,  or  Tables  XX, and  XV);  hence,  the  step  input  yields  partial 
fractions  involving  i/s'5,  i/s  ,  and  i/s,  corresponding  to  parabolic,  ramp, 
and  constant  components  of  the  step  response  time  history.  These 
components  are  primarily  responsible  for  the  differences  between  the 
responses  for  level  flight  and  descent.  The  components  depend  on  the 
vector  distances  from  the  poles  at  the  origin  to  other  uncancelled  poles 
and  zeros,  and  are  markedly  affected  by  the  change  in  the  right  half  plane 
zero  from  s  =  0.074  for  level  flight  to  s  =  0.268  for  steep  descent. 

Closed- Loop  Control  of  H>ight  Prior 

The  adverse  effect  of  the  shift  in  the  nonminimum  phase  zero  from 
0.078  rad/sec  for  level  flight  to  0.268  rad/sec  for  70=  -  4l .8  degrees 
can  be  demonstrated  by  constructing  a  root  locus.  The  simplest  feed¬ 
back  to  consider  is  a  "pure  gain"  feedback  of  fi  to  collective.  This  is 
illustrated  on  Figure  43,  which  shows  the  7  =  -  4l  .8-degree  case. 

Note  that  even  at  low  gains,  a  root  moves  close  to  the  zero  at  s  =  0.268. 

For  example,  for  a  feedback  gain,  G./n,  0.08  radian  per  100-foot 

height  error,  the  system  is  unstable  with  a  root  at  8  =  0.13*  i.e., 
doubling  amplitude  every  3.3  seconds.  This  divergence  is  too  rapid 
to  be  acceptable. 

For  level  flight*  the  z^ro  is  still  present  but  is  reduced  in 
magnitude  from  s  -  0.268  sec"  to  s  -0.078  sec”1.  For  the  same  gain 
as  used  in  descent  (0  /n  *  0.08/100  radian  per  foot);  the  closed- locp 
time  constant  is  appr8ximately  0.037;  hence  the  time  to  double 
amplitude  is  increased  to  12  seconds. 
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.  Sikorsky  3-t>8,  Level  Flight,  20  Knots,  S.  L. 
Components  of  Height  Error  Time  History  per  1 .0 
Radian  Step  Increase  in  Collective  Pitch. 
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Figure  42.  Sikorsky  S-58:  20  Knots,  S.L.  Steep  Descent 
V  =  -  4 1.8  Degrees  Components  of  Height  Error 
Time  History  per  1.0  Radian  Step  Increase  in 
Collective  Pitch. 

145 


M 

o 


146 


i 


i 


Figure  1+3.  Sikorsky  S-58:  Height  Error/Collective  Root  Locus 

U  =  20  Knots  y=  -  4l .8  Deg. 


Of  course,  the  "  pure  gain"  closure  discussed  above  is  the 
simplest  form  of  control,  and  the  instability  can  be  cured  by  more 
complicated  compensation,  plus  "wrong-way"  gain  (i.e.,  increasing 
collective  to  compensate  for  being  above  the  desired  flight  path) . 
However,  little  can  be  done  to  alter  the  position  of  the  nonminimura 
phase  zero,  which  (as  shown  below)  is  the  basic  cause  of  the  unsatis¬ 
factory  tracking  capability. 

The  reason  why  the  zero  is  so  tenacious  is  that,  at  low  speeds, 
collective  pitch  is  the  only  effective  way  of  controlling  height  error. 
This  can  be  seen  by  going  to  the  limit,  i.e.,  hover,  when  cyclic  pitch 
produces  essentially  no  Z- force.  Thus  there  is  no  real  possibility 
for  multiloop  control  of  height  error,  (in  passing,  note  that  this 
possibility  does  exist  on  the  AH-5&A  discussed  in  the  next  chapter 
because  the  propeller  pitch  can  effectively  control  the  flight  path.) 
Whatever  feedbacks  are  employed,  the  overall  transfer  function  relating 
height  error  to  collective  retains  the  nonminimum  phase  zero. 

Summarizing  the  above  discussions:  the  h/0  right-half  plane 
zero  is  objectionable  because: 

( 1 )  It  causes  the  initial  step  response  to  be  in  the 
opposite  direction  to  the  final  response. 

(2)  It  rn^kes  it  difficult  to  obtain  a  stable  closure  of 

the  h-*9  loop, 
o 

(3)  It  cannot  be  moved  appreciably  by  closing  other  loops 

(e.g h-^B.  )  unless  excessively  high  gains  are  used. 

I  s 

Granted  all  the  above  points,  some  further  explanation  is  still 
required.  This  is  because  right-half  plane  zeros  are  quite  common  in 
aircraft  transfer  functions  but  do  not  normally  cause  a  major  deteriora¬ 
tion  in  handling  qualities.  For  example,  the  he ight-to-e levator 
transfer  function  of  a  conventional  aircraft  flying  above  its  minimum 
drag  speed  contains  a  "far-out"  right-half  plane  zero.  However,  neither 
"far-out"  or  "very  close-in"  zeros  cause  a  significant  deterioration 
in  the  step  response.  What,  then,  is  so  bad  about  the  S-58  h/0Q 

transfer  function  at  20  knots,  y  =  41.8  degrees,  and  the  associated 
step  response  shown  on  Figure  40? 

The  answer  can  be  found  rigorously  by  considering  the  effects  of 
the  relative  magnitude  of  the  nonminimum  phase  zero  and  the  gust  break 
frequency  on  the  performance  of  an  optimal  regulator.  This  is  done  in 
Reference  39»  but  the  essential  point  can  be  grasped  by  the  following 
reasoning . 
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Consider  the  task  of  a  pilot  flying  a  helicopter  in  a  gusty 
environment,  approximated  by  a  series  of  randomly  spaced  square  wave 
w-gusts  of  random  amplitude.  Let  the  average  time  between  successive 
gusts  be  Tg  seconds.  Suppose,  the  helicopter  has  a  stable  step  response 
to  collective  which  is  "up"  for  sec  and  "down"  thereafter  (i.e.,  a 
nonminimum  phase  characteristic).  nIf  t  is  very  much  less  than  t  ,  the 
optimal  control  requires  negative  feedblck,  i.e.,  reduced  collective 
pitch,  to  counter  an  up-gust.  On  the  other  hand,  if  the  gusts  are 
very  widely  spread,  the  integrated  absolute  value  of  the  perturbation 
from  the  desired  flight  path  is  minimized  by  positive  feedback.  Between 
these  two  cases  a  condition  exists  where  the  optimal  policy  is  zero 
feedback.  At  this  condition,  there  is  nothing  that  the  pilot  (or  S.A.S.) 
can  do  to  improve  the  accuracy  with  which  the  helicopter  is  following  the 
desired  flight  path.  In  such  circumstances  the  pilot  will  regard  the 
helicopter’s  handling  qualities  as  unsatisfactory. 

It  can  be  shown  (Reference  39)  that  with  a  single  nonminimum 
phase  zero  located  at  s  =  n  rad/sec,  and  standard  gust  power  spectrum 
of  the  form  0  =  k/(s2  +  cu2),  the  above  condition  occurs  when 

ww  u 

n  =  o^,  i.e.,  when  the  zero  is  of  equal  magnitude  to  the  input  break 

frequency.  For  hover  and  low-speed  flight,  is  normally  assumed  to 

be  =  0.3  rad/sec  (see  Appendix  VIII).  Hence,  the  zero  located  at 
s  =  +  0.268,  as  predicted  for  the  yQ  --  4 1.8-degree  case,  will  adversely 
affect  handling  qualities. 

Strictly,  the  above  discussion  applies  only  to  single-loop 
control.  However,  at  low  speeds,  control  of  ft  is  virtually  limited 
to  collective  pitch,  so  this  restriction  is  realistic.  A  zero  cannot 
be  moved  provided  only  one  control  is  used.  Hence, no  cure  is  possible 
through  stability  augmentation  in  the  collective  loop. 

Iffect  of  Stability  Augmentation 

For  simplicity,  the  foregoing  analysis  has  considered  only  the 
unaugmented  helicopter.  The  effect  of  stability  augmentation  on  the 
nonminimum  phase  zero  is  not  large,  as  explained  above,*  however,  some 
benefit  ij-s  obtained,  since  there  is  a  limited  effect  of  cyclic 
pitch  on  h.  This  is  illustrated  by  Table  XXI,  which  shows  the  effect 
of  a  typical  S.A.S.  feedback  on  the  poles  and  zeros.  With  a  feedback 
gain  just  adequate  to  stabilize  the  phugoid,  the  h/9o  zero  remains  about 
the  same  magnitude  as  in  the  unaugmented  case.  The  time  histories 
(see  Figure  44)  show  similar  objectionable  "wrong-way"  characteristics 
to  the  unaugmented  case  (Figure  4o)j  the  stability  augmentation  merely 
damps  the  phugoid. 
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TABLE  XXI.  EIGENVALUES  AM)  HEIGHT  ERROR  NUMERATORS  FOR  THE  S-*8  WITH  q—  B  STABILITY  AUGMENTATION 
GIVING  TOTAL  Mq.  EQUAL  TO  POUR  TIMES  THAT  OF  THE  UNAUGMENTED  AIRCRAFT. 
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Figure  Mu  Sikorsky  S-58:  20  Knots  Airspeed,  S.L.  Stability  Augmented  by 

Pitch  Rate  Feedback  to  B1  ,  Giving  Total  Pitch  Damping  Equal  to 
Four  Times  That  of  the  Unaugmented  Airframe.  Response  to  0.1 
Rad  Step  Increase  in  Collective  Pitch. 


At  higher  speeds,  e.g.,  ^0  knots,  cyclic  pitch  has  more 
effect  on  the  flight  path,  so  q  — ►  B.  feedback  is  effective  in 

reducing  the  wrong- way  characteristic. 

Sunatry 

For  single- rotor  helicopters,  as  typified  by  the  S-58,  the 

effects  of  longitudinal- lateral  cross-coupling  on  the  eigenvalues 

are  relatively  small.  However  the  cross- coupling  is  large,  as 

shown  by  the  eigenvectors  and  such  transfer  functions  as  p/B-  and 

I  s 

q/A1g.  Except  at  very  low  speeds,  the  variation  of  the  eigenvalues 

with  rate  of  descent  is  small.  The  cross- coupling  also  changes 
little  between  level  flight  and  the  upper  boundary  of  the  vortex¬ 
ring  state. 

The  major  reason  why  control  is  degraded  in  steep  descents 
is  shown  by  the  transfer  function  relating  collective  pitch  to 
the  distance  of  the  helicopter  from  the  unperturbed  flight  path. 
This  transfer  function  contains  a  right-half  plane  zero  which 
cannot  easily  be  removed  by  multiloop  feedbacks.  The  effect  of 
the  zero  is  to  cause  the  step  response  to  collective  pitch  to 
be  largely  in  the  wrong  direction  when  performing  steep  descents 
at  low  speeds. 

Because  of  the  large  number  of  derivatives  involved  in 
the  coupled  longitudinal-plus- lateral  equations  of  motion,  it  has 
not  yet  been  possible  to  identify  the  significant  derivatives 
causing  the  nonminiraum  phase  zero.  The  question  arises  as  to 
whether  this  is  some  special  peculiarity  of  the  S-58  or  whether 
it  is  indeed  representative  of  single-rotor  helicopters  in  general. 
Conclusive  proof  will  require  a  more  general  analysis;  however, 
in  the  next  chapter  it  is  shown  that  the  same  nonminimum  phase 
effect  occurs  on  the  AH-56A.  It  is  therefore  probable  that  the 
effect  is  typical  of  single- rotor  helicopters. 
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CHAPTER  VIII.  COMPOUND  HELICOPTER  DYNAMICS 


VHI.1 


This  section  discusses  the  open-loop  and  closed-loop  dynamics  of 
a  typical  compound  helicopter.  The  principal  features  which  differen¬ 
tiate  such  vehicles  from  "single-rotor"  helicopters  are  the  wing  and 
the  auxiliary  thrusting  device,  which  may  be  a  jet  engine  or  a  propeller. 
Both  these  features  are  intended  to  unload  the  rotor  at  high  speeds . 

One  might  expect,  therefore,  that  at  low  speeds,  there  would  be  few 
significant  differences  between  the  dynamics  of  the  "  single- rotor "  and 
"compound"  configurations.  This  proves  to  be  the  case;  the  eigen¬ 
values  and  eigenvectors  display  similar  variations  with  forward  speed 
and  descent  angle  for  both  types  of  helicopters.  For  this  reason 
frequent  reference  will  be  made  to  the  preceding  discussion  of  single¬ 
rotor  helicopter  dynamics,  and  Chapter  VII  should  be  read  bef  re  this 
chapter. 

The  example  vehicle  chosen  to  represent  compound  helicopters 
is  the  Lockheed  AH-56A.  As  originally  flown,  this  helicopter  incor¬ 
porated  an  all-mechanical  stability  augmenter  system  (the  "gyro"). 
Difficulties  were  experienced  in  flight  as  a  result  of  coupling 
between  the  gyro  and  blade  modes  other  than  the  first  flapping  mode. 

At  the  time  of  writing,  alternative  S.A.S.  schemes  are  under  con¬ 
sideration,  and  for  our  purposes,  it  appears  logical  to  consider 
the  basic  AH-56A  with  gyro  removed.  The  resulting  configuration 
should  be  more  representative  of  the  dynamics  of  compound  helicopters 
in  general  than  the  AH-56A  with  gyro,  and  the  data  presented  here 
are  thus  suitable  for  preliminaiy  design  of  S.A.S.  systems  and  for 
prediction  of  handling  qualities. 

The  detailed  AH-56A  data  presented  in  Volume  IV  is  similar 
to  that  previously  discussed  for  the  S-58.  That  is,  it  includes 
derivatives  and  transfer  functions  relating  u,  v,  w,  p,  q,  r  to  Alg, 

Bls,  e0TR>  e0p  (coHective  pitch  of  the  pusher  propeller).  Residues 

and  eigenvectors  were  also  calculated  and  are  quoted  in  this  chapter 
where  appropriate,  but  they  are  not  listed  in  Appendix  VII . 

Hie  above  information  is  presented  for  the  AH-56A  at  sea  level, 
wHh  true  airspeeds  of  0,  20,  4o,  60,  and  100  knots.  For  each  speed, 
descent  rates  of  0,  9.6,  19.2,  28. 8  fps  are  considered.  At  low  speeds, 
the  highest  descent  rate  is  almost  at  the  boundary  of  the  vortex-ring 
state . 
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Two  cases  were  re-run  with  the  cyclically  varying  tip  loss 
factor  set  to  zero.  This  caused  noticeable  shifts  in  the  eigenvalues 
but  no  radical  change  in  the  general  nature  of  the  vehicle  response. 
That  is,  inodes  which  were  stable  with  tip  losses  included  remained 
stable  with  B=1,  and  the  residues  were  of  the  same  order  of  magnitude 
with  and  without  tip  losses.  For  example,  at  Uo  knots  and  28.8  fps 
descent  rate  (7  =  -25.1  deg),  the  eigenvalues  changed  as  indicated 
below. 


AH-56A  Uo  Knots,  8te«p  Descent,  7  •  -25. 1  Deg 


TIP  LOSS 

INCLUDED 

OMITTED 

-.255  +  .758J, 
-234  +  .763. 

-9.9, 

-8.36. 

-.063, 

-.061*. 

-.526, 

-so?. 

RSUH 

The  above  changes  are  rather  less  than  might  be  expected 
by  considering  the  effect  of  tip  losses  on  the  derivatives.  The 
explanation  is  that  the  magnitudes  of  some  derivatives  (e.g.,  Zy) 
are  decreased  by  the  tip  loss,  whereas  others  (e.g.,  Mn,  Lp)  are 
increased  in  magnitude  because  tip  loss  gives  increased  flapping. 

In  all  the  cases  discussed  below,  tip  losses  were  included. 

VIII. 2  EIGENVALUES 

Table  XXII  summarizes  the  AH-56A  eigenvalues.  In  le/el 
flight,  the  eigenvalues  are  generally  similar  to  those  calculated 
for  the  S-58,  at  the  corresponding  speeds.  The  major  differences 
are 

(1)  The  AH-56A  roll  subsidence  root  is  about  twice  as  large 
as  that  of  the  S-58,  due  to  the  increase  in  "effective 
hinge  offset"  of  the  hinge less  rotor. 

(2)  The  AH-56A  develops  longitudinal  static  stability  at 
lower  speeds  than  the  S-58.  This  is  due  to  the  negative 
My  contributions  of  the  tail  and  the  pusher  propeller. 
The  result  is  that  the  phugoid  is  stable  at  speeds 
above  50  knots,  approximately.  In  addition,  the  short- 
period  degenerate  modes  combine  to  form  an  oscillatory 
mode  at  100  knots. 

In  descent,  the  eigenvalues  generally  do  not  change  signif¬ 
icantly  from  their  values  in  level  flight  at  the  same  airspeed.  The 
exceptions  to  this  statement  are  noted  below. 

At  the  higher  speeds,  descent  increases  the  phugoid  damping 
and  causes  some  reduction  in  frequency.  This  would  be  expected  from 
the  increase  in  jXul*  At  100  knots  and  the  most  rapid  descent  rate 
(which  corresponds  to  7  =  -  9.8  deg)  this  trend  progreases  so  far 
that  the  phugoid  degenerates  into  a  pair  of  aperiodic  roots.  One 
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TABUS  mi.  LOCUSED  AH- 56A  KICBtVAUJES 


15U 


of  these  roots  is  slightly  unstable  at  s  =  0.073*  This  "tuck" 
node  instability  is  caused  by  the  change  in  ^  that  occurs  between 
60  and  100  knots  at  28.8  fps  descent  rate.  (At  6o  knots, ^=+0.336  x  1CPH 
rad/sec2;  at  100  knots, Mq—17‘6  x  10*1*  rad/sec2.) 

Apart  fron  this  phugoid  change,  the  only  other  significant 
effect  of  descent  on  the  eigenvalue*  occurs  at  vertical  and  near¬ 
vertical  descents.  For  such  conditions,  two  nodes, which  at  hover 
are  aperiodic,  coalesce  to  for*  a  well-danped  low-frequency  oscil¬ 
lation.  The  eigenvectors  (discussed  later)  indicate  that  this 
behavior  is  due  to  This  derivative,  referred  to  body  axes, 
increases  by  a  factor  of  approximately  30  going  fro*  hover  to  the 
wax! nun  vertical  descent  rate  at  the  boundary  of  the  vortex-ring 
state.  This  narked  change  occurs  because  the  tail  rotor  thrust 
becones  no  re  sensitive  to  perturbations  in  vertical  velocity  as  the 
rate  of  descent  is  increased.  A  similar  effect  was  noted  on  the 
S-58. 


VIII. 3  CEOSS-COUPLIHG  EFIXCTS 


As  with  the  single- rotor  helicopter,  lateral- longitudinal 
cross-coupling  is  severe.  This  is  illustrated  by  Table  XXIII  which 
presents  AH-5&A  eigenvectors  at  hover.  Note  that  u,  v,  w,  p,  q,  r, 
are  referred  to  stability  axes  which,  in  hover,  are  arbitrarily  defined 
such  that  the  x-axis  is  horizontal  and  the  y-axis  is  nomal  to  the 
plcne  of  the  symmetry  of  the  fuselage,  which  is  almost  but  not 
exactly  vertical.  (The  print-out  gives  the  stability  axis  Euler 
angles  in  radians  for  each  flight  condition  • )  In  Table  XXIII  the 
choice  of  the  unit  eigenvector  varies  arbitrarily  from  mode  to  mode. 
As  in  Table  XVII,  only  the  upper-half  plane  root  of  each  complex  pair 
ia  shown. 


Examination  of  Table  XXIII  reveals  the  dominant  degrees 
of  freedom  of  each  mode,  and  the  names  of  each  mode  are  chosen 
accordingly.  However,  with  the  exception  of  the  roll  subsidence, 
each  mode  displays  strong  lateral-longitudinal  coupling,  so  the 
names  should  not  be  taken  too  literally.  The  resemblance  of  each 
mode  to  the  corresponding  uncoupled  mode  is  not  necessarily  close. 
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(a,  v,  v.  p.  i,  r,  in  stability  axes) 


-.235 

-.205 

.035 
♦  .M.5J 

-.018 

♦•59*0 

-10.0 

-2.01 


.014 

.829 


1.0 


loop 

lOOq 

lOOr 

Mode  Name 

.00067 

-.039 

.569 

”  Plunge” 

.0058 

-.102 

1.85 

”  Yaw  ” 

.015 
-  •  175J 

.518 

-.117J 

.231 

♦.4l7j 

"Fhugoid” 

-.48 

-.394 

.681 

”  Dutch 
Roll” 

♦.009J 

-.0967J 

-1.27J 

100.0 

6.3 

-1.95 

”RoU 

Subsl 

-dence” 

-2.02 

-10.56 

-24.2 

” Speed 

-Pitch” 

TABLE  XXIV.  AK-56A  EIGENVECTORS  IN  28.8 


-.181 

•  391 

+  .09U 

--348J 

-.039 

♦.439J 

-1.0 

-.08 

-.449 

♦.362,} 

•  229J 

-9.8 

-.079 

-1.91 

-1.0 

.012 
.171  .0045 


lOOq 

lOOr 

.00135 
+  .0076J 

1.32 

-4.38J 

-.57 
- . 119 J 

.237 

♦.l66j 

-.20k 

-.O09J 

-.717 

♦1.48J 

-2.85 

-.824 

-21.8 

♦.77 

lOOq 

-lOOp 

e  Naae 
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The  cross-coupling  hardly  changes  with  descent  rate,  provided 
that  u,  ▼,  w,  p,  q,  r  are  referred  to  body  axes  rather  than  stability 
axes.  This  point  is  important  for  very  steep  descents.  For  example, 
in  vertical  descent  the  body  x-axis  and  the  stability  x-axis  differ 
by  90  degrees.  This  must  be  considered  when  comparing  Table  XXIII  with 
Table XXIV,  its  counterpart  for  28.8  fps  vertical  descent. 

At  the  foot  of  Table  XXIV  are  listed  the  stability  axis  quantities 
corresponding  to  the  body  axis  quantities  at  the  head  of  the  table. 

Thus,  for  example,  to  trace  the  effect  of  vertical  descent  on  the 
**dutch  roll"  mode,  note  that, at  hover,  the  ratio  of  |u|:  |v|:  |w|  is 
approximately  0.7?:  1.0: 0.045.  For  vertical  descent  the  corresponding 
ratio  in  body  axes  is  0.5:  1*0:  0.10.  Qy  comparing  the  eigenvectors 
in  this  fashion,  it  is  seen  that  there  is  not  much  difference  between 
the  hover  and  descent  modes  with  the  exception  of  the  plunge  and  yaw 
modes  which  couple  together  in  vertical  descent.  This  is  caused  by 
the  Increased  sensitivity  of  the  tail  rotor  to  perturbations  in 
vertical  descent  velocity.  The  tail  rotor  is  not  immersed  in  the 
main  rotor  wash  at  hover  and  in  vertical  descent.  At  hover,  by 
symmetry,  the  tall  rotor  thrust  is  equally  sensitive  to  "up"  and 
"down"  w  perturbations;  hence,  its  contribution  to  the  cross-coupling 
derivative  1^  is  negligible.  However,  as  the  descent  rate  Increases, 

the  unperturbed  flow  conditions  at  the  tall  rotor  become  asymmetric, 
and  the  derivative  Increases  markedly,  causing  the  "yaw"  and  the 

"plunge"  modes  to  couple. 

Cross-coupling  is  pronounced  at  forward  speed.  For  example, 
the  short-period  oscillatory  mode,  s-  -1.58  *  0.556J,  which  emerges 
at  100  knots  Involves  about  one- third  as  much  roll  as  pitch.  This 
is  Illustrated  by  Table  XXV,  which  presents  some  typical  residues  for 
pitch  and  roll  response  to  B^g.  Table  XXV  shows  that  the  cross- coupling 
does  not  change  much  with  shallow  angles  of  descent. 

The  pitch  rate/B.  transfer  function  also  changes  rather 
little  from  level  flight  to  steep  descent,  as  shown  by  Figure  45. 
Numerous  near-cancellations  occur,  leaving  essentially  a  second-order 
numerator  and  a  third-order  denosdnator. 
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TABUS  XXV.  EFFECT  OF  DESCENT  OH  CROSS  -  COUPLING  j  LOCKHEED  AH-‘6A 

AT  100  KNOTS.  SEA  1EVEL.  NO  S .A.S. 

|  Roots 

-q  Residues  to  Bl6 

-p  Residues 

i 

o 

y  •  -9.8 

y  •  o 

V-  -9.8* 

y  «  0 

y  -  -9.8° 

-.1*01 

♦1.5U 

-.1*31 

♦1.49J 

.035  . 

<  -152.8 

.034  . 
<-128.9 

.429  . 
<-137.5 

•36  . 

<  -173-0 

-10.6 

-9*56 

j  1.90 

1.526 

36.4 

3*.  95 

-.037 

-.025 

.104 

-.0079 

.0777 

.00317 

-1.58 

♦-556J 

-1.50 

♦•58j 

13-7  . 

<  -130.6 

11.7 . 

<  -135.8 

5-5*. 
c  -105.0 

• 

• 

V 

-.0**1 
♦  .l88j 

-.129 

•56  . 

<  -89.6 

-517 

.537  . 

<  +98.6 

.265 

.0733 

-.287 

-.32* 

As  noted  previously,  the  AH-*>6A  hss  s  large  roll  subsidence 
root,  (typically  about  10  rad/sec),  due  to  the  hingeless  rotor.  This 
is  well  separated  Aron  the  other  roots,  which  should  lead  to  little 
coupling  with  other  nodes  in  the  p/A-  transfer  function.  This  is 
confirmed  by  Figure  1*6,  which  shows  tnf  transfer  functions  for  p/Aja 
at  100  knots  in  level  flight  and  22.8  fps  descent.  Note  that  only 
the  roll  subsidence  pble  is  not  cancelled  by  a  nearby  zero. 

The  above  trends  are  typical.  For  ease  of  reference, 

Table  XXVI  has  been  extracted  fron  the  data  of  Appendix  VII  to  show 
the  effect  of  descent  angle  on  the  #/Aja  and  0/Bja  numerators  at 
various  forward  speeds.  Fron  this  table  and  Table  XXV,  it 
is  evident  that  the  effect  of  descent  angle  on  pitch  and  roll  control 
is  small.  Similar  comments  apply  to  yaw  control  by  tail  rotor 
collective  (see  Table  XXVII).  Table  XXVII  shows  that  the  yav/0_ 
numerators  change  little  with  descent  angle,  provided  that  the 
yaw  is  referred  to  body  axes  (i.e.,  overall  vehicle  axes)  rather 
than  stability  axes. 
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AH-SS.A 

100  KNOTS  ,  LEVEL  PLISHT 

ROOT  LOCUS  SAIN  *15.3 


Figure  4^.  Compound  Helicopter  Pitch  Rate/Longidutinal  Cyclic 

Transfer  Function  for  Level  Flight  and  for  a  9.8-Degree 
Descent  Angle  at  Sea  Level. 
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a  1 

Figure  *‘6.  Compound  Helicopter  Roll  Rate/Lateral  Cyclic 

Transfer  Functions  for  Level  Flight  and  9.8-Degree 
Descent  at  Sea  Level. 
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VIII. 4  HEIGHT  CONTROL 


As  with  the  S-58  (discussed  in  the  previous  chapter),  it  is 
anticipated  that  the  AH-56X  will  experience  difficulty  in  flight  path 
control  in  slou,  steep  descents .  This  can  be  inferred  from  Table  XXVIII , 
which  presents  numerators  for  the  transfer  function  relating  height 
error  to  collective  pitch.  (See  Figure  59  for  the  definition  of 
"height  error",  n.)  Table  XXVIII  shows  that, at  the  lower  speeds  and 
steeper  descent  angles  a  nonminiraura  phase  zero  appears  at  0.1  <  s  <0.6 
rad/sec.  As  shown  in  Chapter  VII,  such  a  zero  degrades  the  handling 
qualities.  Typical  effects  of  this  zero  are  shown  ^.n  the  step  responses 
of  Figure  h'J ,  which  display  awkward  "wrong-way"  characteristics  for 
the  descending  flight  conditions. 

The  difficulties  of  removing  this  right-half  plane  zero  by 
feedbacks  to  cyclic  pitch,  explained  in  Chapter  VII  for  the  S-58,  apply 
equally  to  the  AH-56A.  However,  the  AH-56A  does  have  the  extra  possibility 
of  flight  path  control  through  the  pusher  propeller,  as  described  below. 

Control  Ifflsctlvtneis  of  tikm  Pusher  Propeller 

Assuming  that  the  aircraft  is  trimmed  so  that  the  pusher  propeller 
is  not  operating  in  the  vortex- ring  state  (a  condition  which  is  satisfied 
for  the  flight  conditions  presented  in  thip  report),  the  principal 
effects  of  0  relate  to  control  of  u  and  h.  The  appropriate  numerators 
for  level  flight  and  steep  descent  are  summarized  in  Table  XXIX. 

Completing  the  transfer  functions  with  the  denominators  given  in  Table 
XXII  permits  some  preliminary  assessment  of  the  effectiveness  of  the 
pusher  propeller  as  a  control  device. 

At  100  knots  in  level  flight,  many  approximate  cancellations 
occur  between  numerator  and  denominator,  leaving  the  following  approximate 
transfer  function  for  speed  control: 


u  .  -5^.0 

0Qp  =  (s  -  0.075) 

This  indicates  that  the  pusher  propeller  controls  speed  by 
exciting  the  unstable  "tuck"  mode  of  the  degenerate  phugoid  pair.  This 
would  be  unsatisfactory  for  x(= /udt)  control  (e.g.,  stationkeeping) 
since  the  x-*0  locus  lies  entirely  in  the  right  half- plane.  The 
situation  wouLcrbe  improved  by  stability  augmentation  involving  feed¬ 
backs  of  0  to  B.8>which  would  stabilize  both  phugoid  roots.  Descent 
improves  the  u/6®  transfer  function:  at  100  knots  and  7  =  -  9.8  degrees, 
it  becomes  approximately 
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TABLE  XXVIII.  HEIGHT  ERROR/COLLECTIVE  injMERATORS  FOR  THE  AH-,6A 
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u _  -34.1  (»  ♦  0.095) 

6  “  (a  ♦  0.04  ♦  0.1B8J) 

op  - 


This  provides  stable  x-*0q^  closures  provided  that  excessive  lags  are 
avoided. 

At  lower  speeds,  fewer  approximate  cancellations  occur.  For 
example,  at  4o  knotr  In  level  flight* 


u _  .  -18.2  (s  +  0.162  ±  0.0623) _ 

0  *  (s  -  0.072)  (u  ♦  2.07)  («  ♦  0.052  +  O.30J) 

op  - 


whereas  at  4o  knots 


in  descent  at  y 


-  25.1  degrees. 


u _  .  -18.14  (a  4  0.174)  (s  -  0.037  ±  0.312J) 

0  '  (s  +  0.526)  (s  ♦  1.09)  (•  +  0.041  ♦  0.2b2J) 

op 


The  appearance  of  a  nonminimua  phase  zero  in  descent  is  not 
important  here  because  of  the  proximity  of  the  phugoid  mode,  which 
will  yield  only  a  small  residue  and  will  largely  cancel  out  the 
wrong- wiy  characteristics  of  the  step  response. 

Simultaneously  with  u- perturbations,  3  also  induces  h 
perturbations;  hence  it  is  necessary  to  consider  these  also  in 
assessing  the  control  effectiveness  of  the  pusher  propeller.  The 
appropriate  transfer  functions  are  given  in  Table  XXIX;  some  rather 
gross  cancellations  have  been  made  in  the  approximate  transfer  functions 
presented  below  in  order  to  Amplify  the  discussion. 


1.62  (s  +  0.475)  (a  -  0.349)  (s  -  4302) 
s(s  +  I.58  +  0.356J)  (s  +  0.041  +  O.lBfcj) 


-3.13  (s  +  3.97)  (s  -  2.12) 
s(s  -  0.073)  (s  +  1.5  +  oT^BjfT 
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1_ 

0 

op 

y  =  -  9.8  degrees 
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The  change  in  the  sign  of  the  root- locus  gain  is  not  significant, 
because  of  the  change  in  the  number  of  right-half  plane  zeros.  The 
zero  at  s  »  +  0.3^9  is  likely  to  induce  a  poor  step  response  in  level 
flight,  but  this  should  be  less  troublesome  in  descent.  However,  the 
presence  of  these  zeros  will  cause  difficulty  in  closing  the  n-»() 
loop,  so  at  100  knots  the  pusher  propeller  is  not  a  desirable  p 
alternative  to  the  usual  and  n-*  0^  loop  closures. 

At  low  speeds  the  pusher  propeller  shows  better  promise  for 
control  of  height  error.  For  example,  at  knots  and  y  =  -  2*.l 
degrees,  the  n/0  transfer  function  is,  approximately. 


h  ^  -8.6 

0  3  s(s  ♦  0*>26) 

op 


This  avoids  the  awkward  nonminimum  phase  characteristics  of  the 
n/collectivc  transfer  functions  listed  on  Table  XXVIII.  Hence  the 
pusher  propeller  may  provide  a  useful  method  of  flight  path  control 
at  low  speeds. 


The  control  characteristics  of  compound  helicopters  have  >cn 
studied  by  examining  the  AH-56A  with  the  all-mechanical  stability 
au@nenter  system  (the  "gyro")  removed.  From  calcinations  of  derivatives, 
transfer  functions,  and  residues,  the  AH-96A  open- loop  behavior  is  seen 
to  be  generally  that  of  a  typical  single-rotor  helicopter.  Strong 
coupling  occurs  between  longitudinal  and  lateral  motions,  and  there  is 
not  much  effect  of  descent  angle  on  the  eigenvalues  except  at  very  low 
speed.  In  slow, steep  descents,  control  of  flight  path  with  collective 
is  difficult  because  of  the  appearance  of  nonminimum  phase  zeros  in  the 
appropriate  transfer  functions.  This  effect  is  similar  to  that 
described  for  the  single- rotor  helicopter.  However  the  compound  has  a 
possible  alternative  method  of  flight  path  control,  through  variations 
of  collective  pitch  of  the  pusher  propeller.  This  circumvents  the  non- 
minimum  phase  effects  and  appears  to  provide  a  satisfactory  method  of 
controlling  the  flight  path. 
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cam  XX.  COMCLUSIOMS  AMD  RECOHMCMDATIONS 


The  principal  results  of  this  report  are  contained  in  the  tables 
of  derivatives  and  transfer  functions  given  in  the  Appendixes.  A  large 
part  of  the  value  of  the  report  lies  in  these  detailed  nwerical 
results  and  in  the  presentation  of  the  analytic  and  computational 
methods  used  to  obtain  them,  e.g.,  the  M06TAB  program.  However,  certain 
general  conclusions  have  been  reached,  and  these  are  presented  below, 
together  with  reconae Delations  for  further  work. 

1 .  A  new  method  for  predicting  the  boundaries  of  the  vortex- ring 
state  is  presented,  using  only  simple  concepts  but  giving 
good  agreement  with  experiment. 

2.  The  tip  losses  on  a  rotor  in  vertical  descent  can  be  predicted 
by  r;onentun  and  bl»*de-clenent  theory,  using  n  flow  model  which 
incorporates  a  sheath  of  vortlcity  surrounding  the  slipstream. 
The  method  gives  good  agreement  with  experiment  in  vertical 
icscentf  no  suitable  experimental  data  are  available  for  inc 
Inclined  descent  to  check  the  theory  for  that  condition,  for 
which  it  predicts  a  cyclically  varying  tip  loss  factor. 

3*  In  calculating  helicopter  derivatives,  the  perturbations  in 
interference  velocities  induced  by  each  element  on  each  other 
element  must  be  included.  All  lateral- longitudinal  cross¬ 
coupling  derivatives  should  be  incited  on  single-rotor  and 
compound  configurations. 

it.  Derivatives  calculated  for  steep  descents  are  generally 
appreciably  different  from  derivatives  calculated  for  the 
same  airspeed  and  level  flight.  The  use  of  level- flight 
derivatives  for  simulations  of  moderate  and  steep  descents 
can  give  seriously  misleading  results. 

For  tandem- rotor  configurations,  the  principal  effects  of 
descent  are  to  cause  the  dutch  roll  to  degenerate  into  an 
aperiodic  pair  of  roots  and  to  worsen  the  nonminimum  phase 
characteristics  of  the  roll-to-lateral-cyclic  numerator. 

The  net  effect  is  to  worsen  the  open- loop  dynamics, 
increasing  the  need  for  multiloop  feedbacks.  Longitudinally 
the  effects  of  descent  angle  are  generally  minor,  except  for 
changes  in  the  speed/longitudinal  cyclic  numerator  referred 
to  cockpit  axes.  This  could  be  important  for  stability 
augnenter  systems  using  air  data  sensors. 
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•  For  a  ingle- rotor  helicopters,  lateral- longitudinal  cross- 
coupling  effects  are  important  in  the  transfer  function 
numerators  but  have  little  effect  on  the  denominators .  The 
most  significant  effect  of  descent  angle  occurs  at  low 
speeds,  where  control  of  the  flight  path  by  collective  pitch 
becomes  seriously  degraded  due  to  nonminimum  phase  effects, 
causing  "wrong-way"  responses.  This  is  believed  to  be  the 
principal  cause  of  the  difficulties  encountered  by  helicopters 
attempting  to  follow  steep  approach  paths  at  low  speeds. 
Remedies  involving  multiloop  feedback  or  elaborate  compensa¬ 
tion  of  single-loop  feedbacks  do  not  appear  feasible. 

7*  The  behavior  of  compound  helicopters  is  similar  to  single¬ 
rotor  helicopters  at  low  speeds,  where  the  wing  has  little 
influence.  The  propeller  provides  an  alternative  method  of 
flight  path  control  which  shows  promise  of  alleviating  the 
above-mentioned  difficulties. 


Despite  the  widespread  use  of  helicopters,  and  the  vast 
production  runs  of  certain  types,  research  on  helicopter  dynamics 
consumes  only  a  tiny  percentage  of  the  total  aerospace  bulge t. 

It  is  wise  to  recognise  this  situation  in  presenting  a  list  of 
topics  deserving  further  research.  For  example,  from  a  purely 
technical  viewpoint,  a  good  case  could  be  made  for  a  method  of 
calculating  derivatives  based  on  a  "  free- wake"  analysis, 
making  the  minimum  number  of  assumptions  regarding  slipstream 
geometry  and  interference  effects.  However,  the  computer  time 
required  to  use  such  a  program  renders  it  less  desirable  than  a 
simpler  alternative  uring  "assumed  wake”  representations.  The 
technical  areas  listed  below  have  been  selected  because  of  their 
practical  importance,  and  the  recommended  approaches  to  solving 
the  problems  have  been  chosen  with  a  realisation  of  the  probable 
financial  constraints. 

The  scope  of  the  items  listed  below  is  confined  to  the  area 
of  helicopter  steep  approaches,  but  some  of  the  recommendations 
are  sufficiently  general  to  apply  to  other  phases  of  helicopter 
operation  or  even  other  types  of  aircraft,  such  as  tilt-rotor 
configurations . 
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Tbs  ▼cctjfcHn 

The  vortex-ring  state  is  important  as  a  limiting  condition 
on  low-speed  descent/deceleration  capabilities  and  as  the  possible 
cause  of  accidents.  A  better  understanding  of  this  state  is 
required  to  define  it  more  precisely  and  to  indicate  methods  of 
alleviating  the  unsteady  flow.  In  addition,  some  unsteadiness  is 
known  to  occur  even  at  shallow  angles  of  descent,  and  a  knowledge 
of  this  self- induced  turbulence  Is  required  to  fully  define  the 
helicopter  turbulence  environment. 

Full-scale  flight  tests  are  recomended,  plus  more  refined 
theoretical  investigations  using  more  sophisticated  analytic 
techniques  than  those  employed  in  Chapter  II.  The  result  would 
be  a  reliable  data  base  (available  data  is  very  incomplete),  plus 
a  nonecipirlcal  explanation  of  the  causes  of  unsteady  flow  in 
steep  descents. 


Virtually  no  published  experimental  data  exist  on  the  forces 
and  moments  experienced  by  helicopter  rotors  in  vertical  and  in¬ 
clined  descent.  Model-scale  experiments  are  needed  on  Isolated 
rotors,  measuring  steady  and  unsteady  loads  and  the  Induced 
velocities  at  moderate  rates  of  descent.  The  flow  model  presented 
in  Chapter  III  contains  empirical  assumptions,  but  there 
seems  little  point  in  producing  a  more  complicated  theory 
without  experimental  data  which  are  suitable  for  correlation. 

Such  experiments  should  be  instrumented  to  measure  the  tip  losses, 
including  possible  cyclic  variations.  The  mnber  of  independent 
variables  and  the  desirability  of  minimizing  interference 
effects  virtually  dictate  that  the  tests  be  conducted  at  model 
scale  in  a  large  wind  tunnel  on  an  isolated  rotor.  Tests  on  a 
complete  configuration  would  be  disadvantageous  because  inter¬ 
ference  effects  from  other  components  would  complicate  the  results. 
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Stability  PriTmtlTM 

As  shown  in  this  report,  narked  differences  exist  between 
stability  derivatives  in  level  flight  and  descent,  causing 
important  changes  in  the  dynamics  of  the  helicopter.  It  is  there¬ 
fore  misleading  to  use  stability  derivatives  calculated  for  level 
flight  to  analyze  descent  problems.  Given  adequate  flow  models 
for  descending  rotors,  derivatives  can  be  calculated  with  ease  by 
M06TA.B  or  sini  ar  programs.  The  prime  question  is:  how  accurate 
are  these  calct lations?  Since  derivatives  describe  only  the 
lav*- frequency  behavior  of  helicopters  In  small  perturbations 
from  a  fixed  trim  condition,  accuracy  of  ±  20f>  should  suffice. 
However,  it  is  not  known  whether  even  this  accuracy  is  being 
obtained.  Published  flight  test  measurements  are  few,  and 
systematic  correlations  of  estimated  and  measured  derivatives 
in  descending  flight  are  conspicuous  by  their  absence. 

To  remedy  this  situation,  a  research  progrin  involving 
correlation  of  estimated  derivatives  with  existing  flight  test 
data  should  be  undertaken.  The  scope  should  be  broad,  including 
single- rotor,  compound,  and  tandem  con  figurations  in  level  flight 
and  descent.  The  results  should  be  used  to  refine  the  derivative 
calculation  process,  particularly  the  estimation  of  interference 
velocities,  which  at  present  is  more  of  an  art  than  a  science. 
Particular  attention  should  be  paid  to  low- speed  descending  fligit 
because  the  requirements  for  flight  path  control  on  the  approach 
ore  particularly  stringent,  requiring  the  most  accurate  derivatives 
for  analysis. 

TTlMftr  rUMttflM 

Transfer  functions  provide  an  accurate  description  of 
helicopter  dynamics  over  a  limited  frequency  and  amplitude  range. 
The  upper  1  requency  limit  is  set  by  considerations  of  blade 
modes  other  than  first  flapping.  Such  modes  are  essentially 
"nuisance"  modes,  i.e.,  they  have  no  value  for  flight  control. 

Thus  they  can  be  ignored  in  flight  control  system  preliminary 
design,  but  must  be  included  in  detail  design  of  system  hardware 
to  avoid  unwanted  coupling  effects.  The  lower  frequency  limit 
is  set  by  time- variations  of  the  nominal  operating  condition 
around  which  the  perturbations  are  performed.  This  limit  is  less 
restrictive  for  descent  than  for  level  flight,  since  the  available 
deceleration  is  less,  because  the  helicopter  is  closer  to  (d/l) 

Thus  there  is  a  broad  range  of  frequencies  for  which  transfer  max 
functions  are  useful.  The  following  recommendations  apply 
for  this  frequency  range. 
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Transfer  functions  including  lateral-longitudinal 
cross- coupling  should  be  used  as  the  basis  f*or  future  studies 
of  mjimiiii  and  automatic  control.  Because  of  computational 
difficulties,  cross-coupling  has  been  ignored  in  many 
previous  studies;  however  as  shown  here  cross- coupling  effects 
are  quite  large  and  significantly  influence  closed-loop  behavior. 
Currently  the  effects  of  cross-coupling  are  countered  on  an 
ad  hoc  basis  by  offsetting  the  A.  ,  B.  servos.  The  cross-coupling 
phenomenon  mist  be  more  clearly  defined  to  permit  a  systematic 
study  of  its  importance  to  handling  qualities  and  to  determine  the 
best  method  of  counteracting  It.  The  present  report  used 
manufacturer's  derivatives  for  the  tandem- rotor  helicopter,  which 
neglect  cross-coupling.  Although  this  appears  to  be  more  Justi¬ 
fiable  for  the  tandem  configuration  than  for  single-rotor  heli¬ 
copters,  cross-coupling  may  still  be  significant.  Flight  test 
data  should  be  examined  to  verify  this. 

The  significant  features  of  certain  transfer  functions 
(i.e.,  the  height  error  to  collective  msec  ra  tor  discussed  in 
Chapter  VII)  should  be  related  to  the  helicopter  geometry.  The 
usual  way  of  doing  this  is  through  the  use  of  literal  approxi¬ 
mate  transfer  (Unctions  (e.g.,  approximating  dutch  roll  natural 
frequency  as  <i>d  '  )•  Similar  expressions  should  be  developed 

for  the  critical  right-half  plane  soro  of  the  height  error  to 
collective  transfer  function  and  other  important  parameters.  The 
task  is  complicated  because  it  is  necessary  to  include  lateral- 
longitudinal  cross-coupling)  however,  it  is  essential  to  be  able 
to  trace  the  cause  of  significant  dypamic  characteristics  to 
desi&i  the  helicopter  and  its  stability  au&aenter  together  for 
good  flying  qualities. 

Tims -Varying  Dynamics  Re-presentations 

The  need  for  a  good  method  of  representing  the  dynamic 
interaction  of  the  control  system  with  modes  other  than  first 
flapping  is  generally  recognized.  No  special  emphasii  on  steep 
descents  seems  necessary  here,  since  development  of  a  method 
suitable  for  level  flight  should  meet  the  requirements  for  descent. 
Of  more  special  relevance  to  approach  is  the  question  of  the 
validity  of  small- perturbation  equations  about  a  time- invariant 
trim  condition. 

For  steep  descents  the  deceleration  capability  is  smaller 
than  in  level  flight  since  the  helicopter  is  operating  closer  to 
(D/L)  •  Accordingly,  the  effect  of  time  variations  in  the  trim 

condi??ons  should  be  less  important.  The  net  recorunendation  is 
that  this  area  is  less  pressing  for  steep  descents  than  for  level 
flight. 
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flying  Quilltlw 

Most  published  research  cn  helicopter  L lying  qualities  has 
concentrated  either  on  hover  or  on  level  flight.  Significant 
differences  have  been  shown  to  exist  between  open- loop  dynamics  in 
level  flight  and  descent,  and  correspondingly  large  differences  will 
occur  in  flying  qualities.  Ibis  area  requires  investigation  by 
simulation  in  conjunction  with  analytic  studies  of  the  man-ln-the- 
loop.  For  example,  flight  tests  indicate  that  handling  qualities 
are  significantly  degraded  in  slow,  steep  descents  (see  Chapter  II); 
the  cause  has  been  Identified  in  Chapters  VII  and  VIII  as  the  "wroob 
way"  height  control  characteristics.  More  work  involving  detailed 
loop  closures,  backed  up  by  simulation,  is  required  to  verify  this 
hypothesis • 

To  obtain  specific  results,  it  was  necessary  to  consider  only 
three  helicopters,  in  this  report.  These  are  believed  to  be  represen¬ 
tative  to  their  general  types,  but  no  information  has  been  generated 
on  the  effects  of  helicopter  else  per  se.  One  would  like  to  know 
how  descent  angle  affects  the  flying  qualities  of  the  Heavy  Lift 
Helicopter.  It  is  possible  that  for  very  large  helicopters,  body¬ 
bending  modes  may  be  of  sufficiently  low  frequency  tc  interact  with 
handling  qualities.  Derivatives  should  be  calculated  and  assessments 
made  of  handling  qualities  problems  in  steep  descents  for  such 
vehicles • 
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